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ABSTRACT 


is  thesis  investigates  the  effects  of  atmospheric 
refraction  on  military  ground-force  tactical  systems  which 
utilize  EM  energy  propagating  through  the  troposphere. 
Anomalous  propagation,  mirages,  and  atmospheric  turbulence 
are  the  three  major  atmospheric  refractive  conditions 
discussed.  EM  energy  from  visible  light  to  VHF  radio  is 
considered . 

The  refractive  characteristics  of  the  atmosphere  and  the 
conditions  that  lead  to  significant  refractive  effects  are 
reviewed.  The  means  of  predicting  atmospheric  refraction 
and  the  usefulness  of  such  predictions,  along  with  the 
probability  of  occurrence  of  atmospheric  refraction  is 
also  discussed. 

Specific  examples  of  atmospheric  refractive  effects  are 
provided  in  two  forms;  computer  simulations  and  qualitative 
studies.  A  scenario  is  simulated  showing  the  wide  range 
of  atmospheric  refractive  effects  a  typical  ground  force 
might  anticipate  in  a  desert  environment.  Finally,  the 
thesis  provides  conclusions  and  recommendations  for  further 
studies . 
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I.  INTRODUCTION 


In  order  to  enhance  and  extend  surveillance  and 
communications  capabilities,  modern  ground  warfare  has 
become  increasingly  dependent  on  the  entire  electromagnetic 
(EM)  spectrum.  Uses  of  electromagnetic  energy  include 
airborne  and  ground  target  detection,  weapons  system 
tracking  and  guidance,  command  and  control,  and  electronic 
warfare.  Utilization  of  this  energy  requires  an  under¬ 
standing  of  the  phenomena  that  affect  the  propagation  of 
EM  waves  as  they  travel  through  the  atmosphere.  After 
nearly  50  years  of  study,  it  is  now  widely  accepted  that 
dispersion,  extinction,  refraction,  reflection,  and 
diffraction  all  affect  EM  propagation  significantly. 
Dispersion,  extinction,  and  refraction  are  properties  of 
the  atmosphere,  while  reflection  and  diffraction  are 
properties  both  of  the  atmosphere  and  of  the  terrain  of 
the  earth’s  surface.  For  radio/radar  frequency  (RF)  and 
electro-optical  (EO)  systems,  the  physical  factors  listed 
impact  on  system  performance  above  and  beyond  the  more 
obvious  effects  of  precipitation,  wind,  and  temperature 
extremes . 
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The  primary  region  of  Army  interest  is  the  troposphere- 
the  lower  10  to  20  kilometers  of  the  atmosphere.^  Within 
the  troposphere,  the  main  factor  affecting  RF  propagation  is 
refraction.  EO  propagation,  while  greatly  affected  by 
aerosols  and  absorption  in  the  troposphere,  is  also  affected 
by  atmospheric  refraction.  Thus,  in  order  to  best  employ 
his  assets,  it  is  important  that  the  ground  forces  comma  'er 
be  able  to  forecast  and  predict  refractive  conditions  tl- 
degrade  or  enhance  RF  and  EO  system  performance. 

It  is  the  purpose  of  this  thesis  to  examine  the  effe 
of  atmospheric  refraction  on  ground  forces’  tactical  systems 
Specifically,  it  will  address  the  environmental  conditions 
leading  to  the  three  major  forms  of  atmospheric  refraction 
(anomalous  propagation,  mirage  formation,  and  atmospheric 
turbulence),  and  discuss  their  effects  on  EM  waves 
propagating  through  the  atmosphere.  It  will  also  address 
prediction  and  probability  of  occurrence  of  these  three 
forms  of  refraction  and  it  will  end  with  a  discussion  of 
their  effects  on  existing  ground  systems,  (either  through 
simulation  or  by  presenting  empirical  test  results). 
Recommendations  and  direction  for  future  research  directions 
will  be  found  in  the  conclusion  of  the  thesis. 

^  The  specific  height  of  the  troposphere  is  a  function  o 
location  and  atmospheric  conditions,  on  the  average  it  is 
the  thickest  at  the  equator  and  narrowest  at  the  poles. 


II.  ATMOSPHERIC  REFRACTIVE  EFFECTS 


In  this  section  the  basic  theory  of  EM  propagation  will 
be  introduced.  Since  the  purpose  of  this  thesis  is  to 
describe  the  effects  of  atmospheric  refraction  on  ground 
forces,  this  discussion  will  be  limited  to  the  troposphere. 

A.  EM  WAVE  PROPAGATION  IN  THE  ATMOSPHERE 


The  electromagnetic  wavelengths  of  interest  to  U.  S. 
ground  forces  extend  from  those  only  micrometers  (lO”  meters) 
in  length  to  those  several  tens  of  meters  long.  This 
encompasses  the  full  spectrum,  from  HF  to  v'sible  -  an 

vST 

EM  region  which  is  described  in  Figure  1, 

Electromagnetic  waves  travel  in  a  straight  line  in 
free  space  (vacuum),  but  in  the  atmosphere  they  are  bent  by 
several  forces.  The  specific  atmospheric  effect  and  layer 
of  the  atmosphere  that  dominate  in  causing  this  bending 
depend  on  the  propagating  frequency.  At  wavelengths  longer 
than  very  high  frequency  (VHF),  approximately  30  MHz,  the 
ionosphere  is  the  primary  influence.  Since  this  thesis  is 
concerned  with  the  lower  atmosphere,  HF  frequencies  will, 
therefore,  not  be  specifically  addressed.  Above  30  MHz, 
refraction  in  the  troposphere  is  the  controlling  influence, 

^  and  it  will  be  these  frequencies,  from  VHF  to  visible,  that 

this  thesis  will  address. 
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Figure  1.  Electromagnetic  Spectrum 


There  is  also  a  relationship  between  regions  of  the 


EM  spectrum  and  the  forms  of  r-  ‘^■^action  dominant  in  each 
region.  The  dominant  effects  at  microwave  frequencies  are 
long  range  changes  in  refraction.  At  optical  frequencies, 
absorption  and  small-scale  refractive  changes,  in  the  form 
of  mirage  formations  and  atmospheric  turbulence,  are  the 
dominant  effects.  In  the  millimeter  wave  (MMW)  region  the 
atmosphere  is  largely  opaque  due  to  molecular  absorption  by 
carbon  dioxide  and  water  molecules,  except  for  a  few  regions 
(or  "windows”)  of  higher  transmittance.  Although  little 
used  until  recently,  much  current  research  is  focused  on 
expanding  the  use  of  this  region  and  in  understanding  the 
effects  of  atmospheric  refraction  on  MMW  propagation. 

2 ,  Index  of  Refraction  (n) 

Electromagnetic  waves  travel  through  matter  via 
absorption  and  re-emission  of  electromagnetic  energy  by  the 
atomic  and  molecular  constituents  of  the  medium.  Thus,  the 
speed  at  which  these  waves  travel  is  affected  by  the  density 
of  the  medium  -  the  denser  the  medium,  the  slower  the  speed 
of  EM  propagation  through  it. 

The  "index  of  refraction  (n)"  is  the  ratio  of  the 
speed  of  electromagnetic  waves  in  free  space  to  their  speed 
in  the  medium.  Thus, 


where  c  is  the  speed  in  vacuum  (speed  of  light)  and  v  is  the 
speed  in  the  medium.  We  note  that  v  is  dependent  on  the 
constituents  of  the  medium  -  in  this  case,  the  atmosphere. 
Note  also  that  n  =  1  in  vacuum,  and,  therefore,  n  is  greater 
than  or  equal  to  1  for  all  media.  For  air, 

"air  = 

In  the  troposphere,  n^^^^  is  a  function  of  temperature  (T), 
water  vapor  content  (e),  and  pressure  (P),  all  of  which  are 
related  to  the  density  of  the  atmosphere  and  all  of  which 
vary  systematically  with  climate  and  elevation.  In  general, 
the  index  increases  with  increasing  density. 

3 .  Snell *s  Law 

The  deflection  of  an  electromagnetic  wave  across  a 
discontinuity  or  boundary  separating  regions  of  differing 
values  of  refractive  index  is  shown  in  Figure  2,  and  is 
described  by  Snell’s  law: 

Sin  9^/Sin  02  =  n2/n^ 

Note  in  Figure  2  that  the  ray  is  bending  (retarded)  towards 
the  higher  n  value,  so  that  refractive  bending  is  toward 
media  of  higher  density. 

Since  the  density  of  air  generally  decreases  with 
altitude,  the  value  of  n  will  similarly  decrease.  Thus,  in 
a  ’’standard  atmosphere”  (a  term  used  to  describe  average 
meteorological  conditions),  EM  waves  do  not  travel  in  a 


straight  line  but  bend  slightly  downward  towards  the  earth's 
surface,  the  region  of  higher  n.  A  list  of  the  factors 
which,  according  to  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  [Ref.  1:  p.  20],  constitute  "standard" 
conditions  is  provided  in  Table  1,  while  Figure  3  is  a 
graphic  illustration  of  the  effect  of  index  of  refraction 
variation,  with  altitude,  on  EM  propagation.  When 
discussing  changes  of  the  index  of  refraction  in  the 
atmosphere,  it  is  important  to  realize  that  there  are  no 
abrupt  boundaries,  as  is  implied  in  Figure  3,  Changes  in  n 
generally  occur  over  relatively  large  distances,  when 
compared  to  the  wavelengths  affected.  It  is  these  gradual 
changes  that  create  the  smooth  bending  effect  shown  in 
Figure  4,  as  opposed  to  the  sharp  angles  shown  in  Figures  2 
and  3.  Figure  4  also  shows  another  of  the  effects  of 
refraction  in  the  troposphere,  i.e.,  extension  of  the 
horizon.  In  essence,  an  EM  ray  directed  above  the  geometric 
horizon  is  refracted  and  bent  slightly  around  the  surface  of 
the  earth,  striking  the  surface  at  a  distance  greater  than 
the  geometric  horizon.  For  example,  J.  Lake  [Ref,  2:  p.  22] 
states  that,  for  microwave  rays  propagating  in  a  standard 
atmosphere,  the  horizon  is  extended  beyond  the  geometric 
horizon  by  approximately  30%. 

Should  the  decrease  in  n,  with  altitude,  be 
sufficiently  great,  a  non-standard  atmosphere  will  exist, 
where  EM  rays  will  bend  downward  more  strongly  than  in  a 
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Figure  2.  Snell’s  Law. 


TABLE  1 

METEOROLOGICAL  CONDITIONS  FOR  THE  STANDARD  ATMOSPHERE 

At  altitude  z  =  0  (sea  level) 

Pressure(P)  =  1013.25  mb 

Temperature  (T)  =  15oc  (288. 15K) 

2 

Acceleration  (g)  =  9.8  m/s 
Density  (p)  =  1.225  Kg/m^ 


Direct  Path 


Refracted  Ray 
Path 


Surface  of  Earth 


Figure  3.  Path  of  microwave  through  the  troposphere. 

The  radar  transmitter  is  located  on  the 
earth’s  surface  and  transmitting  along  path 
A.  The  atmosphere  is  assumed  to  consist  of 
concentric  shells  each  with  a  constant 
refractive  index. 


Refracted  Ray  Path 
—  Direct  Path 


Figure  4.  Extension  of  the  radar  horizon  in  the 

presence  of  refraction.  A  is  the  radar  hori¬ 
zon  in  the  absence  of  refraction  (the 
geometric  horizon)  and  B  is  the  radar 
horizon  in  the  presence  of  refraction. 


standard  atmosphere.  This  can  lead  to  the  most  dramatic 
form  of  atmospheric  refraction,  called  "trapping",  where  a 
horizontally  launched  ray  is  bent  sufficiently  so  that  it 
remains  parallel  to  the  earth's  surface,  extending  the 
horizon  to  great  distances.  If  we  consider  the  atmosphere 
to  consist  of  layers  of  varying  refractive  index,  n,  there 
is  a  critical  angle  of  incidence  needed  to  refract  an  EM 
wave  completely  back  from  a  layer  of  lower  n  value  to  a 
layer  of  higher  n  value.  At  all  angles  of  incidence  smaller 
than  the  critical  angle ,  EM  rays  are  "trapped"  by  the  layer 
of  higher  n.  For  all  angles  larger  than  the  critical  angle, 
the  rays  will  be  bent,  but  not  so  sufficiently  as  to  be 
trapped.  Thus,  a  region  just  above  the  trapped  rays  is 
formed,  where  fewer  EM  rays  exist  than  if  standard  atmos¬ 
pheric  conditions  prevailed.  This  region  is  referred  to 
as  the  radio  or  radar  hole,  the  basic  concept  of  which  is 
shown  in  Figure  5.  According  to  John  Beach  [Ref.  3:  p.  82], 
the  use  of  the  word  "hole"  is  not  entirely  accurate,  since, 
under  normal  circumstances,  there  is  no  distinct  critical 
angle,  and  leakage  of  EM  energy  into  this  region  can  occur 
by  means  other  than  direct  propagation  of  the  signal.  There 
is,  however,  significantly  less  energy  present  in  these 
radar/radio  holes  than  would  be  present  if  no  variation  in  n 
existed. 

The  critical  angle  for  trapping,  9^,  as  derived  from 
Snell’s  law,  is  that  angle,  9,  which  causes  02=90.  Thus, 
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Sin  9  =  n-Zn^ 

c  i  1 


(2.3) 


Since  n  only  varies  from  about  1.00025  to  1.004  in  the 
troposphere;, n 2 Zn^  =  .9998  and  9^  *  89°  Hence,  as  shown  in 
Figure  6,  EM  rays  must  be  almost  parallel  to  the  layer  in 
order  to  be  trapped. 

The  index  of  refraction  of  the  atmosphere  does  not 
always  decrease  with  altitude.  Should  n  increase  with 
altitude,  EM  rays  will  be  bent  upward  instead  of  downward. 

More  will  be  said  about  this  latter  case,  in  discussions  on 
subrefraction. 

B.  REFRACTIVITY 

1.  Refractivity  (N) 

The  tropospheric  index  of  refraction,  n,  is  a 
characteristic  of  lower  atmospheric  properties.  Specifically, 
at  radio  and  microwave  frequencies ,  n  is  a  function  of 
temperature,  humidity,  and  pressure,  and  is  related  to  these 
parameters  by; 

(n-l)xio®  =  (77.6 PZT  -  5.6  eZT  +  3.75x10^  eZT^)  =  N, 

(2.4) 

where  P  is  the  barometric  pressure  in  millibars,  e  is  the 
partial  pressure  of  water  vapor  in  millibars,  T  is  the 
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absolute  temperature  in  Kelvin,  and  N  is  the  refractivity , 
a  paTcuneter  defined  for  the  sake  of  convenience. 

The  second  term  in  Equation  2.4  is  often  neglected, 
since  its  contribution  to  the  calculation  is  very  small 
compared  with  that  of  the  other  two  terms,  and  the 
coefficient  on  the  third  term  is  usually  modified  so  that 
the  actual  equation  is  generally: 

N  =  (77.6  P/T  +  3.73x10^  e/T^)  (2.4A) 


Beach  [Ref.  3:  p.  75]  states  that  since  the  value  of 
n  at  the  earth's  surface  is  usually  between  1.00025  and 
1.0004,  the  value  of  N  generally  lies  between  250  and  400. 

At  optical  frequencies  the  calculation  of  refrac¬ 
tivity  is  modified  to  reflect  the  fact  that  water  vapor  has 
a  negligible  effect.  Thus,  for  optical  frequencies: 


N  =  77.6P/T 


(2.5) 


The  significance  of  refractivity  is  not  in  the 
individual  value  of  M  for  a  given  altitude,  but,  rather,  in 
the  gradient  of  N,  (dN/dZ),  where  Z  is  height.  This  is  due 
to  the  fact  that  the  electromagnetic  wave  path  curvature 
depends,  inversely,  on  the  gradient  of  N.  For  microwave 
frequencies  and  a  standard  atmosphere,  M  decreases  at  a 
fairly  steady  rate  of  approximately  4x10  M  per  meter  of 
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altitude.  It  is  this  characteristic  that  makes  it  convenient, 
in  computations  involving  a  standard  atmosphere,  to  replace 
the  actual  earth's  radius  with  a  4/3  earth  radius.  Using 
this  geometry  EM  waves  can  be  shown  to  propagate  in  straight 
lines  rather  than  in  curved  or  refracted  lines,  as  shown  in 
Figure  7.  The  value  of  using  the  4/3  earth  radius  lies  in  its 
simplification  of  distance  computations  above  the  earth's 
surface,  and  in  ray  tracing.  This  geometry  has  been  widely 
used  in  the  modeling  of  radio  communication  and  radar  propa¬ 
gation.  However,  it  must  be  remembered  that  the  4/3  earth 
representation  is  only  an  approximation  of  conditions  that 
vary  with  meteorological  conditions.  The  actual  value  may 
vary  from  6/5  to  4/3  earth  radius. 

A  vertical  refractivity  gradient  causes  EM  rays  to 
be  bent  with  a  radius  of  curvature  given  by: 

r  =-10®/(dN/dZ) ,  (2.6) 

where  the  curvature  is  in  a  vertical  plane. 

Horizontal  N  gradients  are  normally  not  significant 
in  EM  propagation  and  are,  therefore,  usually  assumed  to  be 
zero.  In.  some  circumstances,  however,  bodies  of  water  either 
partially  or  completely  surrounded  by  land  may  create 
horizontal  N  gradients  which  will  have  some  effect  on  EM 
transmissions  through  the  atmosphere  above  the 
land-water  boundary  region. 
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Earth  radius 


Figure  7a.  Bending  of  antenna  beam  due  to  refraction 
by  the  troposphere. 


4/3  Earth  radius 


Figure  7b.  Shape  of  antenna  beam  in  4/3  earth 
representation. 


2 .  Anomalous  Propagation,  Mirage,  S  Optical  Turbulence 

Refractivity  effects  are  not  the  same  across  the  EM 
spectrum.  As  previously  shown,  the  calculation  of  M  for  the 
radio  and  radar  portion  of  the  spectrum  is  different  than 
that  for  the  optical  portion  (where  the  effect  of  humidity 
is  ignored).  Since  the  amount  of  bending  is  a  function  of 
the  gradient  of  N,  it  is  to  be  expected  that  the  amount  of 
bending  is  different  for  optical  and  microwave  rays.  This 
can  be  seen  in  Figure  8.  Note  the  difference  in  the  bending 
of  the  optical  grazing  ray  and  the  radio  grazing  ray 
vis-a-vis  each  other  and  the  geometric  horizon. 

Regions  of  strong  refractivity  gradients  that  are 
much  larger  in  size  than  EM  wavelengths  cause  anomalous 
propagation.  Smaller  scale  regions  with  large  N  gradients 
also  exist,  in  the  form  of  mirages  and  atmospheric  turbu¬ 
lence.  Certain  forms  of  atmospheric  refraction,  when 
located  within  the  first  several  meters  of  the  surface  of 
the  earth,  and  resulting  from  rapid  changes  in  the 
temperature,  can  lead  to  mirage  conditions  which  have  a 
significant  effect  on  optical  transmission.  Mirage 
formation  is  a  form  of  anomalous  propagation,  but,  for  the 
sake  of  simplicity,  the  two  will  be  discussed  separately. 
They  are  best  differentiated  by  relative  size.  Turbulence 
in  the  atmosphere  creates  fluctuations  in  the  air  index  of 
refraction,  resulting  from  inhomogeneities  that  vary  in  size 
from  those  of  almost  microscopic  dimensions  to  those 
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Geometric 
(straight)  ray 


Grazing  ray 
(optical) 

Grazing  ray 
(radio) 


A  Geometric  horizon 
B  Optical  horizon 
C  Radio  horizon 


Figure  8.  Comparison  of  effect  of  refraction  on 
optical  ray  and  radio  ray. 


hundreds  of  meters  in  diameter.  The  wavelengths  in  the  IR 
and  visible  range  are  most  affected  by  turbulence,  due  to 
the  small  size  of  the  refractive  inhomogeneities  relative  to 
the  wavelengths  of  radio  and  radar  signals.  Each  of  these 
three  conditions  will  be  individually  addressed  in  the 
following  sections,  beginning  with  a  discussion  on  anomalous 
propagation. 

C.  ANOMALOUS  PROPAGATION 

Anomalous  propagation  is  a  misnomer,  referring  to  any 
propagation  other  than  that  in  a  standard  atmosphere. 
Anomalous  propagation  is  generally  considered  to  be  divided 
into  three  forms;  superrefractive ,  subrefractive ,  and 
ducting.  Although  the  name  implies  rare  or  infrequent,  the 
refractive  conditions  listed  above  are  far  from  infrequent. 

C.  G.  Purves  [Ref.  4;  pp,  13-14]  states  that  in  parts  of  the 
world,  and  during  certain  seasons  ducting  or  trapping  occurs 
more  than  75%  of  the  time.  Although  the  frequency  of 
occurrence  of  various  propagating  conditions  will  be 
addressed  in  more  detail,  let  it  suffice  to  say  that 
anomalous  propagation  is  far  from  uncommon. 

Another  point  to  address  is  the  fact  that  superrefrac¬ 
tion,  subrefraction,  and  ducting  occur  on  a  macroscopic 
scale.  They  are  characteristic  of  relatively  large 
seasonal  and/or  atmospheric  conditions,  some  of  which  will 
be  discussed  at  the  end  of  this  section.  This  characteristic 
roughly  separates  anomalous  propagation  from  the  conditions 
leading  to  mirages. 
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1 .  Superrefraction,  Subrefraction,  and  Ducting  Defined 
Equation  2.4  shows  that  if  water  vapor  pressure 
decreases  rapidly  with  height  or  if  the  temperature  increases 
with  height,  N  will  decrease  with  height  at  a  more  rapid  rate 
than  under  standard  atmospheric  conditions.  Thus,  the  value 
of  dN/dZ  becomes  more  negative.  For  small  gradients  a 
superrefractive  condition  exists,  while  when  EM  rays  are 
curved  parallel  to  the  earth's  surface  trapping  occurs  and 
the  phenomenon  of  ducting  exists.  Should  the  index  of 
refraction  increase,  instead  of  decrease  with  altitude,  a 
subrefractive  condition  exists.  The  meteorological  conditions 
leading  to  these  forms  of  anomalous  propagation  will  be 
discussed  in  detail  in  Section  II. C. 4.  of  the  thesis. 

Superrefraction  implies  a  gradient  greater  than  that 
for  a  standard  atmosphere,  but  not  strong  enough  to  form  a 
duct.  EM  waves  propagating  in  this  region  will  be  bent 
downward  more  than  in  normal  (standard  atmosphere)  conditions 
but  not  sufficiently  for  trapping  to  occur.  Superrefractive 
conditions  will  result  in  propagation  distances  beyond 
those  expected  in  a  standard  atmosphere. 

Trapping,  or  ducting,  occurs  when  dn/dZ  is  less  than 
or  equal  to  -157.  Signals  traveling  in  a  duct  propagate  in 
a  manner  similar  to  propagation  in  a  "leaky"  waveguide. 

Ducting  is  categorized  into  three  subgroups:  evaporative 
ducts,  surface  ducts,  and  elevated  ducts.  Evaporative  ducts, 
most  commonly  created  by  evaporation  off  of  the  surface  of 
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water,  are  present  only  to  a  maximum  height  of  about  thirty 
meters  and  can  be  much  thinner.  These  ducts  are  present  over 
the  surface  of  the  oceans  of  the  world  to  some  degree  most  of 
the  time. 

J.  Beach  [Ref.  3:  pp.  78-80]  states  that  surface 
ducts,  which  are  usually  stronger  and  generally  extend  higher 
in  altitude  than  evaporative  ducts,  vary  in  thickness  but  are 
almost  always  less  than  one  kilometer  and  usually  are  300 
meters  or  less  in  thickness.  Surface  ducts  are  characterized 
by  the  lower  boundary  being  the  earth's  surface.  Surface 
ducts,  along  with  elevated  ducts,  occur  much  more  frequently 
over  land  than  do  evaporation  ducts. 

Beach  [Ref.  3;  p.  80]  in  his  discussion  of  elevated 
ducts  states  that  elevated  ducts  begin  at  altitudes  of  from 
close  to  zero  meters,  absolute  altitude,  to  6  kilometers  but 
generally  no  higher  than  3  kilometers.  The  thickness  of 
these  ducts  is  similar  to  the  surface  duct. 

When  the  index  of  refraction  increases  with  altitude , 
a  subrefr active  layer  is  created.  This  type  of  layer  is 
normally  created  by  heating  of  the  surface  and  thus  exists 
in  only  a  thin  layer  close  to  the  earth's  surface  as  shown 
in  Figure  9.  EM  waves  entering  such  a  layer  from  above  will 
be  bent  upward,  as  will  rays  traveling  within  the 
subrefractive  region.  Subrefraction  is  responsible  for  some 
of  the  more  common  forms  of  optical  mirage  (see  Section  D). 


Surveillance 
aircraft  " 


Figure  9.  Description  of  Subrefraction.  A 

subrefractive  layer  bends  EM  energy 
upward,  away  from  the  surface  of  the 
earth. 


2.  Means  of  Describing  Refraction 


As  shown  previously,  the  refractivity ,  N,  is  a 
parameter  that  normally  decreases  with  height  and  is 
therefore  not  very  useful  in  depicting  refraction 
graphically.  Two  indices  that  are  more  commonly  used  for 
this  purpose  are  B  and  M  units. 

B  is  defined  by  K.  Davidson  [Ref.  5:  pp.  4-2  to  4-4] 
by  adding  40  N  units  per  kilometer  to  all  N  values 

B  =  N  +  40Z,  (2.7) 

z 

where  is  the  value  of  refractivity  at  height  Z  in 
kilometers.  Thus,  the  B  gradient  is  zero  for  standard 
atmospheric  conditions. 

Davidson  goes  on  to  describe  the  "modified 
refractivity"  M,  as  being  defined  by  adding  157  N  units 
per  kilometer  to  M,  thus 

M  =  N  +  157Z  (2.8) 

z 

The  M  gradient  is  zero  when  EM  ray  curvature  equals  the 
earth’s  curvature.  Since,  in  a  standard  atmosphere,  M 
increases  with  height,  layers  In  which  trapping  can  occur 
will  be  readily  apparent  by  the  negative  M  gradient.  Figure 
10  shows  the  relationship  between  N,  B,  and  M  units. 
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Figure  11  shows  the  Navy  Integrated  Refractive  Effects 
Prediction  System  (IREPS)  classification  of  refractive 
conditions.  See  also  Appendix  A  for  more  information  on 
IREPS. 

3 .  Duct  Characteristics 

Before  describing  the  causes  of  ducting  it  is 
necessary  to  explain  some  of  the  characteristics  of  ducts. 

To  begin,  an  important  distinction  must  be  made  between  the 
trapping  layer  and  the  duct.  As  shown  in  Figure  12,  the 
trapping  layer  is  the  region  where  dM/dZ  is  less  than  zero. 
Consequently  the  EM  rays  are  bent  downward  with  a  curvature 
greater  than  that  of  the  surface  of  the  earth,  causing  them 
to  be  trapped  between  the  layer  and  the  bottom  of  the  duct . 

A  duct  is  the  region  below  the  top  of  the  trapping  layer 
where  the  EM  rays  are  confined  or  trapped  and  includes  the 
trapping  layer  itself.  The  EM  rays  are  trapped  because  they 
are  partially  confined  between  the  top  of  the  trapping  layer 
(duct)  and  the  bottom  of  the  duct.  The  trapping  layer  is 
always  at  the  top  of  the  duct.  The  duct  bottom  can  be 
located  on  the  earth’s  surface  (surface  duct),  or  can  be 
elevated  above  the  ground  (elevated  duct). 

The  top  of  the  duct  is  where  the  slope  of  dM/dZ 
changes  from  negative  to  positive.  At  this  height  M  is  no 
longer  decreasing  with  height  but  begins  to  increase  in 
value,  as  shown  in  Figure  13,  for  surface  and  elevated  ducts 
The  bottom  of  the  duct  is  located  by  dropping  a  line 


t! 


32 


Elevated  Duct 


vertically  from  the  location  of  the  top  of  the  duct,  in 
Figure  13,  to  the  earth's  surface  (for  a  surface  duct)  or 
to  the  intersection  with  the  M  curve  (for  an  elevated  duct). 
The  optimum  coupling  height  into  the  duct  is  that  point  where 
dM/dZ  goes  negative.  This  is  the  point  where  a  radar/radio 
treinsmitter  will  achieve  the  greatest  extension  of  range  due 
to  the  trapping  effect  of  the  duct.  The  strength  of  the 
duct  is  determined  by  the  most  negative  value  that  dM/dZ 
attains  and  by  the  thickness  of  the  negative  layer.  The 
more  negative  the  value  the  stronger  the  duct.  Note  in 
Figure  12  that  a  surface  duct  can  be  caused  by  a  surface 
trapping  layer  or  an  elevated  trapping  layer.  The  bottom  of 
the  duct  is  the  location  where  the  downward  bent  EM  rays  are 
turned  back  upward  again.  It  is  apparent  from  this 
discussion  that  ducts  can  vary  considerably  in  thickness 
and,  depending  on  the  slope  of  the  trapping  layer,  they  can 
also  vary  considerably  in  intensity  (strength  of  the  duct). 

As  will  be  described  later,  not  all  EM  frequencies 
are  trapped  within  the  duct  and  only  those  EM  rays  entering 
the  duct  or  traveling  through  the  duct  at  a  very  shallow 
angle  will  be  trapped.  Note  that  for  the  EM  ray  to  be  most 
effectively  trapped  the  transmitter  must  be  within  the  duct 
although  EM  signals  can  be  coupled  into  the  duct  and  trapped 
within  it  from  an  antenna  located  near  but  outside  the  duct. 
In  addition,  note  that  the  trapped  EM  ray  is  either 
refracted  back  upward  after  reaching  the  bottom  of  an 


elevated  duct  or  reflected  upward  off  the  surface,  as  in  the 
case  of  a  surface  or  evaporation  duct.  In  both  of  these 
cases  at  least  part  of  the  energy  of  the  EM  ray  is  lost  due 
to  surface  absorption  and  diffuse  scattering  (surface  and 
evaporation  ducts)  or  lost  through  leakage  through  the 
bottom  of  the  duct  (elevated  duct).  Further,  part  of  the  EM 
ray  is  lost  by  leakage  through  the  top  of  all  ducts. 

Diffuse  scattering  and  absorption  are  greater  over 
land  than  over  water.  M.  L.  Meeks  [Ref.  6:  pp.  13-26]  states 
that  these  two  surface  characteristics  result  in  the  surface 
reflection  coefficient,  the  value  of  which  varies  as  a 
function  of  frequency,  angle  of  incidence  of  the  reflected 
ray  to  the  surface,  and  the  terrain  surface  itself.  The 
lower  the  angle  of  incidence,  the  smoother  the  terrain,  the 
more  moist  the  surface,  and  the  lower  the  frequency,  the 
greater  the  amount  of  reflection  of  the  EM  ray.  The 
reflection  coefficient  for  water  varies  from  approximately 
.85  to  .99  and  is  largely  a  function  of  the  sea  state.  The 
reflection  coefficient  over  land  is  much  more  variable, 
ranging  from  as  low  as  .15  to  higher  than  .9.  The 
reflective  properties  of  various  land  surfaces  and  types  of 
vegetation  are  still  not  completely  understood. 

The  amount  of  EM  energy  lost  as  a  signal  travels 
through  a  duct  is  a  complicated  function  of  frequency,,  duct 
characteristics,  and  in  the  case  of  surface  and  evaporation 
ducts,  the  terrain’s  reflective  characteristics.  Due  to 
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these  losses  and  those  resulting  from  atmospheric  scattering 
and  absorption  it  can  be  seen  that  there  is  a  limit  to  just 
how  far  an  EM  signal  can  be  propagated  in  the  atmosphere. 

It  is  apparent  from  this  discussion  that  ducting  of  signals 
is  less  effective  in  increasing  EM  signal  range  over  land 
than  over  water. 

As  mentioned  earlier,  not  all  frequencies  are  trapped. 
It  is  the  thickness  of  the  duct  itself  that  is  the  determining 
characteristic  as  to  which  frequencies  are  trapped.  Duct 
thickness  is  therefore  used  to  calculate  the  lower  limit  of 
those  frequencies  that  are  strongly  trapped,  as  shown  in  the 
following  formula  from  Davidson  [Ref.  5:  pp.  4-15]: 


(3.6033  X 


(2.9) 


Frequencies  below  f^  are  still  affected  by  the  duct  to  a 

somewhat  lesser  extent.  Table  2  shows  values  of  f  for 

m 

various  duct  thicknesses. 

It  is  important  at  this  point  to  reiterate  that  not 


only  does  the  thickness  of  the  duct  limit  the  frequencies 


trapped  within  a  duct,  but  so  too  doe^  xhe  angle  of 
incidence  of  the  EM  rays  entering  the  trapping  layer.  As 
mentioned  earlier,  only  EM  rays  incident  at  very  shallow 


angles  relative  to  the  horizon  will  be  refracted  enough  to 
be  trapped.  In  actuality,  angles  above  1  degree  -from  the 
earth’s  horizon  are  not  trapped  and,  according  to  Skolnik 
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[Ref.  7:  p.  450],  refraction  is  troublesome  only  at  angles 
of  elevation  of  5  degrees  or  less  from  the  horizon  in  most 
microwave  applications.  Figure  5  provides  a  simplified 
example  of  the  effect  of  angle  of  incidence  on  refraction, 
trapping,  and  creation  of  a  radar/radio  hole. 

^ .  Causes  of  Superrefraction,  Ducting,  S  Subrefraction 

Superrefractive ,  subrefractive  and  ducting  conditions 
are  created  by  variations  in  water  vapor  and/or  temperature 
in  the  atmosphere,  brought  on  by  a  multitude  of  climatic  and 
topographic  conditions. 

Evaporation  ducts  are  created  above  the  surface  of 
large  bodies  of  water  by  the  rapid  decrease  in  humidity  above 
the  water's  surface.  The  humidity  of  the  air  immediately 
above  the  water  is  at  the  saturation  point,  but  decreases 
quickly  with  height.  Within  the  first  thirty  meters,  an 
ambient  level  of  humidity  is  reached,  dependent  on  the 
climatic  conditions  prevailing.  This  rapid  decrease  in 
humidity  causes  M  to  decrease  with  altitude  to  a  point  where 
M  will  reach  a  minimum  (the  top  of  the  duct)  and  then 
increase  with  altitude.  Meeks  [Ref.  6:  p.  9]  states  that 
over  land  a  similar  situation  is  created  after  a  rainfall, 
when  the  water  evaporates  from  the  ground.  Snow,  and 
particularly  melting  snow,  can  also  lead  to  the  creation  of 
an  evaporation  duct  over  land.  Overland  evaporation  ducts, 
however,  are  generally  of  short  duration.  Evaporation  ducts 
will  also  occur  over  the  continental  land  masses  above 


larger  bodies  of  water.  For  the  Navy,  the  occurrence  and 
effects  of  evaporative  ducts  are  very  important  and  largely 
responsible  for  the  development  of  the  IREPS  program  which 
models  the  effects  of  superrefraction  and  ducting  on  EM 
communications  and  on  surveillance  systems . 

Superrefractive  conditions  and  surface  and  elevated 
ducts  are  created  by  the  same  climatic  conditions .  Which 
type  of  propagation  condition  occurs  depends  primarily  on 
where  the  transition  between  two  differing  air  masses  forms 
and  on  the  strength  of  the  resulting  refractivity  gradient. 
Temperature  inversions  are  the  primary  cause  of  superrefraction 
and  ducting.  As  would  be  expected,  superrefractive  conditions 
are  far  more  common  than  ducting. 

A  temperature  inversion  is  a  region  of  the  atmosphere 
where  there  is  an  increase  in  temperature  with  height.  As 
shown  in  Equation  2.4,  a  temperature  inversion  must  be  very 
strong  to  cause  superrefraction  by  itself.  If,  however,  a 
water  vapor  gradient  is  combined  with  the  temperature 
inversion,  superrefraction  and  ducting  will  more  easily 
result.  Thus,  it  is  no  surprise  to  find  that  superrefraction 
and  ducting  are  more  common  over  v/ater  (particularly  oceans) 
than  over  land. 

Ducting  is  characterized  by  the  upper  air  being  much 
warmer  and  drier  than  the  underlying  air.  There  are  several 
meteorological  conditions  that  can  lead  to  this  situation. 

Over  land,  the  most  common  ducting  is  caused  by  radiation  of 
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heat  from  the  earth’s  surface  on  clear  calm  nights.  This 
surface  cooling  creates  a  temperature  inversion  and  a  sharp 
decrease  in  moisture  with  height  which  will  often  lead  to  a 
surface  fog  in  summer.  Refractivity  gradients  are  made  even 
stronger  if  the  ground  is  moist,  a  situation  most  common  in 
summer  and  in  winter  over  snow.  Meeks  [Ref.  6:  p.  9] 
elaborates  on  the  fact  that  in  desert  regions  the  temperature 
inversion  created  at  night  can  be  significant  enough,  in  and 
of  itself,  to  create  superrefractive  conditions  and  even 
ducting.  Ducting  is  therefore  most  common  over  land  at  night, 
and  will  usually  weaken  or  disappear,  along  with  the  weakening 
and  disappearance  of  the  temperature  inversion,  during  the  day. 

Subsidence,  the  sinking  of  air  from  above,  also  is  a 
common  cause  of  ducting.  As  the  air  sinks  it  is  warmed  by 
compression  capping  the  cooler  air  below  and  a  strong 
inversion  is  formed.  If  the  lower  level  air  has  a  high 
moisture  content,  a  strong  humidity  gradient  will  also  be 
formed.  Since  subsidence  is  most  often  associated  with  high 
pressure  areas,  this  type  of  ducting  is  generally  a  fair 
weather  phenomena,  characterized  by  light  winds,  and  clear 
or  low  stratus  skies.  The  presence  of  fair  weather  low 
stratus  clouds  that  show  a  marked  discontinuity  at  the  top 
of  the  cloud  mass  will  often  indicate  the  presence  of  an 
inversion  and  a  good  possibility  that  superrefraction  or 
ducting  is  present.  Figure  14  shows  temperature  inversions 
created  by  radiation  of  heat  and  by  subsidence. 


Initial 


Initial 


Surface  Fog 


Stratus 

naensation 


Level 


Figure  14. 


Types  of  Inversion:  (A)  Subsidence  reaching 
the  surface;  (B)  Subsidence  above  a  mixed 
surface  layer;  (C)  Ground  Inversion  caused 
by  cooling  of  surface  air;  (D)  Inversion 
caused  by  cooling  and  mixing  of  surface 
lave 


In  temperate  regions  in  the  low  and  middle  latitudes, 
ducting  occurs  most  often  in  summer  and  does  not  usually 
occur  during  bad  weather.  Storms  are  characterized  by  the 
lower  atmosphere  being  mixed  by  precipitation  and  winds, 
resulting  in  a  more  homogenous  atmosphere.  In  arid  regions 
humidity  will  have  a  lesser  impact  on  the  formation  of 
anomalous  propagation  conditions  but  in  temperate  regions 
the  air  temperature  may  be  high  enough,  particularly  in 
summer,  for  the  water  vapor  component  of  the  atmosphere  to 
be  important.  Seasonal  variations  in  refractive  conditions 
are  strongest  over  large  land  masses  and  generally  weakest 
over  large  bodies  of  water  (oceans). 

As  described  earlier,  ducts  over  land  generally  are 

weakest  or  break  up  entirely  during  the  warmest  part  of  the 

2 

day.  This  diurnal  characteristic  of  the  land  does  not, 
however,  occur  significantly  over  large  bodies  of  water. 

W.  Thompson  [Ref.  8:  p.  3]  additionally  notes  that  in 
temperate  regions ,  diurnal  variations  are  more  pronounced 
during  summer  months  and  less  so  in  the  winter.  C.  Samson 
[Ref.  9:  p.  6]  states  that  due  to  the  rapidly  changing 
thermal  structure  of  the  air  near  the  earth's  surface  at 


The  changes  in  the  thermal  structure  are  due  to  the 
effects  on  the  atmosphere  of  the  heating  and  cooling  of  the 
earth's  surface. 
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sunset  and  sunrise,  the  most  extreme  changes  in  refractivity 
will  occur  at  these  times. 

Skolnik  [Ref.  7;  p.  452]  states  that  the  main 
exception  to  the  fair  weather  rule  for  ducting  is  the 
creation  of  ground  ducts  by  thunderstorm  downdrafts.  These 
will  usually  last  on  the  order  of  1  hour  or  less  and  are 
small  in  spatial  extent. 

Coastal  regions  develop  their  own  unique  ducting 
conditions.  As  described  in  Davidson's  notes  [Ref.  5: 
pp.  7-30  to  7-31]  the  presence  of  ducting  in  coastal  regions 
and  straits  is  quite  common  and  is  due  to  the  jutaposition  of 
marine  and  continental  air  masses  at  the  surface.  Offshore 
flow  of  warm  dry  air  would  lead  to  warm,  dry,  continental 
air  overlying  the  cooler  and  more  moist  marine  air.  Offshore 
flow  can  be  caused-by:  1)  local  land  to  sea  breezes  induced 
by  the  differential  heating  of  land  and  water  areas  and 
sea-breeze  recirculation,  and  by  2)  larger  scale  circulations 
in  which  an  offshore  flow  occurs  with  a  semi-stationary  high 
pressure  area  located  over  land.  Figure  15  shows  both  of 
these  conditions. 

In  the  latter  situation,  extremely  strong  ducts  can 
be  formed,  some  of  which  will  extend  well  inland  from  the 
coast.  The  Santa  Ana  winds  and  monsoon  circulations  are 
characteristic  of  this  type  of  offshore  flow. 

Subrefraction  is  the  rarest  of  the  anomalous 
propagation  categories.  The  only  regions  where  n  increases 
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Figure  15a.  Illustration  of  elevated  layer  occurrence 
due  to  local  differential  land-sea  heating 


Figure  15b.  Illustration  of  elevated  layer  due  to 
synoptic  scale  driven  offshore  flow. 
Note  adiabatic  heating  with  downs lope 
trajectorv. 


with  height  with  any  regularity  are  the  deserts  of  the  world 
and  particularly  those  that  are  near  large  bodies  of  water. 

The  sun  heats  the  desert  surface  causing  it  to  be  very  hot 
and  dry,  the  temperature  to  decrease,  and  the  water  vapor 
to  increase,  with  height.  This  subrefractive  condition  can 
be  made  even  stronger  by  an  onshore  breeze  which  is 
characteristic  of  a  high  pressure  area  located  at  sea.  The 
onshore  breeze  carries  cool  moist  air  that  rises  over  the 
land . 

Subrefractive  conditions  tend  to  follow  a  diurnal 
cycle  and  are  strongest  in  the  afternoon  and  usually  disappear 
at  night.  In  semi-arid  and  desert  regions  it  is  not  uncom¬ 
mon  to  see  strong  subrefractive  conditions  forming  in  the 
afternoon,  due  to  solar  heating  of  the  surface,  and  super- 
refractive  conditions  develop  at  night,  due  to  radiation 
of  heat  from  the  surface. 

5 .  Consequences  of  Anomalous  Propagation 

As  described  earlier,  the  effects  of  anomalous 
propagation  vary  with  the  wavelength  and  the  gradient  of  the 
refractive  index.  Increases  in  the  refractive  gradient 
(due  to  superrefraction  or  ducting)  result  in  an  extension 
of  the  radar  horizon,  causing  an  increase  in  radar  coverage 
beyond  that  expected  for  a  standard  atmosphere,  as  shown 
in  Figure  4.  Refraction  can  also  lead  to  errors  in  the 
measurement  of  range  and  elevation  angle  by  microwave 
systems.  An  elevation  angle  error  is  created  by  the  bending 


47 


of  the  EM  waves  away  from  a  straight  path  as  they  propagate 
to  the  target  and  back  as  shown  in  Figure  16.  Range  error  is 
created  by  the  increase  in  the  effective  path  length  due  to 
refraction  of  the  EM  rays  from  the  radar  transmitter  to  the 
target  and  back  to  the  radar  receiver.  Figure  17  shows 
range  error  as  a  function  of  the  apparent  elevation  angle. 

An  example  of  the  type  of  range  error  that  can  be  created  by 
atmospheric  refraction  is  that  experienced  by  the  FPS-16  air 
search  and  tracking  radar  at  Mount  Lemmon,  Arizona.  For 
this  radar,  range  errors  are  generally  between  30  and  40 
yards  over  a  150,000  to  200,000  yard  distance.  In  many 
instances  this  type  of  error  is  of  little  consequence,  but, 
for  accurate  tracking  of  a  target's  location  some 
compensation  for  the  range  error  must  be  made. 

Communications  signals  are  affected  by  refraction  in 
a  manner  similar  to  radar  signals,  leading  to  extended  ranges 
in  "line  of  sight"  communications.  However,  the  terrain 
over  which  the  signals  propagate  will  often  have  a  greater 
effect  on  communications  range  than  will  refraction  for 
ground  based  systems. 

The  effect  of  refraction  on  Direction  Finding  (DF) 
systems  can,  in  certain  circumstances,  be  significant. 
Harrington  and  Callaghan  [Ref.  10;  p.  3]  state  that  changes 
in  the  refractive  index,  particularly  if  occurring 
horizontally,  will  lead  to  the  same  type  of  error  in 
direction  finding  as  the  elevation  angle  error  in  a  radar, 
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Figure  17. 


Range  error  as  a  function  of  apparent 
elevation  angle,  where  apparent  elevation 
angle  is  defined  as  the  angle  from  the 
horizontal  that  the  radar  is  aimed  in  order 
to  locate  the  target  (as  opposed  to  the 
angle  from  the  horizontal  that  the  target 
is  actually  located).  CRPL  reference 
refractivity  atmosphere  (1958)  is 
assumed.  Note  that  the  closer  the  apparent 
elevation  angle  is  to  the  horizontal,  the 
greater  the  range  error. 


n 


V? 


as  shown  in  Figure  18.  An  example  of  a  situation  that  might 
result  in  a  horizontal  gradient  is  a  signal  transmitted 
across  a  body  of  water  where  the  air  is  more  humid  than  over 
land.  Similarly,  a  signal  traversing  a  wide  valley,  where 
the  air  is  much  drier  than  the  surroundings,  can  result  in 
an  error  in  location  of  a  transmitter. 

Multiple  modes  on  an  EM  signal  propagating  in  a  duct 
can  result  in  variations  in  the  strength  of  the  EM  signal  at 
various  points  along  the  path  due  to  interference.  This 
spatial  fading  can  be  caused  by  multipath,  or,  as  in  the 
above  example,  by  multimode  interference,  both  of  which  can 
result  in  a  condition  where  the  vector  sum  of  several 
components  of  a  signal  are  not  in  phase.  The  variation  in 
phase  is  a  result  of  the  difference  in  path  length  induced 
by  refraction  and/or  path  route.  According  to  Skolnik 
[Ref.  7;  p.  454]  fading  on  the  order  of  20  dB  can  occur  for 
radar  systems  experiencing  multimode  propagation.  Fading 
can  be  greatly  enhanced  by  the  presence  of  ground  and/or  sky 
waves  combining  with  the  refracted  signals  (multipath). 

Beach  [Ref.  3:  p.  75]  states  that  temporal  fading  will 
frequently  occur  near  the  radio  horizon  due  to  changes  in 
the  index  of  refraction  leading  to  large  changes  in  signal 
strength  in  the  horizontal  EM  signal  path. 

In  the  presence  of  ducting,  the  effects  previously 
described  become  more  acute.  The  angle  error  in  a  ground 
based  radar  will  become  greater  with  a  larger  refractive 
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Denser  Air  Pocket 


Figure  18.  Signal  Refraction  because  of  change  of 
humidity  over  a  lake  and  its  effect  on 
D.F.  accuracy. 


gradient  and  can  often  be  more  serious  than  the  range  errors 
induced.  For  those  VHF,  UHF,  and  microwave  signals  propa¬ 
gating  nearly  horizontal  to  the  earth's  surface,  trapping 
can  also  occur,  particularly  if  the  transmitting  antenna  is 
located  within  the  duct.  When  trapping  occurs,  the  range  of 
the  EM  signal  within  the  duct  is  greatly  extended  beyond 
that  expected  in  a  standard  atmosphere  as  shown  in  Figure 
19.  Thus,  for  surface  based  radar  signals  the  coverage  can 
be  extended  well  beyond  the  expected  radar  horizon  and,  for 
communications  signals,  extraordinary  ranges  have  been 
caused  by  the  presence  of  surface  based  ducts.  Examples  of 
VHF  signals  propagating  2,000  miles  without  the  use  of 
satellites  are  well  known.  Likewise,  electronic  support 
measures  (ESM)  detection  of  these  signals  is  greatly  enhanced 
by  the  extension  of  the  range  of  the  propagating  signal.  For 
ground  based  DF  systems,  trapped  signals  can  be  a  mixed 
blessing.  Although  the  range  at  which  a  signal  can  be 
detected  is  extended,  the  accuracy  of  the  direction  finding 
can  be  degraded  by  the  very  nature  of  the  ducted  transmission. 

The  existence  of  radar/radio  holes  is  indicated  by 
regions  in  the  atmosphere  where  EM  signal  strength  is  weaker 
in  the  presence  of  a  duct  than  it  would  be  in  free  space. 

The  variables  that  dictate  the  existence  and  characteristics 
of  a  radar  or  radio  hole  are:  the  radar  or  radio  transmitter 
height  relative  to  the  duct,  the  radar  target  or  ESM 
receiver  height  relative  to  the  duct,  the  operating 
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Figure  20.  Coverage  Diagram  for  a  9000  MHz  airborne 
early  warning  radar  transmitting  within 
an  elevated  duct  located  bewteen  2,590 
feet  and  3,220  feet  above  the  surface. 

Note  the  extended  coverage  in  the  duct  and 
the  partial  loss  of  coverage  just  above 
the  duct.  This  coverage  diagram  like  those 
in  Figure  20  was  produced  by  the  Navy  IREPS. 
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.  Coverage  Diagram  for  a  9000  MHz  airborne 
early  warning  radar  transmitting  within 
an  elevated  duct  located  bewteen  2,590 
feet  and  3,220  feet  above  the  surface. 

Note  the  extended  coverage  in  the  duct  and 
the  partial  loss  of  coverage  just  above 
the  duct.  This  coverage  diagram  like  those 
in  Figure  20  was  oroduced  by  the  Navy  IREPS 


farther  away  from  the  aircraft  until  a  point  is  reached 
where  the  hole  is  no  longer  apparent. 

Electronic  countermeasures  (ECM) ,  ESM,  and  DF  systems 
performance,  as  well  as  air  to  air  communications  and  air  to 
ground  communications  are  also  affected  by  elevated  ducts, 
with  the  effect  being  most  likely  at  higher  frequencies.  It 
has  also  been  noted  that  the  propagation  of  EM  signals  from 
transmitters  located  well  below  the  elevated  duct  is  greater 
in  the  duct  than  expected.  One  explanation  mentioned  by 
Skolnik  [Ref.  7:  p.  453)  postulates  that  the  EM  signals  are 
being  scattered  into  and  out  of  the  elevated  duct  by 
variations  in  the  index  of  refraction  profile  due  to  the 
"wavy”  top  of  the  duct,  causing  it  to  be  similar  to  a 
nonuniforro  waveguide.  Figure  21  shows  examples  of  the 
effect  that  the  position  of  the  aircraft  (whether  radar  or 
jamming)  relative  to  the  elevated  duct  has  on  the  trapping 
of  EM  waves  and  on  the  effect  of  the  radar/radio  hole  on  EM 
waves . 

For  radar  systems,  the  presence  of  surface  or 
evaporation  ducts  can  lead  to  an  increased  surface  clutter 
return,  both  in  magnitude  and  in  coverage  extent,  masking 
targets  over  a  much  greater  range.  An  example  of  this 
effect  is  shown  in  Figure  22.  As  described  by  Beach 
[Ref.  3:  p.  85],  if  the  clutter  return  is  stronger  than  the 
target  return  at  the  same  range,  detection  of  the  target  may 
be  seriously  degraded  or  impossible.  The  strength  of  the 
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Figure  72 


.  Air  search  radar  geometry  showing  possible 
clutter  return  from  rough  surface  at  the 
same  range  as  the  air  target  for  a  surface 
duct. 


clutter  return  is  partially  dependent  on  the  strength  of  the 
duct.  MTI  radars  can  be  seriously  affected  by  the  enhanced 
strength  of  the  clutter  return  which  may  overcome  the  radar's 
ability  to  separate  out  a  moving  target  from  stationary 
clutter . 

The  fading  phenomena  described  earlier  can  also 
become  more  severe  in  the  presence  of  ducting  and  can  lead 
to  the  "skip"  effect  in  the  presence  of  a  surface  duct.  This 
phenomenon  is  characterized  by  periodic  spatial  gaps  in  radar 
coverage  or  communications  signals.  These  spatial  intervals 
give  the  appearance  of  rings  on  the  radar  display.  Figure  23 
shows  the  ray  paths  that  characterize  skip  rings.  Note  that 
the  appearance  of  skip  rings  is  dependent  on  the  surface 
reflectivity  as  well  as  the  strength  of  the  surface  duct 
and  the  angle  at  which  the  EM  ray  is  launched.  Since  water 
generally  has  a  higher  surface  reflectivity  than  land  this 
skip  effect  is  more  significant  over  the  oceans  of  the  world. 

Subrefractive  conditions  can  cause  the  creation  of 
blind  spots  within  the  subrefractive  zone  for  EM  emitting 
aircraft  (airborne  early-warning  and  EW  aircraft)  located 
above  the  zone,  as  shown  in  Figure  9.  Aircraft  traveling 
close  to  the  earth's  surface,  below  the  top  of  the  sub¬ 
refractive  zone,  can  take  advantage  of  this  condition.  In 
the  case  of  airborne  early-warning  systems,  penetrating 
aircraft  will  be  able  to  move  in  much  closer  to  the  radar 
aircraft  before  being  detected.  As  an  example,  Davidson 
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Figure  23.  Skip  Effects.  EM  rays  from  the  transmitter 
will  follow  different  paths  depending  on 
the  angle  they  are  transmitted  at  and  the 
surface  duct  characteristics. 


ERef.  5:  pp,  3-10]  discusses  the  effects  noted  by  an  E-2C 
squadron.  One  of  the  phenomena  noted  was  that  for  aircraft 
flying  at  approximately  25,000  feet  ray  angles  of  as  much  as 
5  degrees  below  the  horizontal  were  required  to  penetrate 
solar  heated  subrefractive  conditions,  seriously  degrading 
the  early  warning  ability  of  these  systems.  In  a  similar 
manner,  air  to  ground  communications  can  also  be  affected  by 
subrefractive  conditions. 

Optical  wavelengths  are  affected  by  refraction  in  a 
manner  similar  to  microwaves,  but  on  a  much  reduced  scale 
due  to  the  relative  ineffectiveness  of  water  vapor  gradients. 
Figure  8  shows  that  refraction  does  extend  the  optical 
horizon  beyond  the  geometric  horizon  but  not  nearly  to  the 
extent  that  the  radar/radio  horizon  is  extended.  Since 
there  is  usually  little  measurable  difference  in  the 
refractive  effects  experienced  by  IR  and  visible  wavelengths, 
which  together  form  the  optical  region,  they  will  generally 
be  considered  together  in  this  thesis. 

Refraction  in  a  standard  atmosphere  causes  the  horizon 
seen  by  the  observer  to  appear  higher  and  more  distant  than 
the  geometric  horizon.  McCartney  [Ref.  11:  p.  100]  describes 
a  formula,  given  by  Bowditch,  for  calculating  the  distance 
to  the  optical  horizon  in  a  standard  atmosphere  as 
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where  s  is  the  height  of  the  eye  in  meters.  This  relationship 
is  shown  graphically  in  Figure  24.  Superrefractive  and 
ducting  conditions  will  extend  the  optical  horizon  even  more, 
causing  objects  beyond  and  below  the  horizon  to  become  visible. 
Extremely  strong  inversions  near  the  surface  of  the  earth  will 
form  a  mirage,  known  as  looming,  leading  to  greatly  extended 
optical  horizons.  This  subject  will  be  taken  up  in  more 
detail  in  the  section  dealing  with  mirages.  Figure  24  shows 
that  optical  waves  traveling  in  the  atmosphere  are  also 
subject  to  downward  bending  due  to  refraction,  leading  to 
range  and  angle  errors  in  surveillance  and  tracking  systems. 

The  degree  of  bending  is  dependent  not  only  on  the  refractive 
conditions  prevailing  but  also  on  the  zenith  angle  S .  As 
described  by  McCartney  [Ref.  11;  p.  112],  when  the  zenith 
angle  is  75  degrees  or  less  the  effect  of  refraction  in  a 
standard  atmosphere  is  largely  ignored.  For  larger  values 
of  the  zenith  angle  the  effects  become  apparent  particularly 
if  C  is  greater  ■  than  90  degrees  or  if  superrefraction  or 
ducting  conditions  exist.  When  the  beam  path  lies  at  such 
an  angle,  and  is  also  tangent  to  the  horizon,  the  refractive 
effect  in  bending  and  lengthening  the  path  can  be  significant. 

The  effects  of  refraction  on  millimeter  wave 
propagation  in  both  a  standard  and  non-standard  atmospheres 
is  an  area  not  clearly  understood  at  this  time.  The 
millimeter  wave  region  of  the  EM  spectrum  has  only  recently 
been  actively  considered  for  application  to  communications. 
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surveillances  and  weapons  system’s  guidance.  The  majority  of 
the  information  currently  available  is  largely  theoretical. 

In  a  standard  atmosphere,  there  is  some  reduction  of  signal 
strength  and,  like  microwave  signals,  there  are  errors 
induced  in  range  and  angle  estimation  of  radar  targets. 

These  effects  appear  to  be  rather  small  in  standard  conditions 
T.  Patton  et.al.  [Ref.  12:  p.  170]  states  that  in  non¬ 
standard  atmospheric  conditions  (anomalous  propagation) 
there  is  some  evidence  of  trapping  of  millimeter  waves  by 
surface  ducts ,  but  there  has  been  little  investigation  of  the 
effects  of  elevated  ducts.  Also,  the  importance  of  a  water 
vapor  gradient  in  the  atmosphere  on  millimeter  waves  has  not 
been  adequately  addressed  although  it  does  appear  to  be  more 
significant  in  MMW  propagation  than  in  microwave  propagation 
and  less  so  in  optical  propagation.  D.  Snider,  in  private 

3 

communications  ,  states  that  due  to  the  difference  in  the 
responses  to  humidity  of  optical  and  millimeter  wave  energy, 
if  there  exists,  simultaneously,  a  temperature  inversion, 
created  by  heat  radiation,  and  a  slight  increase  in  humidity 
with  altitude,  optical  rays  would  bend  downward,  due  to  the 
temperature  inversion,  while  MMW  rays  would  bend  upward  due 
to  the  humidity  gradient. 

3 

The  discussion  with  D.  Snider  was  in  regards  to  the 
difference  in  the  effects  of  atmospheric  refraction  on 
optical  and  MMW  systems . 


6,  Methods  of  Predicting  Refractive  Effects 


Calculating  the  refractivity  of  the  atmosphere, 
detecting  the  presence  of  anomalous  propagation  and 
determining  the  effect  of  refractivity  on  EM  wave  propaga¬ 
tion  is  very  complex  and  as  yet  inexact.  Characterizing  the 
variation  of  atmospheric  refractivity  with  altitude  is 
usually  done  by  making  measurements  of  temperature,  pressure, 
and  humidity  as  a  function  of  altitude  using  a  radiosonde 
balloon.  Based  on  these  measurements,  and  utilizing 
Equations  2.4  and  2.5,  a  determination  of  the  refractivity 
and  its  gradient  at  various  elevations  can  be  made .  This 
information  is  generally  sufficient  to  determine  the 
presence  of  anomalous  propagation  phenomena. 

The  use  of  a  radiosonde  for  measurement  has  a  major 
drawback;  it  is  generally  a  slow  response  instrument  and 
somewhat  inaccurate.  Consequently  the  data  collected  by 
the  radiosonde  greatly  limits  the  ability  to  accurately 
calculate  refractive  index  gradients.  Ongoing  efforts  are 
underway  to  develop  and  deploy  smaller  and  more  accurate 
means  of  measuring  refractivity,  directly  or  indirectly. 
Kaelin  [Ref.  13:  p.  124]  in  his  thesis  states  that  examples 
of  new  systems  are  the  mini-refractionsonde  and  airborne 
dropsonde  being  developed  by  the  Naval  Air  Development 
Center,  Warmister,  Pennsylvania.  A  direct  and  accurate,  but 
not  commonly  used  method  of  measuring  atmospheric 
refractivity,  is  the  airborne  microwave  refractometer . 
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A  thorough  and  accurate  table  of  climatological 
records  that  describe  local  geographic  conditions  will  also 
provide  a  valuable  means  of  at  least  probabilistically  cal¬ 
culating  the  expected  refractive  profiles  for  a  region. 
Forecasts  of  anticipated  weather  conditions  can  also  be  used 
to  develop  a  general  idea  of  expected  refractive  conditions. 
These  two  methods,  one  for  prediction  and  the  other  for 
forecasting  of  refractive  conditions  will  both  be  discussed 
further  in  Section  E  of  the  thesis. 

To  determine  the  effect  of  the  existing  refractive 
conditions,  the  frequency  of  the  EM  wave  of  interest  must  be 
known.  The  Navy's  IREPS  will  take  the  condition  of  the 
atmosphere  and  the  wavelength  of  interest  and,  using  a 
desk-top  computer,  graphically  display  the  refractive  index 
calculated  and  the  effect  this  condition  has  on  various 
radar  and  radio  systems  working  in  the  VHF,  UHF,  and 
microwave  region  of  the  EM  spectrum.  Likewise,  the  Navy  is 
developing  a  computer  based  model  to  assess  atmospheric 
effects  (including  refractivity )  on  electro-optical  system's 
performance.  This  system,  described  by  F.  Snyder  [Ref.  1*+], 
is  called  PREOS  (Prediction  of  Performance  Range  for 
Electro-Optical  Systems).  The  Navy  eventually  intends  to 
unite  IREPS  and  PREOS  into  a  single  prediction  system. 

Figure  25  provides  an  example  of  a  PREOS  output  and  Figures 
19  and  20  are  examples  of  IREPS  outputs. 
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Figure  25.  Example  of  PROCAL-PREOS  Output 


Before  methods  can  be  used  to  compensate  for,  correct 
or  take  advantage  of  refractive  effects  on  EM  propagation, 
two  things  are  necessary.  First,  accurate  prediction  and 
forecasting  are  needed.  Some  prediction  means  have  been 
described  here  and  in  a  later  section  prediction  methods  and 
forecasting  methods  will  be  described  in  detail.  Secondly, 
an  in-depth  understanding  of  the  effects  of  refraction  on 
specific  wavelengths  and  even  specific  systems  that  utilize 
propagating  EM  energy  is  needed.  In  the  second  half  of  the 
thesis  several  systems  will  be  used  to  provide  examples  of 
refractive  effects'  wavelength  dependence. 

D.  MIRAGE 

One  form  of  refraction  of  visible  light  that  occurs  in 
the  lower  atmosphere,  with  which  almost  everyone  is 
familiar,  is  the  mirage.  Before  discussing  the  effects  of 
the  several  types  of  mirages  on  EM  energy  transmitted 
through  the  lower  atmosphere,  a  general  description  of  the 
conditions  leading  to  mirage  formation  will  ’■^e  addressed. 

The  majority  of  the  information  used  in  this  discussion 
comes  from  an  article  by  A.  B.  Fraser  and  W.  H.  Mach  in  a 
1980  issue  of  Scientific  American  [Ref.  15:  pp.  29-37]. 

1.  General  Description  of  the  Mirage  Phenomena 


is  largely  dependent  on  temperature,  as  seen  in  Equation  2.5. 
The  greater  the  temperature  gradient,  the  stronger  the 
gradient  of  refraction  and  the  greater  the  amount  of  bending 
of  light.  The  portion  of  the  atmosphere  responsible  for 
forming  mirages  is  located  in  the  first  several  meters  above 
the  surface.  The  exact  form  of  the  mirage  is  highly  variable 
and  is  dependent  on  the  solar  absorption  at  the  surface,  the 
temperature  gradient,  the  position  of  the  object  relative  to 
the  temperature  gradient,  and  the  position  of  the  observer. 

Light  traveling  through  the  shallow  lower  layer  of 
the  atmosphere  will,  due  to  the  existing  temperature 
gradient,  travel  in  a  parabola.  The  curvature  of  the  ray  can 
be  upward  or  downward,  depending  on  whether  there  is  a 
decrease  or  increase  of  temperature  as  altitude  increases. 

A  general  way  to  describe  this  phenomenon  is  to  draw  two 
light  rays  from  an  object,  one  moving  directly  to  the 
observer's  eye  and  the  other  following  a  parabolic  path, 
as  shown  in  the  two-image  inferior  mirage  in  Figure  26. 

The  first  ray  presents  the  light  directly  from  the  object. 

The  second  ray  path  is  bent  upward  (due  to  a  decrease  in 
temperature  at  higher  altitude)  and  thus  appears  to  be  light 
from  an  object  which  is  located  vertically  below  the  actual 
object . 

The  simplest  types  of  mirage  to  distinguish  between 
are  superior  and  inferior.  In  a  superior  mirage  the 
temperature  increases  with  height,  as  in  superrefraction  or 
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Figure  26,  General  Mirage  Effect.  Light  ray  A  from 

object  P  is  above  the  subrefractive  condition 
and  travels  in  a  straight  line  to  the  observer 
C,  Light  ray  B  from  object  P  travels  through 
the  subrefractive  condition  and  thus  creates 
a  virtual  image  P’  which  is  seen  by  the 
observer  C . 


ducting,  and  a  horizontal  surface  will  appear  to  e  concave 
upward,  presenting  an  impression  to  the  observer,  of  being 
inside  a  shallow  bowl.  Figure  27A  is  a  ray  diagram  of  a 
superior  mirage,  showing  that  the  image  of  the  observed 
object  is  displaced  upward.  According  to  McCartney  [Ref.  11 
p.  100],  although  most  mirages  only  involve  relatively  short 
distances  (1-5  km),  a  superior  mirage  created  by  a  very 
strong  inversion  can  make  objec+s  located  far  beyond  and 
below  the  optical  horizon  visible  in  a  phenomenon  known  as 
looming. 

An  inferior  mirage  is  created  by  a  gradient  in  which 
temperature  decreases  with  height  (subrefraction).  As 
described  by  Fraser  and  Mach  [Ref.  15:  p.  35],  a  horizontal 
surface  altered  by  an  inferior  mirage  will  appear  to  be 
convex  upward.  The  impression  is  also  of  a  bowl,  but  this 
time  the  bowl  is  inverted  and  the  observer  is  on  top  of  it. 
There  is  thus  an  optical  horizon  created  beyond  which  the 
surface  being  observed  disappears.  This  phenomenon  is 
referred  to  as  sinking.  The  decrease  of  temperature,  with 
altitude,  -which  leads  to  this  profile  is  common  over 
enclosed  bodies  of  water  in  the  early  morning  and  in  desert 
regions  in  the  late  afternoon.  Note  that,  with  an  inferior 
mirage,  light  rays  from  the  bottom  of  the  object  may  not 
reach  the  observer  if  this  portion  of  the  object  is  beyond 
the  optical  horizon,  as  shown  in  Figure  27b. 
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Figure  27a. 


Superior  Mirage.  The  virtual  image  p*'  is 
seen  by  the  observer  A  as  being  above  the 
actual  object  P. 
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Figure  27b, 


This  time  the  virtual  image  P’  is  sieen  by 
the  observer  A  as  being  below  the  actual 
object  P. 


If  the  temperature  gradient  is  constant  with  height, 
no  distortion  of  the  image  will  occur  in  a  mirage,  since 
displacement  of  all  portions  of  image,  upward  or  downward, 
will  be  the  same.  However,  such  a  situation  is  not  the 
general  case,  since  some  magnification  or  reduction  of  the 
image  will  usually  occur.  When  an  image  is  reduced  in 
height,  the  condition  is  referred  to  as  stooping.  This 
occurs  if  the  temperature  gradient  decreases  as  the  tem¬ 
perature  increases  with  altitude.  Such  a  condition  often 
occurs  over  lakes  or  inlets  on  a  sunny  afternoon,  where 
warm  air  from  the  land  is  carried  over  the  cooler  water.  A 
superior  mirage  is  thus  formed  with  the  bottom  of  the  image 
displaced  upward  more  than  the  top  of  the  image.  The  result 
to  the  observer  is  a  squashed  image  of  the  object,  and  this 
image  will  appear  to  fluctuate  in  vertical  dimension,  as  the 
day  progresses. 

A  magnification  of  the  image,  or  towering  occurs, 
according  to  Fraser  and  Mach  CRef.  15:  p.  35],  in  an  inferior 
mirage  when  the  image  is  displaced  downward ,  but  the  bottom 
of  the  image  is  moved  downward  more  than  its  top  is .  This 
happens  when  the  temperature  gradient  increases  as  the 
temperature  decreases  with  altitude  and  occurs  over  enclosed 
bodies  of  water  in  the  early  morning  or  over  sun-heated 
ground  later  in  the  day. 


The  familiar  mirage  of  the  desert  is  the  two-image 
inferior  mirage.  The  "water”  seen  by  the  observer  is 
actually  an  inverted  image  of  the  sky  seen  below  the  horizon 
The  same  atmospheric  phenomena  that  create  the  inferior 
mirage  and  towering  create  the  two-image  inferior  mirage, 
i.e.,  the  temperature  and  temperature  gradient  both  decrease 
with  height.  The  difference  is  that  the  temperature  profile 
must  be  stronger,  thus  resulting  in  the  ray  being  bent  more 
significantly.  Since  the  object  rises  so  far  above  the 
surface  that  it  is  near  the  top  of  the  temperature  gradient, 
rays  will  travel  both  directly  from  the  object  and  in  a 
parabolic  manner,  similar  to  the  simplified  diagram  in 
Figure  26.  Because  the  ray  of  light  passing  through  the 
f  region  of  strong  temperature  gradient  is  so  strongly  bent, 

it  will  not  join  the  eye  with  the  bottom  of  an  object  but 
will  join  the  eye  with  the  top  of  an  object  to  provide  a 
second,  but  inverted,  image.  As  Fraser  and  Mach  [Ref.  15: 
p.  35]  state: 

"The  image  is  inverted  because  as  the  observer  lifts  his 
gaze  slightly  he  is  looking  through  a  region  of  the 
atmosphere  that  has  a  weaker  temperature  gradient,  so 
that  the  ray  is  less  strongly  curved.  It  will  therefore 
join  the  eye  to  a  point  lower  on  the  object  rather  than 
higher,  as  would  usually  be  expected.” 

An  object  moving  away  from  the  observer  will  first  appear  as 

a  single  image,  then  as  a  double  image,  and  then  vanish 

from  the  bottom  up.  The  presence  of  a  double  image  will 

make  the  ground  between  the  observer  and  tte  object  appear 
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to  be  a  water  surface.  The  distance  to  the  point  where  the 
"water”  begins  is  the  location  of  the  optical  horizon.  As 
the  temperature  cools,  the  optical  horizon  is  extended  and 
the  "water"  recedes  and,  finally,  vanishes. 

According  to  Fraser  and  Mach  [Ref.  15:  p.  36],  the 
most  common  form  of  the  three-image  mirage  is  caused  by 
temperature  increasing  with  height,  for  a  distance,  and  then 
decreasing.  This  condition  will  occur  over  enclosed  bodies 
of  water  on  warm  afternoons,  and  results  in  the  surface  of 
the  water  appearing  as  a  "large,  flattened  letter  S"  -  thus, 
the  surface  appears  to  fold  over.  If  the  inflection  point 
is  gentle,  the  surface  will  appear  to  rise  up  rather  than 
fold  over.  An  object  at  the  distance  where  this  occurs  will 
be  magnified,  creating  an  example  of  the  "Fata  Morgana". 

Fraser  and  Mach  [Ref.  15:  p.  36]  feel  that  the 
occurrence  of  gravity  waves  in  the  atmosphere,  along  with  the 
conditions  leading  to  a  three-image  mirage,  can  result  in 
the  creation  of  the  classic  "Fata  Morgana",  which  will 
sometimes  appear  as  a  hugh  wall  rising  up  above  the  surface 
or  as  a  mountain  range  in  the  distance.  This  combination  of 
stratified  atmosphere  and  gravity  waves  causes  a  periodic 
oscillation  of  the  refractive  index.  Energy  imparted  from 
surface  winds  and  a  countering  force  imparted  by  the  earth’s 
gravity  are  responsible  for  the  waves.  The  oscillating 
bands  of  varying  refractive  indices  will  lead  to  the 
appearance  of  bright  pilasters  and  dark  windows  in  the 
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wall  created  by  the  "Fata  Morgana”,  and  are  actually  the 
image  of  the  surface  below  the  mirage. 

Although  there  are  other  less  common  forms  of  mirage , 
those  described  above  show  some  of  the  great  variety  of 
types  and  the  atmospheric  conditions  leading  to  their 
creation.  Such  conditions  can  potentially  disrupt  the 
operation  of  optical  systems. 

2 .  Effects  of  Mirages  on  Optical  Devices 

As  previously  stated,  there  are  a  variety  of  forms 
of  mirage,  all  dependent  on  the  conditions  of  the  atmosphere 
near  the  surface  of  the  earth.  The  resulting  refractive 
conditions  are  similar  to  those  of  anomalous  propagation 
but  are  on  a  much  smaller  scale,  though  the  line  separating 
anomalous  propagation  from  mirage  is  often  unclear.  The 
basic  difference  between  these  two  forms  of  refraction  are 
the  local  nature  of  the  atmospheric  conditions  leading  to 
mirage  formation  and  the  larger  scale  atmospheric  conditions 
leading  to  anomalous  propagation. 

The  effects  of  mirages  are  most  apparent  on  optical 
systems.  Generally,  there  is  little  difference  in  the 
effects  of  mirages  on  IR  and  visible  wavelengths,  even  in 
the  presence  of  strong  refractive  conditions.  The  effects 
on  optical  tracking  systems  can  be  quite  significant, 
resulting  in  extension  (superior  mirage)  or  shortening 
(inferior  mirage)  of  the  optical  horizon  and  distortion  of 
the  image  (towering  and  sinking).  Mirages  such  as  the  three 


77 


image  mirage  and  the  ’’Fata  Morgana”  can  appear  very  solid 
and  real,  and,  therefore,  be  quite  misleading.  Just  as 
light  from  the  object  to  the  observer  is  strongly  bent,  so, 
too,  will  be  any  visible  energy  directed  from  the  observer 
to  a  target  in  the  mirage. 

The  effects  of  mirages  can  be  significant  for  optical 
trackers,  visible  surveillance,  and  laser  systems  operating 
at  low  elevation  angles  within  the  mirage. 

In  a  similar  manner,  optical  systems  using  the  IR 
spectrum  will  also  be  affected  by  refraction  in  mirage 
conditions.  The  PREOS  model  [Ref.  14]  incorporates  the 
ability  to  predict  the  amount  of  image  distortion  due  to 
mirage  effects. 

The  effect  of  mirages  on  radio  wave  and  microwave 
systems  will  depend  not  only  on  the  general  type  of  mirage 
(inferior  or  superior),  but  also  on  the  strength  of  the 
refractivity  gradient  and  the  mirage  thickness.  The  effects 
will  be  the  same  as  those  predicted  for  surface  based 
anomalous  propagation,  whether  superrefractive ,  ducting,  or 
subrefr active.  The  effect  on  microwave  systems,  such  as 
radar,  will  only  occur  at  very  low  elevation  angles  and  at 
frequencies  high  enough  to  be  above  f^,  (Equation  2.9). 

Prediction  of  mirage  occurrence  is  difficult,  since 
it  is  dependent  on  terrain,  large  scale  atmospheric,  and 
local  atmospheric  conditions.  Mirages  are  generally 
associated  with  partially  or  fully  enclosed  large  bodies 
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of  water,  such  as  lakes,  sounds,  or  bays,  and  with  hot  desert 
regions.  The  probability  of  occurrence  of  mirage  conditions 
and  examples  of  the  effects  of  mirages  on  specific  systems 
will  be  addressed  later. 

E .  TURBULENCE 

Anomalous  propagation  and  mirage  conditions  are  created 
by  relatively  large  scale,  steady  state  variations  in  the 
refractive  index  of  the  atmosphere.  Suck  refractive 
conditions  affect  EM  signals  at  both  the  long  and  short  end 
of  the  spectrum,  from  radio  frequencies  to  visible  light. 

The  air  is,  however,  never  entirely  motionless.  Thus,  much 
smaller  scale  variations  in  the  index  of  refraction  also 
exist.  These  random  fluctuations  of  the  index  of  refraction 
of  the  air  are  created  by  turbulent  motion  in  the  atmosphere 
coupled  with  temperature  and  humidity  gradients.  This 
turbulent  motion  includes  inhomogeneities  in  the  atmosphere 
varying  in  size  from  hundreds  of  meters  in  extent  of  those 
on  the  order  of  1  mm  in  length.  The  effects  of  turbulence 
are  largely  confined  to  the  shorter  wavelengths  (MMW,  IR, 
and  visible),  due  to  the  size  of  the  refractive 
inhomogeneities  relative  to  the  size  of  the  wavelengths. 

1 .  General  Description  of  Turbulence 

The  inhomogeneities  in  the  atmosphere  consist  of 
varying  sized  regions  (discrete  cells)  which  fluctuate  in 
temperature  and  water  vapor  content  in  time  and  in  space. 
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According  to  Zuev  [Ref.  16:  p.  348],  variation  of  1  degree 

centigrade  will  result  in  a  change  in  the  refractive  index 

*•6 

of  approximately  1  x  10"  .  Fluctuations  of  the  temperature 

of  the  air  at  a  given  location  may  be  as  high  as  tenths  of 

a  degree  and  can  vary  over  periods  of  milliseconds  to  seconds. 

Zuev  further  states  that  the  variation  of  temperature  over  a 

horizontal  path  can  be  several  degrees  for  points  on  the 
2  3 

order  of  10  to  10  m  apart. 

These  inhomogeneities  are  themselves  constantly 
changing  in  size  and  location.  They  are  created  and  modified 
by  wind  shear  and  convection  arising  from  vertical  temperature 
gradients  and  from  vertical  water  vapor  gradients  over  lakes 
and  oceans.  As  stated  by  F.  Lutomirski  et.  al.  [Ref.  17: 
p.  1] 

’’The  most  important  factors  affecting  daily  variations  (in 
the  index  of  refraction)  are  short  wave  solar  radiation  and 
long-wave  terrestial  and  atmospheric  radiation,  the  mean 
wind  speed,  type  of  soil  and  ground  cover,  and  ground 
surface  roughness.” 

The  overall  effect  of  the  inhomogeneities  present  in 
a  turbulent  atmosphere  is  to  modify  the  characteristics  of 
optical  propagation  through  the  atmosphere.  The  nature  of 
the  refractive  effect  is  dependent  on  the  operating  wave¬ 
length,  pathlength,  and  nature  of  the  turbulence.  Specifi¬ 
cally,  the  size  of  the  inhomogeneities  is  important  to  the 
overall  effect  of  turbulence  on  optical  and  MMW  transmissions. 
The  size  of  the  inhomogeneities  created  by  wind  shear  and 


convection  will  vary  greatly  from  those  hundreds  of  meters 


in  size  (outer  turbulent  scale)  L  ,  to  those, 

o 

to  Zeuv  [Ref.  16:  p.  349],  on  the  order  of 


according 


(2.11) 


where  1^  is  the  inner  turbulent  scale,  v  is  the  kinematic 
viscosity,  and  e  is  the  rate  of  dissipation  of  turbulent 
kinetic  energy.  The  inner  scale  1^,  is  the  smallest  size 
inhomogeneity  that  can  exist  due  to  viscosity  while  is 
the  largest  such  region  that  can  exist  in  which  temperature 
can  be  expected  to  remain  nearly  constant.  According  to  A. 
Cooper  [Ref.  18:  p.  A-40],  between  these  two  scales  there 
often  exists  a  region  of  well  established  turbulence  which 
can  be  described  by  the  Koloraogorov  spectrum  as 

0  (K)=K"^^^^  ,  (2.12) 

n  ’ 


where  K  =  2Tr/l  is  the  "turbulence  wave  number”  representing 
the  size  of  the  eddy  or  inhomogeneity  and  9j^(K)  is  the 
"intertial  subrange”  of  well  established  turbulence. 

The  structure  function  of  the  refractive  index  between 
two  points  (r  and  r^)  of  a  medium  with  random  inhomogeneities 
is 


=  [n(r,  +  r)  -  n(r, ) 3 


D  (r) 


2 


(2.13^ 


where  r  is  the  distance  between  the  two  points.  In  a 
similar  manner  Zuev  [Ref.  16:  p.  349]  defines  the  temperature 
structure  function 


( r )  = 


(2.14) 


Atmospheric  turbulence  is  generally  described  by  the 
Kolomogorov  model  for  fluid  velocity  turbulence.  For  the 
Kolomogorov  turbulence  model  the  two  Equations  2.13  and  2.14 
can  be  written  equivalently  for  small  scale  turbulence,  as 


D^(r)  =  r^/^ 


(2.15) 


D  (r)  =  C  ^ 
n  n 


(2.16) 


2  2 
where  C  is  the  refractive  index  structure  value  and  C_  is 
n  ^ 

the  temperature  structure  value.  These  relationships  in 
Equations  2.15  and  2*16  only  hold  for  limited  separation 
values ,  when 


1  <<  r  <<  L 

o  o 


(2.17) 


The  two  structure  constants  are  related  by 


C  ^  =  (79P/T^xlO"®)^  , 

n  T 


(2.18) 
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and  are  the  usual  means  of  describing  the  turbulence  level. 

2 

According  to  F.  Hall  [Ref.  19:  p.  TuCl-1]  can  also  be 

described  as  the  difference  in  refractive  index  between  two 
points,  squared,  averaged  and  divided  by  the  2/3  power  of 


the  distance  between  them,  or 


^2  ,  ^2  ,  2/3 

C  =  (n,  -  n„)  /r 
n  12 


(2.19) 


where  r  is  the  distance  between  the  two  points.  Note  that 
2 

has  the  dimensions  of  length  to  the  -2/3.  Generally 

2 

IS  the  only  parameter  needed  to  describe  small  scale 

turbulence  and,  for  this  reason,  is  the  parameter  most  often 

used  in  studying  optical  turbulence  effects.  Actual 

2 

measurements  of  can  be  made  by  determining  the  difference 

in  temperature  between  two  sensitive,  fast  response  thermo¬ 
meters  located  some  distance  apart  or  from  the  temperature 

power  spectrum  from  one  thermometer. 

2 

Values  of  in  the  lower  atmosphere  generally  lie 

—16  —1 2 

between  lO”  cm  and  10  cm  and  vary  widely  with  time  and 
2 

location.  C  values,  and  therefore  turbulence  conditions, 
n  ’ 

-15  -13 

are  considered  to  be  light  around  10  and  severe  around  10 

2 

Variation  in  the  value  of  C  occurs  not  only  along  a  hori- 

n 

zontal  path,  but  also  along  a  vertical  path  from  the  earth’s 
surface  to  the  top  of  an  inversion.  Turbulence  will  generally 


be  strongest  at  the  surface  and  at  the  top  of  an  inversion. 

2 

For  such  a  vertical  path  normally  decreases  in  value  as 
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altitude  increases  until  the  inversion  region  is  reached 
2 

where  values  increase  dramatically.  Above  the  inversion 

there  can  also  be  significant  turbulence  depending  on  the 

wind  and  large  scale  weather  conditions. 

The  usual  result  of  the  many  factors  affecting 

turbulence  formation  is  the  turbulence  maximum  in  the  early 

or  midafternoon  when  heating  of  the  surface  is  greatest  and 

winds  are  light,  and  minima  at  sunrise  and  sunset  when  little 

heating  occurs.  Turbulence  will  occur  at  night,  particularly 

if  an  inversion  develops,  but  only  if  a  light  wind  is  present 

and  will  generally  be  weaker  than  during  the  day.  Thus, 

turbulence  shows  a  characteristic  diurnal  variation  as  shown 

in  Figure  28,  although  the  diurnal  variations  will  be  weaker 

when  the  sky  is  overcast.  Table  3  provides  a  rule-of-thumb 

for  predicting  the  occurrence  of  turbulence  given  existing 

meteorological  and  surface  conditions.  Figure  29  through  31 

2 

are  actual  plots  of  C  values  for  three  different  climatic 

n 

conditions.  Figure  29  is  a  desert  region  in  summer.  Figure 
30  is  a  treeless  area  over  snow,  and  Figure  31  is  a  wet 
grassland  typical  of  many  regions  in  Central  Europe.  Note 
the  distinct  diurnal  cycle  in  each  plot  and  that  the 
strongest  turbulence  occurs  over  the  desert  region. 
Significant  turbulence  is  recorded  in  the  other  two  regions 
as  well,  particularly  over  the  wet  grassland. 


Diurnal  Variation  of  Turbulence ,  The  graph 
shows  the  measured  values  (solid  lines)  of 
optical  turbulence  at  1.5  meters  above  a 
grassy  surface  in  a  semi-arid  region  and 
those  from  bulk  meteorological  measurements 
(dashed  lines). 


TABLE  3 


INTERIM  RULES-OF-THUMB 


Minimum  thermal  turbulence  in  the  boundary  layer  is 
observed  with  the  following  atmospheric  and  earth-surface 
conditions : 


o  Under  a  heavy  overcast,  day  or  night. 

o  Over  water  or  wet  ground,  day  or  night. 

o  With  winds  exceeding  *  8  m  sec”^  (*  15  kn),  nearly 

regardless  of  surface  heating  or  cooling,  day  or  night. 

o  During  clear  weather,  for  periods  of  10  to  30  minutes 
shortly  after  sunrise  and  before  sunset. 

o  During  generally  unstable  days ,  for  brief  periods  in 
the  shadows  of  transient  clouds  (the  same  effect  maybe 
obtained  duriiTg  nights  that  are  generally  stable). 

o  In  descending  volumes  of  air,  specifically  between 

convective  plumes  in  daytime,  free-convection  conditions 

o  On  clear  nights,  with  absolutely  still  air  (any  air 
motion  can  drastically  change  this  situation) . 

Maximum  thermal  turbulence  is  generally  found  with  the 

conditions: 


o  Over  relatively  dry  ground,  with  strong  solar  heating 
and  light  winds . 

o  On  clear  nights  with  light  but  nonzero  winds  (the  most 
unpredictable,  sporadic  turbulence  is  found  under  these 
conditions, 

o  In  the  upper  atmosphere,  in  stable  strata  associated 
with  vertical  wind  shears. 
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Figure  29.  C  va  mean  solar  time  at  .5,  2,  8,  and  16  meters  (top  to 
bSttom)  for  a  dry  desert  in  summer. 
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Figure  31.  C  vs  mean  solar  time  at  .5,  2,  8,  and  16  meters  (top  to 
bottom)  for  a  wet  green  short-grass  grassland. 


.  Principal  Effects  of  Turbulence  on  Optical  Propagation 


Turbulence  in  the  atmosphere,  along  the  optical  path, 

affects  the  optical  properties  of  focused  beams  and  imaging 

2 

systems.  Point  measurements  of  along  an  optical  path 

often  show  wide  variations  in  values  due  to  temperature  and 

water  vapor  variations.  Due  to  this  spatial  variation  in 
2 

,  a  path  integrated  value  is  usually  determined  for  the 

overall  path.  Often  in  calculations  of  the  effects  of 

2 

turbulence  on  optical  systems  the  value  of  used  will  be 

weighted  for  the  calculation.  These  weighting  functions  will 

emphasize  the  effects  of  certain  parts  of  the  path  (the 

center  or  ends )  and  may  totally  ignore  other  regions .  The 

value  of  weighting  will  vary  with  the  type  of  turbulence 

effect  being  considered,  thus,  care  must  be  exercised  in 

2 

using  "weighted”  values  of  in  calculations. 

Instantaneous  deviations  in  n  due  to  turbulence  can 
lead  to  a  number  of  effects  on  optical  systems  to  include: 
scintillation,  wave  front  tilt  or  fluctuations  in  angle  of 
arrival,  beam  spread,  beam  steering  or  beam  wander,  and 
degradation  of  the  Modulation  Transfer  Function  (MTF)  or 
loss  of  resolution.  Each  of  these  effects  will  be  considered 
separately,  although  it  will  be  seen  that  many  are  closely 
related.  E.  Crittenden  [Ref.  20;  p.  2]  states  that  H.T.  Yura 
in  his  analysis  of  turbulence  effects  shows  that  the  effects 
on  an  active  focused  beam  system  will  be  applicable  to  those 


-■•v-;.'  imaging  systems  as  well,  although  the  manifestation  of 

the  effect  may  be  somewhat  different, 
a.  Scintillation 

Scintillation  or  beam  breakup  is  due  to  the 
fluctuation  of  the  intensity  of  a  focused  beam  or  image 
which  has  propagated  through  a  turbulent  medium.  This 
characteristic  sometimes  referred  to  as  "image  boil"  is 
caused  by  random  bending  of  optical  rays  resulting  in 
destructive  and  constructive  interference  patterns.  These 
patterns  vary  with  position  and  time.  Figure  32  shows  an 
example  of  scintillation  of  a  laser  beam.  The  twinkling  of 
stars  viewed  from  the  earth  is  another  example  of 
scintillation. 

The  variance  of  the  log  amplitude  intensity  of 

a  laser  beam,  as  seen  through  a  pinhole  at  the  receiver,  can 

2 

be  related  to  ,  the  propagating  wavelength,  and  the 
pathlength,  according  to  E.C.  Crittenden  et.  al.  [Ref.  20: 
p.  16]  by 


2 
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(2.20) 


where  k  =  lTtl\  and  is  the  propagation  pathlength.  As 
expected,  the  greater  the  level  of  turbulence  or  the  longer 
the  pathlength,  the  greater  the  scintillation  observed.  It 
should  also  be  noted  that  the  shorter  the  optical  wavelength 
■5?^  the  greater  the  scintillation  effect.  Equation  2.20  is  only 
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accurate  for  small  amplitude  fluctuations,  for  8 . 

For  higher  turbulent  levels.  Equation  2.20  is  used  only  to 

2 

infer  the  approximate  relationships  between  ,  pathlength, 
wavelength,  and  log  amplitude  variance. 


Turbulence  levels  generally  vary  along  the 
propagation  path,  rather  than  being  homogeneous  as  implied 
by  Equation  2.20.  This  fact  must  be  considered  in  calculating 
the  actual  log  amplitude  variance,  and  Equation  2.20  can  be 
rewritten  (for  spherical  waves)  as 


2 

where  w  is  a  weighing  function  that  emphasizes  the 

values  in  the  middle  of  the  optical  path  and  ignores  those 

values  at  the  path  ends. 

2 

C  values  can  also  be  calculated  from  measure- 
n 

ments  of  scintillation.  According  to  Crittenden  [Ref.  20: 
p.  16],  however,  in  cases  of  nonuniform  turbulence  levels 
along  the  path  of  propagation  such  calculations  may  not  be 
of  any  use  in  calculating  other  turbulent  effects  because 
of  the  different  weighting  functions  involved, 
b.  Wave  Front  Tilt  and  Beam  Spread 

Wave  front  tilt  or  fluctuations  in  the  angle  of 


arrival  of  a  beam  or  image  are  caused  by  bending  of  the 
entire  wave  front  of  an  image  or  beam  as  it  passes  from  one 
inhomogeneity  to  another.  As  would  be  expected,  the  inner 


scale  turbulent  inhomogeneities  have  little  effect  on 
creating  wave  front  tilt,  but  rather,  it  is  the  larger, 
outer  scale,  turbulence  regions  that  lead  to  this  phenomenon 
The  result  of  this  effect  is  to  cause  "image  dance",  a  move¬ 
ment  of  the  apparent  position  of  an  object’s  location  in  an 
imaging  system.  This  is  the  phenomenon  that  leads  to  the 
false  apparent  location  of  stars  viewed  from  the  earth’s 
surface. 

One  characteristic  of  beam  spread  is  the 
increasing  of  the  width  of  an  initially  narrow  beam  as  it 
propagates  through  the  atmosphere.  This  expansion  of  the 
beam  width  and  resultant  decrease  in  intensity  is  due  to  the 
random  refraction  of  the  beam  as  it  passes  through  a  turbu¬ 
lent  atmosphere.  In  essence,  relatively  small  scale  wave 
front  tilt  is  occurrinF  ^  many  different  locations  across 
the  beam  profile.  Long  term  beam  spreading  can  be  separated 
into  two  components:  beam  expansion  as  described  above  and 
beam  wander.  Beam  wander  is  the  movement  of  a  beam  from  a 
straight  path  due  to  turbulent  inhomogeneities  that  are  much 
greater  in  length  than  the  beam  width,  thus,  refracting  the 
entire  beam.  As  expected,  beam  wander  is  a  function  of  the 
beam  diameter  and  the  turbulence  scale.  Lutomirski  [Ref.  17 
p.  115]  states  that  beam  wander  is  closely  related  to  wave 
front  tilt  and  in  cases  where  the  only  significant  effect  on 
a  beam  propagating  through  the  atmosphere  is  wave  front  tilt 
the  instantaneous  beam  wander  angle  and  wave  front  tilting 
angle  will  be  equal. 


94 


c.  Resolution 

A  cumulative  effect  of  the  various  turbulence 
effects  is  an  overall  degradation  of  the  image  resolution  of 
an  object  seen  through  the  atmosphere.  Characteristics  of 
an  object  are  said  to  be  resolved  if  their  images  are 
discernible  as  being  separate.  This  loss  of  resolution  is 
characterized  by  a  degradation  of  the  Modulation  Transfer 
Function  (MTF)  of  the  atmosphere  by  turbulence.  The  MTF 
describes  the  detail  of  an  image,  transmitted  from  one  point 
to  another,  in  terms  of  the  spatial  frequency.  The  MTF  for 
a  total  optical  system  can  be  calculated  by  determining  the 
MTF  for  the  components  (lens,  mirrors,  and  atmosphere)  and 
forming  the  product  of  these  component  values .  Knowledge 
of  the  MTF  of  the  atmosphere  is  important  in  predicting  the 
effects  of  turbulence  on  low,  medium  and  high  power  laser 
beams  leading  to  beam  wander  and  modification  of  the  beam 
profile.  However,  in  high  power  lasers  the  heating  of  the 
air  must  also  be  considered.  To  be  entirely  accurate,  when 
discussing  the  MTF  as  it  relates  to  optical  svstems,  it  is 
the  Optical  Transfer  Function  (OTF)  that  is  discussed.  The 
OTF  describes  both  the  phase  and  amplitude  changes  in  the 
atmosphere.  When  only  the  amplitude  changes  are  considered, 
the  OTF  is  called  the  MTF. 

Knowledge  of  the  MTF  can  also  aid  in  predicting 
the  turbulence  contributions  to  resolving  ability,  image 
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broadening,  and  image  wander  of  imaging  systems.  Using 

weighting  functions  that  emphasize  the  ends  of  the  path,  MTF 

2 

curves  can  be  developed  from  measured  values  of  and  X. 

For  imaging  systems,  the  Mutual  Coherence  Function 
(MCF)  is  the  equivalent  of  the  OTF  and,  in  fact,  the  atmos¬ 
pheric  MCF  is  the  same  function  as  the  OTF  but  expressed  in 
terms  of  different  variables.  The  MCF  is  of  considerable 
concern  in  coherent  detection  systems.  It  reflects  the  loss 
of  transverse  coherence  of  an  initially  coherent  wave  as  it 
passes  through  the  turbulent  atmosphere.  Knowledge  of  the 
MCF  can  also  be  used  to  provide  information  on  atmospheric 
effects  on  laser  propagation. 

The  OTF  for  a  propagation  system,  along  a  path 
where  the  turbulence  varies  is  described  by  Crittenden 
[Ref.  10:  p.  14]  by 

M(f)  =  exp{-57.64  x"^^^  f^^^  /  °  C  ^  (z/z  dz}  , 

Q  n  o 

(2.22) 

where  z^  is  range,  f  is  the  angular  spatial  frequency, 
defined  as  f  =  1/L,  where  L  is  the  spatial  period  of  the 
distribution  of  the  modulation  of  the  amplitude  (intensity) 
of  the  optical  beam,  and  X is  the  wavelength.  The  same 
function  M(f)  is  used  to  calculate  the  MCF  of  an  imaging 
system  but  using  different  variables. 


The  most  common  parameter  used  to  describe  FLIR 
propagation  through  the  atmosphere  is  the  Minimum  Resolvable 
Temperature  Difference  (MRTD) .  The  MRTD  implies  the  use  of  a 
standard  four-bar  pattern  with  an  individual  bar  length  to 
bar  width  ratio  of  1:7,  The  MRTD  is  the  minimum  temperature 
difference  between  the  bars  and  the  background  for  which  the 
bar  pattern  can  be  resolved  at  the  50%  probability  of  detec¬ 
tion  level  for  a  given  FLIR.  Figure  33  shows  a  four-bar  MRTD 
test  pattern  such  as  that  discussed  above.  As  with  other 
means  of  measuring  resolution,  the  MRTD  for  a  FLIR  system 
can  be  greatly  degraded  by  turbulence . 

Figure  34  provides  a  comparison  of  line  spread 
functions  measuring  the  atmospheric  OTF  for  two  different 
levels  of  turbulence.  Note  that  several  of  the  effects  of 
turbulence  are  occur ing  simultaneously.  The  cumulative 
effect  of  scintillation,  image  broadening,  and  image  wander 
(wave  front  tilt)  is  to  reduce  the  resolution  of  an  image  or 

degrade  the  shape  and  move  the  end  point  location  of  a  beam. 

2 

As  mentioned  earlier,  C  can  be  calculated  from 

n 

measurements  of  scintillation.  Likewise,  Crittenden  et.  al. 

[Ref.  21:  pp.  13-1  to  13-2]  states  that  measurements  of  the 

2 

OTF  and  beam  wander  can  also  be  used  to  calculate  C  and 

n 

with  more  ease  than  using  scintillation. 


Figure  33.  Standard  Four-bar  MRTD  test  pattern  where 
the  bars  represent  a  blackbody  source  of 
temperature  T+AT  and  the  spaces  represent 
the  background  temperature  T, and  AT  is  the 
difference  between  the  background  temperature 
and  the  targets  temperature. 
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Figure  34. 


Line  spread  spectrum  for  two  levels  of 
turbulence.  These  photos  show  that  associated 
with  an  increase  of  C_^  is  a  broadening  of  the 
image  related  to  the  OTF,  increased  detailed 
structure  (scintillation)  and  an  increase  in 
the  Winder  of  the  center  (image  wander). 


3 ,  Manifestation  of  Turbulence  in  Optical  Systems 
In  beam  forming  optical  systems,  such  as  laser 
systems,  scintillation  may  result  in  a  few  dB  to  upwards  of 
40  dB  fluctuation  above  and  below  the  average  irradiance 
value.  For  a  laser  communications  system  this  effect  can 
be  significant  in  degrading  the  bit  error  rate.  Scintillation 
is  less  important  at  IR  and  longer  wavelengths  than  at 
visible  wavelengths.  According  to  Lutomirski  [Ref.  17:  p.  230] 
recent  studies  have  shown  that  at  the  longer  wavelengths  (IR) 
a  coherent  detection  laser  communications  system  is  preferable, 
while  at  visible  wavelengths,  an  incoherent  detection  system 
is  best  when  turbulent  effects  are  concerned.  Beam  wander  can 
cause  significant  enough  movement  of  a  laser  beam  to  degrade 
detection  in  an  optical  communications  system  or  to  sweep  the 
beam  off  the  target  in  a  laser  designation  or  laser  weapon 
system.  Wave  front  tilt  will  also  cause  a  degradation  of 
detection  in  an  optical  receiver,  although  this  may  be  the 
simplest  turbulence  effect  to  correct.  Beam  spread  can 
result  in  a  reduction  of  the  irradiance  of  a  laser  weapon  and 
a  sensitivity  reduction  in  laser  communications  systems.  The 
MCF  can  be  used  to  determine  the  limiting  mean  irradiance 
from  an  initially  coherent  beam  and  the  turbulence  limited 
S/N  ratio  of  an  optical  heterodyne  (coherent)  receiver. 

Table  4  consists  of  two  tables  listing  several  types  of  army 
laser  systems  and  the  types  of  turbulence  induced  degradation 


each  system  can  suffer  from. 


TABLE 
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Turbulence  will  create  the  same  effects  on  imaging 
systems  but  will  often  be  manifested  in  a  different  manner 
than  the  effects  on  beam  producing  optical  systems. 
Scintillation  will  create  variations  in  the  intensity  of  the 
image  of  the  object  being  viewed  and  will  thereby  degrade 
detection  and  identification.  Image  wander  can  cause  a 
degradation  of  detection  ability,  dependent  on  the  type  of 
imaging  system  used,  specifically  the  nature  of  the  detection 
method  used  by  the  imaging  system.  Image  broadening  may 
degrade  both  detection  and  identification.  Wave  front  tilt 
can  present  a  false  apparent  location  of  the  object.  Tur¬ 
bulence  induced  degradation  of  the  resolution  of  an  image 
can  be  quite  severe,  affecting  detection,  discrimination 
and  identification.  The  limiting  resolution  due  to 
turbulence  can  be  determined  from  the  MCF. 

It  is  apparent  from  this  discussion  that  turbulence 
effects  on  optical  systems  are  varied,  but  generally  result 
in  a  degradation  of  the  system's  operation.  Some  means  have 
been  developed  to  compensate  for  turbulence  effects,  a  few 
of  which  will  be  mentioned. 

Lutomirski  [Ref,  17:  p.  65]  states  that  perhaps  the 
easiest  form  of  turbulence  effects  to  correct  is  wave  front 
tilt  to  the  entire  beam  or  image,  if  it  occurs  largely  by 
itself.  Mechanical  or  electronic  tracking  of  the  image  can 
be  used  in  a  system  to  correct  the  displacement  in  the  image 
plane,  producing  a  nondegraded  image. 


other  forms  of  turbulence  effects  are  not  so  easy  to 
correct,  particularly  when  several  occur  at  once.  Because 
of  the  small  scale  of  turbulence  effects,  turbulence 
correction  requires  a  resolution  element  related  to  the  size 
of  the  turbulence  scale.  These  elements  must  be  small  and 
response  times  very  short,  due  to  the  rapidly  changing 
nature  of  turbulence  effects.  This  is  the  basis  of  adaptive 
optics,  the  purpose  of  which  is  to  minimize  turbulent 
effects  on  optical  systems,  both  imaging  and  beam  forming. 

The  adaptive  optics  system  senses  phase  shifts  induced  in 
the  optical  beam  by  turbulence  and  introduces  a  compensating 
phase  shift  to  correct  for  this.  In  particular,  the  Coherent 
Optics  Adaptive  Technique  (COAT)  is  basically  a  self-phasing 
phased  array  system.  Errors  in  the  phase  of  the  incoming 
signal  are  sensed  across  the  wavefront  and  corrected  by  the 
COAT  system.  The  COAT  output  aperature  is  divided  into  many 
small  segments  whose  phase  can  be  controlled  independently 
using  mirrors  or  splitters,  providing  the  means  to  correct 
phase  errors  in  the  incoming  beam.  This  technique  has  been 
shown  to  improve  to  some  degree  the  degrading  effects  of 
scintillation,  beam  spread  and  beam  or  image  wander  but  can 
not  as  yet  fully  improve  the  system's  performance  to  that  at 
low  turbulence  conditions. 


.  Turbulence  and  MMW  Propagation  in  the  Atmosphere 

Although  turbulence  has  been  largely  studied  because 
of  its  performance  degradation  of  optical  systems,  it  has 
been  theorized  for  some  time  that  it  may  also  adversely 
affect  millimeter  waves.  Recent  tests  have  confirmed  this 
fact  and  have  also  shown  agreement  with  theory  in  that  there 
is  a  stronger  dependence  on  the  atmosphere ’ s  hum  ity  at  MMW 
than  at  optical  wavelengths.  Thus,  in  calculati  s  involving 
the  effects  of  turbulence  on  MMW  propagation,  be  ,,  the 
temperature  structure  value  and  C^,  the  humidity  structure 
value,  must  be  considered.  Additionally,  the  cross  corre¬ 
lation  of  these  two  parameters,  C^g>  must  be  included  in 
determining  turbulence  effects.  Scintillation  and 
fluctuations  in  the  angle  of  arrival  of  the  beam  in  radar 
applications  have  been  shown  to  be  the  primary  turbulence 
effects  at  millimeter  wavelengths.  McMillan  and  Bohlander 
et.  al.  [Ref,  22:  pp.  32-39],  in  their  1982  article  on  MMW 
propagation,  provide  a  detailed  discussion  of  the  formulas 
used  in  calculating  intensity  fluctuations  (scintillation) 
and  angle  of  arrival  fluctuations. 

As  described  by  McMillan  and  Bohlander  [Ref.  22: 
pp.  33-36],  results  of  tests  at  White  Sands  Missile  Range 
in  New  Mexico  show  that  scintillation  effects,  though  present 
are  generally  not  expected  to  be  significant  on  MMW  radar 
systems.  However,  fluctuations  of  the  angle  of  arrival  of 
MMW  beams  are  on  the  Scime  order  of  expected  angular  accuracy 


and,  thus,  can  be  significant  in  degrading  angle  tracking 
ability  as  shown  in  Figure  35.  The  results  of  the  testing 
seem  to  indicate  that  turbulence  effects  on  angle  tracking 
may  be  smaller  than  those  due  to  overall  anomalous 
propagation  effects.  Although  more  data  on  this  subject  is 
necessary  before  the  full  nature  of  turbulence  effects  on 
MMW  is  known,  the  results  to  date  appear  to  indicate  that 
turbulence  should  be  considered  when  discussing  MMW 
applications . 


F.  PREDICTION  AND  PROBABILITY  OF  OCCURRENCE  OF  ANOMALOUS 

PROPAGATION,  MIRAGE  FORMATIONS,  AND  TURBULENCE 

The  effects  of  atmospheric  refraction  on  EM  wave 
propagation,  regardless  of  severity,  are  of  little  real 
consequence  unless  these  conditions  occur  with  some 
frequency.  Furthermore,  to  respond  to  these  effects,  it 
is  necessary  to  be  able  to  forecast  their  occurrence  and/or 
determine  their  existence.  Some  discussion  on  meteorological 
conditions  and  characteristics  of  terrain  that  lead  to  the 
three  major  forms  of  atmospheric  refraction  has  already  been 
dealt  with  and  will  not  be  reit..rated  here,  except,  as 
needed,  for  clarity.  Rather,  this  section  will  attempt  to 
look  at  the  larger  scale  causes  and  characteristics  of 
atmospheric  refraction  and  focus  on  the  probability  of 
occurrence  of  anomalous  propagation,  mirage  conditions,  and 
turbulence  in  various  regions  of  the  world.  Finally,  two 
specific  regions.  Central  Europe  and  the  Mid-East  will  be 
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Figure  35.  Calculated  peak-to-peak  angle-of-arrival 
fluctuations  for  millimeter  wave  signal. 


looked  at  in  some  detai''  ,  not  only  because  of  the  obvious 
military  interest,  but  also  to  compare  these  two  distinctly 
different  areas. 

Anomalous  propagation  will  be  examined  first,  for  two 
reasons:  first,  because  a  larger  body  of  information  exists 

on  this  subject;  second,  because  much  of  what  is  discussed 
concerning  anomalous  propagation  can  be  directly  related  to 
turbulence  and  mirage  forecasting  and  prediction. 

1.  Anomalous  Propagation 

A  great  deal  of  study  has  been  conducted  towards 
understanding  the  causes  of  anomalous  propagation  in  its 
many  forms,  and  in  determining  analytical  methods  of 
predicting  its  existence.  Most  of  these  have  already  been 
discussed.  Forecasting  anomalous  propagation  conditions  can 
also  be  performed,  given  an  understanding  of  the  location 
and  terrain  features  of  the  region  of  interest,  the  time  of 
year  and  time  of  day,  and  existing  large-scale  weather 
conditions  (air  masses,  high  and  low  pressure  regions,  etc.), 
many  of  which  interact  closely. 

The  general  nature  of  the  earth's  surface,  be  it  a 
large  body  of  water  or  a  large  land  mass,  contributes  to  the 
expected  refractive  characteristics.  This  is  due  to  the 
physical  nature  of  the  surface  (mountains,  desert,  shallow 
tropical  sea,  etc.)  and  the  effect  of  the  surface  on  the  air 
masses  above.  The  persistence  of  superrefractive  and  ducting 


conditions  over  the  oceans  of  the  world  is  due  largely  to  the 


constant  presence  of  moisture  at  low  altitudes  and  the 
relatively  slow  changes  of  large  air  masses  over  these 
bodies  of  water.  This  persistence  over  water  tends  to 
decrease  with  distance  (north  or  south)  from  the  equator. 

Figure  36  shows  this  fact  graphically. 

Over  large  land  masses  there  generally  is  no  such 
consistency  in  atmospheric  conditions  as  exists  over  oceans. 

The  humidity  content  and  temperature  of  the  air,  as  well  as 
the  large  scale  weather  conditions  (fronts  and  pressure 
regions)  change  much  more  rapidly  over  land  than  over  sea. 

For  this  reason,  persistent  ducting  is  much  rarer  over  land, 
and  forecasting  of  anomalous  propagation  is  generally  more 
difficult . 

The  presence  of  certain  terrain  features  also  tends 
to  have  an  effect  on  the  occurrence  and  characteristics  of 
refraction.  High  mountains  will  significantly  affect  the 
movement  of  air  masses  across  the  surface  of  the  earth  and 
impact  on  the  creation  of  trade  wind  inversions  and  land-sea 
breeze  recirculation.  C.  Samson  [Ref.  9:  p.  7]  states  that 
large  river  valleys  may  show  patterns  of  refractivity 
differing  markedly  from  surrounding  regions  due  to 
stratification  of  the  atmosphere  by  air  drainage,  above  the 
valleys.  Heating  of  large  arid  surfaces  will  often  result 
in  subrefractive  conditions. 

Large  scale  weather  conditions  are  the  most  significant 
factor  in  forecasting  overland  anomalous  propagation.  Air 
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(a)  February 


(b)  August 


Figure  36.  Regions  of  persistent  oceanic  atmospheric 
duct  layers. 


a 

mass  characteristics,  the  presence  of  high  pressure  regions, 
and  the  movement  of  fronts  and  cyclones  are  all  important 
in  the  formation  of  anomalous  propagation  conditions. 

Air  masses  are  regions  of  air  in  which  the  properties 
of  the  atmosphere  are  relatively  uniform.  These  air  masses 
are  formed  in  source  regions  which  include  large  desert  areas 
tropical  seas,  and  polar  regions.  The  presence  of  a  high 
pressure  system  over  a  source  region  enhances  the  development 
of  the  air  mass.  Source  regions  used  in  the  classification 
of  air  masses  are  polar  ’P*,  tropical  'T*,  continental  'C*, 
and  maritime  ’M’.  Table  5  describes  the  major  air  masses 
and  their  characteristics.  Each  of  these  air  masses  has  a 
unique  source  region  and  movement  pattern  that  varies 
seasonally.  As  would  be  expected,  each  also  has  a  charac¬ 
teristic  refr.  ^tivity  profile  associated  with  it.  Figure  37 
shows  the  average  profiles  associated  with  these  air  masses. 

The  profiles  shown  in  Figure  37  can  be  modified  by 
the  presence  of  subsidence  associated  with  high  pressure 
regions.  Large,  permanent  high  pressure  regions  exist  over 
the  subtropical  and  mid-latitude  oceans  of  the  world,  to 
some  extent  year-round.  These  high  pressure  regions  shift 
toward  the  poles  during  each  hemisphere's  summer  months  and 
tend  to  drive  the  movement  of  frontal  regions  poleward. 
Frontal  regions  and  cyclones,  which  are  common  in  the  mid¬ 
latitudes  and  polar  regions,  tend  to  destroy  elevated 
layers,  and,  thus,  no  significant  ducting  or  superrefraction 
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TABLE  5 

AIR  MASS  CLASSIFICATIONS  AND  GENERAL  DESCRIPTIONS 


Air  Mass  Source  region(s)  Properties  at  Source 


Polar  maritime 
(Pm) 

Oceans  in  latitudes 
greater  than  50° 
(approx . ) 

Cool,  rather  moist; 
unstable 

Polar  contin¬ 
ental  (Pc) 

Continents  in 
vicinity  of  Artie 
Circle;  Antarctica 

Cold,  dry;  stable 

Arctic  or 
Antarctic  (A) 

The  Arctic  Basin 
and  (central) 
Antarctica  in 
winter 

Very  cold  and  dry; 
very  stable 

Tropical 
maritime  (Tm) 

Sub-tropical  oceans 

Warm;  moist,  and 
rather  unstable 
near  surface,  dry 
and  stable  above 

Tropical 

continental 

(Tc) 

Deserts  in  low 
latitudes:  pri¬ 
marily  the  Sahara 
and  Australian 
deserts ,  but  also 
south-west  U.S.A. 
and  Mexico  in 
summer 

Hot  and  dry; 
unstable 
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Figure  37.  Typical  B-profiles  in  the  source 
for  (A)  Summer  and  (B)  Winter, 
ducting  gradient. 


region 
.  is  the 


is  expected  in  the  vicinity  of  fronts  or  of  traveling  low 
pressure  regions  (cyclones).  Due  to  the  persistence  of 
these  traveling  weather  patterns  in  extreme  northern  or 
southern  (polar)  regions,  anomalous  propagation  is  generally 
a  rare  and  intermittent  occurrence  at  and  near  the  poles. 

Coastal  regions  will  reflect  the  variation  in  humidity 
and  temperature  related  co  the  ebb  and  flow  of  maritime  and 
continental  air  masses  by  exhibiting  extreme  horizontal  and 
vertical  refractive  gradients. 

Over  land,  and  to  a  lesser  extent  over  water,  there 
tends  to  be  a  distinct  seasonal  variation  in  the  occurrence  of 
anomalous  propagation.  This  can  generally  be  tied  to  the 
prevailing  large-scale  weather  conditions  and  to  the 
intensity  of  solar  irradiance. 

Many  of  the  smaller  oceanic  islands  show  little 
seasonal  change  in  refractivity  patterns ,  due  to  the  fact 
chat  their  climates  are  dominated  by  the  maritime  air  masses. 
Islands  large  enough  to  affect  these  air  masses  will,  on  the 
other  hand,  often  develop  distinctly  seasonal  refractive 
gradient  patterns. 

Continental  regions,  particularly  those  in  the 
mid-latitudes,  tend  to  show  significant  seasonal  variation 
in  refractive  characteristics,  with  the  greatest  occurrence 
of  superrefractive  and  ducting  conditions  occurring  in  the 
summer  months  and  early  fall.  This  is,  in  part,  due  to  the 
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shifting  of  the  subtropical  high  pressure  regions,  higher 
mean  temperatures,  and  the  greater  amount  of  moisture 
present  in  air  masses  in  summer. 

Seasonal  shifts  in  the  prevailing  wind  patterns,  such 
as  the  monsoon  circulations,  lead  to  distinct  seasonal 
variations  in  the  refractive  patterns  of  regions  thus  affected 

Many  regions  will  also  show  strong  diurnal  variations 
in  refractivity  gradients  due  to  surface  changes  in  tempera¬ 
ture,  humidity,  and,  in  some  regions,  shifts  in  wind  direction 
In  general,  according  to  Samson  [Ref.  9:  p.  9],  continental 
and  coastal  regions  tend  to  show  anomalous  propagation  as 
occurring  most  often  between  2000  and  0700  LST  with  the 
fewest  occurrences  between  1200  and  1500  LST.  Oceanic 
islands  show  just  the  opposite  trend  in  the  occurrence  of 
anomalous  propagation.  Subrefraction  in  desert  regions  also 
occurs  most  often  in  the  late  afternoon.  As  was  mentioned 
earlier  in  this  thesis,  it  is  not  uncommon  in  semi-arid  and 
desert  regions  to  see  subrefractive  conditions  develop  in  the 
afternoon  and  superrefractive  conditions  form  at  night. 

Figure  38  shows  this  graphically,  along  with  the  resulting 
effects  on  IR  rays.  These  diurnal  characteristics,  which 
tend  to  be  strongest  in  the  warmer  summer  and  early  fall 
months,  are  superimposed  on  the  large  scale  changes  in  air 
masses . 

Forecasting  anomalous  propagation  conditions  involves 
many  variables.  However,  studies  of  meteorological  data 


collected  from  weather  stations  located  worldwide  has  led  ro 
the  development  of  distinct  climatic  types,  which  are 
terrain  types  with  similar  climates  found  across  the  earth's 
surface.  Table  6  is  a  list  of  climatic  types  and  their 
associated  meteorological  and  refractive  characteristics. 

Meteorological  data  collected  at  weather  stations 
has  been  used  to  prepare  climatic  charts  reflecting  the 
percentage  of  time  superrefractive  conditions  or  ducting 
occurs,  on  the  average,  across  the  surface  of  the  earth. 
These  charts  provide  some  idea  as  to  the  probability  of 
occurrence  of  superrefraction,  surface  ducting,  and  elevated 
ducting  in  any  region  of  interest  where  sufficient  data  has 
been  collected,  and  also  reflect  the  seasonal  variation  in 
anomalous  propagation  over  most  of  the  earth's  surface. 
Figure  39  and  40  are  representative  of  the  types  of  graphic 
output  climatic  charts  can  present.  These  figures  show  the 
percent  of  occurrence  of  surface  ducts,  elevated  ducts,  and 
superrefractive  layers.  Note  that  Figure  39  is  of  Northern 
and  Central  Europe  and  Figure  40  is  of  the  Mid-East.  It  is 
readily  apparent  that  the  percent  of  occurrence  of  surface 
layers  and  elevated  layers  varies  a  great  deal  between  the 
two  regions,  and  that  coastal  areas  tend  to  show  a  greater 
occurrence  of  ducting  than  do  inland  regions.  Charts  like 
these,  while  helpful,  do  not  show  the  percentage  of  time 
that  subrefractive  conditions  occur,  nor  do  they  explain 
whether  the  anomalous  propagation  that  they  describe  is 
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39B.  Percent  occurrence  of  elevated  layers  for  Northern 
and  Central  Europe. 


Figure  4 OB.  Percent  occurrence  of  elevated  layers  for  the 
Mediterranean  region,  to  include  the  Nidr-east 


Forecasting  turbulent  conditions  is,  perhaps,  more 

difficult  than  other  forms  of  refraction,  since  the 

interplay  of  surface  roughness,  solar  irradiance  and  wind 

strength  is  difficult  to  predict.  Likewise,  the  direct 

2 

measurement  of  is  not  always  possible,  however,  Wesely 

and  Alcaraz  [Ref.  23:  p.  9]  state  that  an  indirect  method 

2 

has  been  demonstrated.  With  this  method,  C  is  calculated 

n 

using  sensible  and  latent  heat  flux  components  of  the 

surface  energy  budgets,  which  are  available  in  meteorological 

literature.  This  method  is  most  accurate  for  heights  of 

2 

4  meters  or  less,  where  comparison  of  values  of  C  determined 

2 

by  the  indirect  method  and  by  direct  measurements  of  C 

n 

appear  to  agree  closely,  and  thus  provide  a  reasonably 
accurate  means  of  predicting  turbulence  over  land.  Just  as 
with  mirage  conditions  and  subrefraction,  turbulence  condi¬ 
tions  occur  most  often,  and  with  the  greatest  strength,  in 
desert  regions  in  the  late  afternoon,  but,  as  seen  earlier, 
significant  turbulence  does  occur  elsewhere  and  at 
different  times  of  day. 

Although  few  actual  climatic  charts  of  turbulence 
conditions  currently  exist,  over  land  solar  heat  flux  charts 
and  irradiance  charts  are  available,  and  provide  a  fairly 
accurate  idea  of  those  regions  where  significant  over-land 
turbulence  can  be  expected.  Figure  41  is  an  example  of  such 
a  radiation  chart.  Note  that,  in  general,  the  further  away 
from  the  equator,  the  less  solar  radiation  is  received. 


Figure  41, 


Solar  radiation  c^art  of  eastern  hemisphere 
measured  in  kj/cm  . 
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3 .  Comparison  of  Mid-East  and  Central  Europe 

The  Mid-East  and  Central  Europe  are  two  distinctly 
different  regions  in  terms  of  both  terrain  and  climate.  The 
Mid-East  generally  is  typified  by  a  hot,  dry  climate  with 
relatively  little  seasonal  variation  in  temperature  and 
humidity  except  for  its  southeastern  areas,  which  are 
affected  by  monsoons.  Much  of  the  interior  is  rocky  desert 
with  coastal  mountain  ranges  rising  up  along  the  Mediterranean 
Sea.  Coastal  regions  on  theMediterranean  Sea  exhibit  mild 
Mediterranean  climates,  while  similar  regions  on  the  Arabian 
Sea  are  significantly  affected  both  by  continental  air  masses 
and  by  maritime  air  masses. 

Central  Europe  has  a  continental  climate,  characterized 
largely  by  cold  winters  and  cool-to-warm  summers.  There  is 
a  moderate  amount  of  precipitation  year-round  and  relatively 
high  humidity  in  the  summer  months.  Coastal  regions  tend  to 
experience  some  modification  of  the  continental  climate  by 
maritime  air  masses.  The  terrain  varies  from  relatively 
flat,  open  land  in  the  north  to  wooded,  hilly,  and 
mountainous  terrain  to  the  south. 

It  should  come  as  no  surprise,  judging  from  the 
differences  in  climate  and  terrain,  that  the  refractive 
characteristics  of  these  two  regions  are  distinctly 
different,  as  shown  in  Figure  39  and  40.  The  Mid-Eastern 
desert  shows  a  significantly  higher  year-round  incidence  of 
both  surface  and  elevated  superrefractive  conditions. 


Before  going  any  further,  it  must  be  noted  that  all 
measurements  for  the  Mid-East  are  from  coastal  areas,  and, 
thus,  do  not  provide  an  accurate  ina.  ation  of  the 
refractive  conditions  in  the  interior.  Figure  42,  however, 
shows  that  the  incidence  of  superrefraction  and  elevated 
ducting  decreases  with  distance  from  the  coastlines  -  in  the 
Mid-East,  and  in  Europe  as  well.  Figure  42  also  highlights 
the  seasonal  nature  of  ducting  and  superrefraction  by 
showing  the  wide  variations  between  the  average  percent  of 
annual  occurrence  of  refractive  conditions  and  the  percent 
of  occurrence  for  each  month,  given  with  the  maximum 
occurrence  of  ducting  and  superrefraction. 

Figures  39  thru  41  provide  a  good  indication  of 
general  refractive  conditions  for  these  two  regions.  From 
these  figures  it  is  apparent  that  superrefraction  and  ducting 
occur  far  more  commonly  in  the  Mid-East  than  in  Central 
Europe.  Coastal  areas  in  both  regions  show  the  greatest 
incidence  of  these  forms  of  anomalous  propagation,  with  the 
Mid-East  coastline  along  the  Arabian  Sea  showing  an 
extremely  high  year-round  incidence.  Europe  has  a  strong 
seasonal  variation  in  the  percent  of  occurrence  of  these  two 
forms  of  anomalous  propagation,  whereas,  in  the  Mid-East, 
this  occurs  only  along  the  Mediterranean  coast. 

Subrefractive  conditions,  mirage  formation,  and 
optical  turbulence  are  far  more  common  in  the  Mid-East 
because  of  its  hot,  dry  climate.  Some  incidence  of  mirage 


Figure  42B.  Percent  occurrence  of  elevated  ducts  during 
the  month  of  maximum  occurrence.  This  month 
of  maximum  occurrence  varies  from  region  to 


Percent  occurrence  of  superrefractive  layers 
(SRLRs)  during  the  month  of  maximum 
occurrence.  This  month  of  maximum  occurrence 
varies  from  region  to  region. 


formation  is  found  in  Europe,  particularly  over  relatively 
large  bodies  of  water  and  in  the  more  arid  south,  but  only 
during  the  warmer  summer  months.  In  the  Mid-East,  these 
conditions  occur  year-round  and  are  fairly  regular  phenomena 
By  referring  to  Figure  41,  it  can  be  seen  that  the  Mid-East 
receives  a  far  higher  level  of  solar  irradiation  annually 
than  does  Europe,  and  this  provides  the  primary  ingredient 
for  the  formation  of  subrefraction,  mirages  and  optical 
turbulence . 

4 .  Conclusions 

Although  it  can  be  shown,  with  some  accuracy,  just 
what  percentage  of  the  time  superrefraction  and  ducting  can 
be  expected  to  occur  worldwide,  and,  to  a  lesser  extent, 
subrefraction,  mirage  formation,  and  optical  turbulence,  it 
is  difficult  to  determine  just  where  the  occurrence  of  these 
refractive  characteristics  becomes  ’’significant".  Such 
decisions  must  be  based  on  each  specific  weapons  system  (or 
component  of  a  weapons  system)  affected,  and  on  the  nature 
of  the  effects,  as  well  as  on  the  expected  probability  of 
occurrence  of  anomalous  propagation.  Much  research  is  being 
conducted  by  the  U.S.  Army  at  the  Atmospheric  Sciences 
Laboratory  at  White  Sands  Missile  Range  and  at  the  U.S.  Army 
Electronic  Proving  Grounds  (EPG)  at  Fort  Huachuca,  Arizona. 
The  research  concerns  the  refractive  effects  for  specific 
regions  of  the  EM  spectrum,  and  on  specific  weapon  systems. 


What  can  be  determined  by  the  results  of  this  dis¬ 
cussion  is  that  the  various  forms  of  atmospheric  refraction 
are  not  rare,  and  may  even  be  the  norm  in  certain  regions  of 


the  world  at  certain  times  of  the  year.  Due  to  the  wide 
regional  variability  in  refractive  forms  and  in  probability 
of  occurrence,  it  becomes  apparent  that  each  region  of  the 
earth  poses  a  unique  refractive  problem.  For  this  reason, 
when  considering  military  operations  in  a  specific  region, 
the  potential  refractive  conditions  should  not  be  approached 
in  generalities,  but,  rather,  should  be  addressed 
specifically  to  the  region  of  interest. 
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III.  ATMOSPHERIC  RETRACTION  AND  GROUND  WARFARE 

A.  INTRODUCTION 

In  Section  II  of  this  thesis,  the  three  major  forms  of 
atmospheric  refraction  were  discussed.  The  char'cteristics 
of  each  of  these  forms  and  their  effects  on  EM  waves 
propagating  through  the  atmosphere  were  elaborated  on  at 
length  in  sub-Sections  C  through  E.  Those  discussions 
showed  the  wide  diversity  of  both  the  forms  of  refraction 
and  their  effects  across  the  EM  spectrum.  Sub-Section  F 
described  the  factors  affecting  the  ability  to  predict  and 
forecast  atmospheric  refractive  conditions,  contained  a 
general  discussion  of  the  probability  of  occurrence  of  the 
various  forms  of  atmospheric  refraction  across  the  earth’s 
surface,  and  concluded  with  a  description  of  the  conditions 
characteristic  of  the  Mid-East  and  Central  Europe. 

The  purpose  of  this  background  was  to  show  that 
atmospheric  refractive  conditions  which  have  a  measureable 
effect  on  EM  signals  should  not  simply  be  considered  as  so 
rare  or  unpredictable  that  they  can  be  ’’written  off"  as 
inconsequential.  Rather,  it  has  been  shown  that,  at  times, 
these  conditions  are  quite  common,  can  have  a  significant 
impact  on  EM  signal  propagation,  and  are  reasonably 
predictable.  Of  course  the  probability  of  occurrence  and 


ability  to  forecast  the  three  major  forms  of  refraction 
varies  from  region  to  region  (and  often  season  to  season) 
across  the  earth. 

With  this  general  background  as  a  foundation,  the 
remainder  of  this  thesis  will  now  discuss  specific  examples 
of  the  effects  of  atmospheric  refraction.  This  will  be  done 
in  two  ways  -  by  use  of  computer  simulations  and  by 
qualitative  study.  The  examples  will  cover  the  full  EM 
spectrum  and  are  intended  to  enhance  the  reader’s  under¬ 
standing  of  this  subject  as  important  for  the  employment  of 
weapon  systems  utilizing  the  EM  signals  propagating  through 
the  atmosphere. 

B.  EFFECTS  OF  ATMOSPHERIC  REFRACTION  ON  TANK  OPTICAL  FIRE 

CONTROL  SYSTEMS 

.  4 

This  example  is  based  on  private  communications  with 
Dr.  Don  Snider  of  the  Atmospheric  Sciences  Laboratory  (ASL) 
in  White  Sands,  New  Mexico.  This  section  discusses  the 
effects  of  superrefractive  and  subrefractive  conditions  on 
the  accuracy  of  tank  fire.  The  tank’s  optical  fire  control 
system  (using  both  visible  light  and  IR  energy)  and  current 


The  communication  with  Dr.  Snider  concerned  the  effects 
of  atmospheric  refraction  on  tank  optical  fire  control 
accuracy . 


tank  gunnery  techniques  (boresighting,  synchronizing, 
zeroing  of  the  main  gun,  and  firing  at  center-of-mass )  are 
addressed . 

Current  tank  gunnery  doctrine,  as  described  in  FM  17-12 
[Ref.  24:  pp.  4-9  -  4-10],  requires  boresighting  and 
synchronizing  both  sights  and  main  gun,  followed  by  zeroing 
the  tank  weapon  systems  prior  to  firing  the  tank’s  main  gun 
and  machine  gun(s).  Boresighting  establishes  a  convergent 
relationship  between  the  axis  of  the  tube  of  the  main  gun 
and  its  direct  fire  sight.  Synchronization  insures  that 
main  gun  and  sights  remain  aligned  when  they  are  elevated 
and  when  they  are  depressed.  The  tank  crew  is  expected  to 
fire  from  boresight  using  telescope  and  periscope  during  the 
day  and,  if  so  desired,  the  IR  night  vision  sight  at  night. 

At  the  same  time  the  crew  must  compensate  for  drift,  cant, 
parallax,  tube  droop  and  winds  on  the  battlefield.  As  would 
be  expected,  there  are  a  great  many  potential  sources  of 
error  involved  in  tank  gunnery.  The  magnitude  of  the 
sources  '^f  error  will  vary  depending  on  internal  and 
external  conditions  (relative  to  the  tank)  and  on  the  weapon 
and  ammunition  being  fired. 

Figure  43  shows  a  typical  tank  fire  control  error  budget. 
It  is  currently  thought  that  the  most  important  cause  of 
error  in  firing  the  kinetic  energy  (KE)  round  is  round-to- 
round  variation  and  loss  of  zero  or  improper  zero.  For  the 
slower  HEAT  round,  crosswinds  and  air  temperature  extremes 
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Figure  43.  Typical  tank  fire  error  budget.  KE  is  the  Kinetic 

Energy  round  and  HEAT  is  the  High  Explosive  Antitank 


are  thought  to  cause  the  greatest  errors  in  azimuth  and 
elevation.  Note  that,  originally,  the  figure  shows  optical 
ray  bending  causing  only  minor  errors  in  elevation. 

It  was  the  contention  of  the  ASL  report  that  the  full 
effects  of  atmospheric  refraction  were  not  being  considered 
in  current  tank  gunnery  doctrine,  particularly  those  effects 
created  in  a  hot/arid  or  semi-arid  environment.  The  ASL 
proposed  that  refraction  was  affecting  tank  gunnery  in  two 
ways.  First,  atmospheric  refraction  could  be  inducing  major 
elevation  errors  at  ranges  beyond  1500  meters.  Second, 
those  elevation  errors  currently  blamed  on  other  sources 
could  actually  be  due  to  refraction. 

1 .  Atmospheric  Conditions  Creating  the  Gunnery  Errors 

As  was  previously  described  in  Section  II,  atmospheric 
refraction  patterns  tend  to  follow  a  distinct  diurnal  pattern 
over  land.  During  the  day,  the  surface  is  heated  up  by  the 
sun,  and,  in  hot/arid  or  semi-arid  regions  this  heating  can 
cause  an  atmospheric  temperature  structure,  such  that  the 
temperature  decreases  rapidly  with  elevation.  This  condition 
can  cause  optical  turbulence,  subrefraction,  mirage  forma¬ 
tion,  or  a  combination  of  all  three.  At  night,  the  earth's 
surface  re-radiates  the  heat  absorbed  during  the  day  into 
the  atmosphere,  often  leading  to  the  creation  of  a 
temperature  inversion.  This  inversion  is  the  exact  opposite 
of  the  temperature  structure  that  characteristically  forms 
in  the  late  afternoon.  About  one-half  hour  before  sunset. 


and  some  time  ^ust  after  dawn,  a  condition  is  formed  when 
the  heat  radiating  from  and  into  the  ground  is  roughly 
equal. 

Local  weather  conditions,  such  as  the  presence  of  a 
high  pressure  region,  can  have  a  significant  impact  on 
atmospheric  refraction.  Depending  on  the  weather  conditions 
present,  the  diurnal  characteristics  described  above  can  be 
affected  in  such  a  manner  as  to  enhance  or  degrade  either 
the  subrefractive  condition,  the  superrefractive  condition, 
or  both.  Thus,  atmospheric  refractive  conditions  can  vary 
significantly  through  the  day,  and  from  one  day  to  the  next. 

The  ASL  report  was  primarily  concerned  with  the  effect 
of  going  from  a  subrefractive  condition  during  the  day  to  a 
superrefractive  condition  (characterized  by  the  temperature 
inversion)  at  night,  then  back  again  to  the  subrefractive 
condition. 

2 ,  Expected  Effects  on  Tank  Gunnery 

Since  optical  signals  bend  toward  the  denser  medium, 
the  direction  of  bending  for  the  diurnal  condition  described 
above  would  be  opposite  for  day  (subrefraction)  and  night 
(superrefraction).  This  is  graphically  shown  in  Figure  38. 

If  a  target  were  located  at  a  given  range,  its  image  would 
appear  to  shift  upward  at  night  and  downward  during  the  day. 


Figures  27A  and  27B  show  such  effects,  created  by  a  very 
strong  temperature  inversion  (superior  mirage)  and  an 
equally  strong  subrefractive  condition  (inferior  mirage). 


It  should  be  noted  that  a  mirage  need  not  be  present  for  the 
upward  and  downward  shifting  of  an  image  to  occur. 

Bending  of  optical  rays  can  cause  a  gunner  to  aim  at 
an  image  point  which  is  high  at  night  or  low  during  the  day, 
as  shown  in  Figure  44.  Since  most  tank  crews  boresight  and 
zero  their  tanks  during  the  day,  the  effects  of  subrefraction 
due  to  surface  heating  are  largely  removed.  But  this  can 
result  in  a  much  greater  error  at  night,  as  shown  in  Figure 
45.  It  should  be  reiterated  here  that  there  is  little 
difference  between  the  amount  of  refractive  bending  of 
visible  light  and  of  IR  energy.  Thus,  the  same  errors  occur 
going  from  an  optical  fire  control  system  using  visible 
light  during  the  day  to  one  that  uses  IR  at  night. 

3 .  Experimental  Results 

An  experiment  was  conducted  by  the  ASL  to  determine 
whether  the  theorized  refractive  effects  do,  in  fact,  occur. 
An  American  tank,  using  currently-deployed  optical  fire 
control  equipment  and  current  gunnery  doctrine,  was  used  in 
the  experiment,  which  took  place  at  White  Sands  Missile 
Range.  An  eight  foot  by  sixteen  foot  target  was  used.  The 
tank  was  boresighted,  synchronized,  and  zeroed  during  the 
day,  and  compensation  was  made  for  as  many  external  and 
internal  sources  of  error  as  possible.  Special  effort  was 
made  to  maintain  the  correct  boresight  while  firing.  If  the 
boresight  was  off  by  .1  mils  or  more  the  gun  was 
re-bore  sighted. 


Once  the  main  gun  had  been  satisfactorily  zeroed,  a 
series  of  5-round  shot  groups  was  fired  over  a  two-day  and 
three-night  period.  The  results  of  these  firings  are  shown 
in  Figure  46.  This  figure  clearly  shows  that  the  centers- 
of-impact  of  the  nightfired  groups  were  higher  on  the  target 
than  those  of  the  day  groups.  The  group  shifts  are  measured 
relative  to  the  center  of  impact  of  the  zero  confirmation 
group . 

The  inversion  and  temperature  lapse  strengths  can  be 
deduced  by  the  relative  distance  between  the  day  and  night 
5-round  groups.  As  can  be  seen,  the  greatest  distance 
between  shot  groups  occurred  between  the  first  night  and  the 
first  day,  and  the  differences  between  the  groups  pro¬ 
gressively  became  less  and  less.  On  the  first  night  the 
air  was  calm  and  clear,  resulting  in  the  creation  of  a 
strong  superre tractive  condition.  The  second  night  the  sky 
was  hazy,  weakening  the  superrefractive  condition.  On  the 
third  night  a  front  moved  in,  bringing  winds  and  clouds 
which  all  but  eliminated  the  superrefractive  condition.  The 
correlation  between  shot  placement  and  weather  conditions 
shows  that  refractive  errors  and  tank  fire  accuracy  are 
greatly  affected  by  weather  conditions. 

A  statistical  analysis  was  performed  by  the  ASL, 
proving  that  the  day  shots  and  night  shots  were  from  two 
statically  different  populations.  The  results  of  this 
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analysis  showed,  at  a  95%  confidence  level,  that  the  distance 
between  the  two  shot  centers  was  at  least  46  cm. 

To  bring  the  effect  of  the  difference  between  shot 
groups  into  perspective.  Figure  47  shows  how  many  of  the 
night  shots  would  have  missed  the  target  if  the  gunner  had 
actually  been  shooting  at  a  silouhette  of  a  T-62  tank  while 
aiming  at  the  base  of  the  turret.  This  figure  is  based  on 
the  results  of  a  vertical  T-test  of  the  5-round  shot  groups, 
which  determined  the  one  sigma  standard  deviation  in 
elevation  and  azimuth  for  each  group. 

Although  the  above  statistical  analysis  does  show 
significant  differences  between  day  and  night  shot  groups, 
it  does  not  prove  that  the  differences  are  due  to  refraction. 
To  prove  this  relationship,  temperature  measurements  at 
various  altitudes  above  the  ground  were  made  concurrently 
with  the  firings.  These  measurements  were  made  at  .5  meter 
invervals  from  O-to-4  meters  above  ground  level  (AGL). 
Apparent  elevation  angle  measurements  using  a  theodolite 
were  also  made  just  prior  to  each  firing.  There  measure¬ 
ments  were  made  so  that  the  line-of-sight  of  the  theodolite 
was  roughly  parallel  to  the  tank’s,  and  was  at  the  same 
height  AGL.  A  regression  analysis  was  then  conducted,  which 
showed  a  definite  correlation  between  the  temperature 
gradient,  the  theodolites  observed  absolute  elevation 
variations,  and  the  shot  pattern  elevation  shifts  during 


firing.  These  interactions  showed  a  significant  relationship 
between  refraction  and  firing  errors,  as  shown  in  Figures 
48  and  49. 

In  order  to  assess  the  importance  of  these  refractive 
effects,  based  on  the  probability  of  occurrence  of  the 
necessary  refractive  conditions,  the  ASL  investigated  the 
frequency  of  occurrence  of  those  atmospheric  conditions 
which  led  to  the  most  significant  firing  errors.  Figures  50 
and  51  give  a  general  idea  of  the  frequency  of  occurrence  of 
clear  skies  and  light  winds  in  southern  Germany  and  in  Iran, 
which  would  roughly  indicate  the  occurrence  of  refractive 
conditions.  Note  that,  for  both  cases,  these  conditions  can 
be  quite  common  at  certain  times  of  the  year. 

4 .  Conclusions 

The  results  of  the  ASL  experiment  showed  that  a 
definite  and  systematic  day/night  shift  in  elevation  angle 
can  occur  for  certain  meteorological  and  surface  conditions . 
The  values  measured  agreed  closely  with  what  was  predicted 
by  theory.  The  magnitude  of  the  measured  effect  was  as 
great  as  .8  mr.  At  longer  ranges  even  larger  effects  would 
be  expected.  Because  the  occurrence  of  temperature  gradients 
most  conducive  to  the  greatest  refractive  effects  is 
characteristic  of  hot/arid  and  semi-arid  regions,  it  is 
expected  that  gunnery  errors  of  this  type  will  occur  most 
frequently  in  the  desert  regions  of  the  world. 


Least  Squares  Slope 


igure  48,  Test  Results  vs  Predictions.  If  the  95%  line  slope  were  zero  or 
negative,  there  would  be  no  correlation  between  variables.  This 
graph,  however,  shows  that  there  is  a  definite  correlation 
between  the  theodolite’s  observed  absolute  elevation  variations 
and  the  shot  patterns  elevation  shifts  during  firing. 
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This  type  of  effect  is  a  characteristic  of  the 
atmosphere  and  not  of  any  particular  sight  design.  The  ASL 
is  developing  two  refractive  prediction  models  for  use  with 
tanks.  The  first  model  would  use  a  micro-computer,  require 
extensive  input,  and  provide  an  accurate  prediction  of 
existing  refractive  conditions.  The  second  model  would  be 
small  enough  to  be  placed  within  the  fire  control  computer 
of  the  tank,  and  would  only  use  input  readily  available  to  a 
tank  commander.  The  output  of  this  model  would  provide  the 
tank  commander  with  the  guidance  needed  to  correct  for 
predicted  refractive  effects.  Another  idea  would  be  to 
develop  simple  rules -of -thiomb  which  could  be  used  by  a  tank 
commander  to  predict  refractive  conditions  and  to  provide 
recommended  compensatory  actions. 

Theory  predicts  a  similar  effect  at  MMW.  If  true, 
these  effects  must  also  be  considered  in  any  research 
currently  under  way  into  using  this  region  of  the  EM 
spectrum. 

C.  COMPUTER  SIMULATION  1  -  THE  EFFECT  OF  REFRACTIVE 
CONDITIONS  ON  U.S.  TACTICAL  SYSTEMS,  USING  THE 
TROPOPLOT  COMPUTER  PROGRAM 

The  purpose  of  this  section  is  to  explore  atmospheric 
ref r activity  effects  on  actual  army  systems  operating  in  the 
RF  spectrum.  A  computer  program  called  "Tropoplot"  was 
used,  which  predicts  long  term  median  tropospheric 
transmission  loss  over  irregular  terrain.  This  computer 


program  was  originally  written  in  FORTRAN  by  the 
Environmental  Science  Services  Administration  (ESSA)  research 
laboratories  for  Control  Data  Computers  [Ref.  25].  The 
program  was  updated  for  use  on  IBM  computers  by  LT  James  M. 
Callaghan  [Ref.  26].  A  copy  of  the  program  is  included  in 
Appendix  B . 

1.  Tropoplot  Program  Description 

The  Tropoplot  program  uses  a  single  surface  refractivity 
value,  and,  although  developed  primarily  to  describe  the 
effects  of  atmospheric  refractivity  on  RF  signals,  calculates 
other  effects,  such  as  knife-edge  diffraction.  The  surface 
refractivity  value  is  related  to  a  refractive  gradient  through 
a  relationship  described  by  Bean  and  Dutton  [Ref.  27:  p.  63], 
However,  the  relationship  does  not  provide  gradient  values 
which  result  in  ducting  and  therefore  the  program  is  not  used 
to  evaluate  systems  operating  in  the  presence  of  ducts. 

Other  phenomena  which  affect  communications  systems,  such  as 
short  and  long  term  fading  adjustments  to  median  attenuations, 
are  accounted  for  in  the  program  by  using  factors  taken  from 
actual  measurements  made  on  various  communications  systems. 

All  variables  except  the  surface  refractivity  were  held 
constant  for  the  data  runs  in  order  to  investigate  only  the 
effect  of  refractivity  over  the  transmission  path. 

An  interdecile  range  is  used  to  characterize  the 
terrain.  The  ESSA  report  and  the  Callaghan  thesis  describe 
the  interdecile  height  correspondence  of  various  types  of 
terrain.  This  relationship  is  shown  below. 


TYPE  OF  TERRAIN 

Very  smooth  plains  or  water 

Smooth  plains 

Slightly  rolling  plains 

Rolling  plains 

IHlls 

Mountains 

Rugged  mountains 

Extremely  rugged  mountains 


ASYMPTOTIC  VALUE 
(meters ) 

CF  THE 

INTERCEDILE  HEIGHT 
0-5 
5-20 
20  -  40 
40  -  80 
80  -  150 
150  -  300 
300  -  700 
700  + 


The  program  is  usable  for  frequencies  from  20  MHz  to 
40  GHz  over  distances  of  up  to  2,000  kilometers.  The  surface 
refractivity  can  be  varied  from  250  to  400  N  units.  The 
ESSA  research  laboratories  [Ref.  25:  p.  1]  have  made  empirical 
measurements  for  comparison  with  the  program  results  and  have 
found  them  to  be  satisfactory. 

2.  AN/FPS-16  Radar 

The  first  system  used  in  the  program  was  the  AN/FPS-16 
air  tracking  radar.  This  high  power  system  can  be  found  at 
permanently  fixed  sites.  The  interdecile  range  for  "smooth 
plains"  was  used.  The  transmission  and  reception  line  losses 
were  not  available,  but  were  estimated  to  be  2  dB.  Many  of 
the  system's  parameters  are  provided  in  the  signal  strength 
data  input  section  of  Appendix  B's  data  run  output. 


A  "generic”  radar  warning  receiver  (RWR)  with 
characteristics  typical  of  fielded  systems  was  assumed  to  be 
located  on  an  aircraft  at  a  radar  range  of  approximately  125 
km  and  an  altitude  of  610  meters.  This  example  was  used  due 
to  its  similiarity  to  the  "standoff"  type  mission  profile 
that  might  be  flown  by  current  army  surveillance  and  signals 
intercept  airborne  platforms. 

During  the  data  runs,  the  surface  refractivity  was 
varied  from  250  N  units  to  400  N  units  in  steps  of  30  N  units. 
The  range  at  which  the  receiver  could  detect  a  signal  varied 
from  117  km  (at  250  N  units)  to  138  km  (at  400  N  units). 

This  represented  a  fairly  significant  increase  in  the 
distance  at  which  an  RWR  will  detect  the  radar  signal. 

The  next  step  involved  use  of  the  IREPS  program  to 
compare  the  same  radar  in  the  presence  of  a  strong  duct  and 
in  standard  atmosphere.  Although  IREPS  is  not  intended  for 
use  over  land,  the  comparison  is  helpful  in  exploring  the 
similarity  of  conclusions  for  Tropoplot  and  IREPS.  The 
outputs  to  the  Tropoplot  program  is  provided  in  Appendix  B. 

The  IREPS  results,  shown  in  Figures  52  and  53,  cannot 
be  meaningfully  compared  with  the  Tropoplot  results  for 
several  reasons,  including  the  fact  that  reflective 
coefficients  are  different  for  land  and  sea  water.  However, 
the  loss  diagram  for  IREPS  in  a  standard  atmosphere  versus 
that  for  a  surface  based  duct  is  interesting.  It  shows  a 


*»**  LOSS  DISPLRY 


FPS-16  TRKING  RDR  <I)ET> 


location:  not  specified 

DfiTE/TIME:  STANDARD 


0  40  B0  120  160  200 


RRNGE  IN  NRUTICRL  MILES 


FREE  SPACE  DETECTION  RANGE  BASED  ON  5054  POD  OF  A  1  3Q  M  AIR  TARGET. 


DASHED  LINE  INDICATES  DETECTION,  COMMUNICATION,  OR  INTERCEPT  THRESHOLD 


FREE  SPACE  RANGE:  163.8  NAUTICAL  MILES 
FREQUENCY*.  5825  MHZ 
TRANSMITTER/RADAR  HEIGHT:  65.0  FEET 
RECEIVER/TARGET  HEIGHT:  2000.0  FEET 


Figure  52.  IREPS  loss  display  for  the  FPS-16  tracking 
radar  in  a  standard  atmosphere. 
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LOSS  DISPLAY  **** 


FPS-16  TRKING  RDR  <DET) 


location:  33N  OOSE 
DATE^TIME:  SB  DUCT  540  FT 


0  100  200  300  400  500 


RANGE  IN  NflUTICRL  MILES 


FREE  SPACE  DETECTION  RANGE  BASED  ON  50^.  POD  OF  A  1  SQ  N  AIR  TARGET. 


DASHED  LINE  INDICATES  DETECTION,  COHMUNICATION,  OR  INTERCEPT  THRESHOLD 


FREE  SPACE  RANGE:  163.0  NAUTICAL  NILES 
FREQUENCY:  5825  MHZ 
TRANSniTTER/RADAR  HEIGHT:  65.0  FEET 
receiver/^targkT  height:  2000.0  feet 


Figure  53.  IREPS  loss  display  for  the  FPS-16  tracking 
.itTV  radar  operating  in  a  surface  based  duct. 
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loss  which  causes  the  radar’s  minimum  detectable  signal 
threshold  to  be  reached  at  about  7  4-  nm  in  a  standard 
atmosphere.  This  is  actually  the  horizon  limited  value  for 
the  given  heights  of  the  transmitter  (65  ft.)  and  target 
(2000  ft.).  In  a  strong  surface  based  duct,  the  path  loss 
is  much  less  and  the  threshold  is  reached  at  about  350  nm. 
The  loss  diagrams  are  included  in  Appendix  B. 

3.  Atf/VRC-12  Radio 

The  second  system  examined  was  a  tactical  radio  with 
parameters  similar  to  the  AN/VRC-12  series  used  extensively 
by  the  U.S.  Army  for  tactical  VHP  communications.  The 
parameters  were  obtained  from  the  U.S.  Army  Intelligence 
School,  at  Ft.  Devens,  Massachusetts,  and  from  FM  24-24 
[Ref.  28:  pp.  1-12]  and  FM  24-25  [Ref.  29:  pp.  2-16]. 
Parameters  that  were  not  available  were  estimated,  and  are 
considered  to  be  close  to  those  which  actually  exist. 

For  this  scenario,  two  radios  are  assumed  to  be 
operating  over  a  fairly  smooth  surface  with  their  antennas 
placed  on  terrain  which  allows  for  an  antenna  height  of  ten 
meters  above  the  immediate  land  in  order  to  enhance 
communications.  The  refractivity  was  again  varied  from  250 
to  400  N  units,  although  the  increments  were  increased  to 
50  N  units. 

In  this  example,  the  results  were  less  impressive 
than  in  the  previous  case.  The  maximum  range  for  reception 
only  varied  from  24  km  to  26.1  km.  Although  not  shown,  a 


simulation  of  a  ground  surveillance  radar  was  run,  using 
actual  parameters,  which  also  resulted  in  little  or  no 
chcinge  in  the  maximum  detection  range  of  3000  meters, 
regardless  of  surface  refractivity . 

The  results  of  these  latter  computer  runs  suggest  that 
refractivity  may  not  be  critically  important  for  shorter 
range,  lower  power  systems,  a  fact  that  is  not  totally 
unexpected, 

4.  Conclusions 

Although  it  does  appear  that  refractivity  plays  a 
critical  role  in  the  quantitative  effects  on  RF  signal 
propagation,  further  study,  using  a  program  for  irregular 
terrain  with  strong  vertical  refractivity  gradients ,  as 
opposed  to  a  surface  refractivity  value,  is  necessary  to 
define  more  precisely  the  critical  ranges  for  refractivity. 

At  present,  there  do  not  appear  to  be  any  programs  which 
have  this  capability. 

D.  COMPUTER  SIMULATION  2  -  THE  EFFECT  OF  HORIZONTAL 
REFRACTIVE  GRADIENTS  ON  AIRBORNE  DIRECTION  FINDING 
ACCURACY 

Airborne  Direction  Finding  (DF)  systems  provide 
critically  important  support  to  U.S.  Army  tactical  commanders. 
Systems  which  will  provide  enhanced  accuracy  are  currently 
being  planned  or  fielded.  In  general,  accuracy  is  dependent 
on  the  system's  ability  to  measure  individual  lines  of 
bearing  (LOB)  to  the  target  emitters  and  on  the  precision 
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with  which  the  navigation  system  can  determine  the  airborne 
platform’s  location.  Five  major  causes  of  error,  which  are 
independent  of  a  DF  system  but  which  affect  its  ability  to 
measure  LOB’s  are  reflections,  re-radiations,  diffraction, 
refraction,  and  fading  signals  [Ref.  10:  p.  2].  The  purpose 
of  this  section  is  to  investigate  the  impact  of  refraction 
on  DF  accuracy  through  the  use  of  computer  simulation. 

It  is  important  to  recognize  typical  situations  where 
refraction  could  cause  problems.  These  include  geographic 
areas  where  rapid  changes  in  temperature  and  humidity  (and 
sometimes  pressure)  occur  in  the  horizontal  plane.  Such 
areas  include  broad  river  valleys,  or  large  bodies  of  water 
surrounded,  at  least  partially,  by  land,  and  coastal 
regions.  In  such  areas  the  geometry  involved  in  taking  LOBs 
may  vary  greatly,  depending  on  where  the  refractive  gradient 
occurs  in  relationship  to  the  emitters  and  the  DF  platform. 
This  simulation  used  a  fixed  situation  from  which  it  is 
possible  to  extrapolate  to  other  circumstances.  This 
procedure  is  explained  further  in  sub-Section  4. 

This  section  is  organized  into  several  sub-sections. 
Section  2  provides  a  brief  general  description  of  DF  and  of 
the  tactical  scenario  used  in  the  simulation.  Section  3  is 
a  basic  description  of  the  program,  including  the 
assumptions  and  algorithm.  A  copy  of  the  program  is 
included  in  Appendix  C.  The  final  section  describes  the 
simulation  runs,  the  results  of  the  simulations,  and  the 
conclusions . 
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It  was  necessary  to  reduce  the  amount  of  output  from  the 
simulation  since  a  large  number  of  LOBs,  fixes,  and  other 
data  was  generated  on  each  "flight"  across  the  tactical 
area.  The  thesis  contains  examples  of  the  pertinent  data 
and  graphic  representation  of  the  output  in  this  section  and 
in  Appendix  C.  In  addition,  it  would  be  possible,  but  not 
practical,  to  do  a  number  of  sensitivity  analyses  by  changing 
the  variables  in  different  combinations.  For  example,  the 
depth  of  the  refractive  area,  the  horizontal  gradient,  or 
the  depth  of  the  emitters  could  be  changed,  either  singly  or 
in  combinations.  However,  the  purpose  here  was  to  make 
representative  changes  in  the  variables  from  which  it  would 
be  possible  to  draw  definitive  conclusions. 

1. 

Airborne  DF  systems  operate,  basically,  like  ground 
based  systems,  except  that  they  must  also  determine  the 
aircraft’s  location  at  all  times.  This  is  normally  done 
using  an  inertial  navigation  system  (INS)  with  some  type  of 
external  update.  The  aircraft  flies  a  preplanned  course, 
often  a  "racetrack"  pattern,  while  taking  LOBs  on  target 
emitters.  These  are  entered  into  emitter  LOB  files  and  are 
eventually  correlated  into  fixes,  a  term  used  to  describe 
the  target’s  calculated  location.  The  emitter  files  may  be 
sorted  by  frequency  or  by  other  parameters.  Most  systems 
use  further  correlation  or  refinement  to  "throw  out"  some  of 


General  Description  of  Direction  Finding  and  the 
Tactical  Scenario 
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the  fixes  which  do  not  fit  the  main  distribution  of  fixes. 
Reports  on  emitters  usually  include  an  elliptical  error 
probable  (EEP)  to  allow  for  the  expression  of  a  confidence 
factor  in  the  fix  accuracy.  The  simulation  computes  an 
"average”  fix  based  on  all  of  the  correlated  fixes. 

A  simple  type  of  DF  system  uses  a  "two-bearing”  cross 
(TBC)  to  produce  fixes.  In  this  System  the  point  where  two 
LOBs  cross  is  a  fix.  The  accuracy  of  the  fix  improves  as 
the  angle  formed  by  the  two  LOB’s  approaches  ninety  degrees. 
In  this  simulation,  a  TBC  routine  is  used  to  produce  fixes. 

The  simulation  used  a  corps  level  DF  mission  along  a 
corps  front  of  200  kilometers.  This  is  reasonable  for  a 
three  division  front.  A  single  aircraft  was  used,  although 
some  systems  use  two  or  more  in  order  to  provide  geometry 
that  will  allow  almost  instantaneous  fix  determination 
without  having  to  fly  long  distances.  The  aircraft  flew  at 
a  cruise  speed  of  200  knots  arid  required  a  little  over  3  0 
minutes  to  fly  the  front.  To  simplify  the  simulation,  the 
aircraft  flew  the  same  baseline  on  each  sweep  of  the  corps 
area. 

Four  emitters  were  used  in  the  scenario  and  their 
distance  from  the  corps  FLOT  could  be  varied  (see  Figure  54). 
The  deeper  locations  are  particularly  applicable  to  the 
tactical  concerns  of  a  corps.  Since  emitters  located  close 
to  the  corps  boundary  would  be  difficult  to  accurately 
fix,  due  to  poor  geometry,  the  emitters  in  the  simulation 
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were  placed  well  within  the  corps  lateral  boundaries.  Of 
course,  in  an  actual  situation,  some  compensation  could  be 
made  for  this  geometry  problem,  such  as  an  interface  with  DF 
assets  of  any  adjacent  corps.  Three  of  the  emitters  (A,  B, 
and  C)  are  set  up  to  evaluate  errors  from  refraction  and 
Emitter  D  is  used  to  compare  other  types  of  error,  such  as 
DF  system  error  or  platform  location  reporting  error,  with 
refractive  error.  Emitter  D,  therefore,  is  not  affected  by 
refractivity .  All  of  the  emitters  can  move  about  within  the 
200  kilometer  area  which  corresponds  to  the  aircraft 
baseline.  The  aircraft  was  assumed  to  be  flying  parallel  to 
an  area  with  a  refractive  gradient  similar  to  that  found 
along  a  coastline,  taking  fixes  on  land  based  emitters.  The 
gradient  area  has  a  uniform  horizontal  refractive  gradi-^nt 
(see  Figure  55).  The  depth  of  this  area  can  be  varies  n 
order  to  evaluate  the  effects  on  lines  of  bearing  which  are 
bent  over  long  or  short  distances.  Any  refractive  gradient 
value  can  be  used,  except  zero. 

With  the  exception  of  Emitter  D,  an  LOB  departs  the 
emitter  in  a  straight  line,  passes  through  an  area  where  it 
is  refracted,  and  departs  that  area  in  a  straight  line 
which,  like  the  previous  straight  line,  is  tangent  to  the 
arc.  The  LOB  arrives  at  the  baseline,  where  the  angle  of 
arrival  and  the  position  are  recorded  (see  Figure  55). 
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Emitter 


Figure  55.  Geometry  of  lines  of  bearing  from  an  emitter 


The  program  does  not  examine  the  effects  of  vertical 
changes  in  refractivity  since  the  ray  tracing  problem 
becomes  extremely  involved  when  working  with  both  horizontal 
and  vertical  ray  tracing. 

2 .  Algorithm  Description 

Throughout  the  program,  a  concerted  effort  was  made 
to  eliminate  as  many  assumptions  as  possible,  or,  at  the 
very  least,  to  make  assumptions  which  closely  resembled  real 
world  circumstances.  In  addition,  a  graphic  output  was 
chosen  over  other  output  forms,  such  as  statistical  tables, 
as  graphs  are  most  effective  in  quickly  evaluating  where  the 
fixes  were  computed. 

The  following  assumptions  are  pertinent  to  the 
understanding  of  the  simulation: 

1.  Flat  earth  geometry. 

2.  LOBs  are  taken  every  15  seconds  and  emitters  are 
selected  using  a  uniform  distribution. 

3.  LOB  angles  are  taken  with  respect  to  the  aircraft's 
longitudinal  axis . 

4.  All  emitters  can  be  received  at  all  times  during  the 
flight . 

5.  The  aircraft  track  is  along  the  Y  =  0  axis.  The 
emitter  depths  assume  a  standoff  range  in  order  to  provide 
aircraft  survivability  in  a  high-threat  air  defense 
environment . 

6.  Emitter  positions  are  fixed  for  each  mission. 
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In  reference  to  Emitter  D,  a  set  of  preset  flags 
controls  the  ”on-off"  condition  of  the  variables.  These  are 
described  in  the  first  part  of  the  program  (see  Appendix  C), 
Basically,  these  flags  control  introduction  of  LOB  error, 
navigation  error  (drift  or  random),  and  navigation  system 
update.  During  the  majority  of  the  simulation  runs,  LOB 
system  error  was  used  only  with  Emitter  D.  This  allowed 
comparisons  with  the  other  Emitters  (A,  B,  and  C)  which  had 
refracted  LOBs,  but  which  were  not  subject  to  system  error. 
Navigation  error  was  not  used  as  a  variable  in  any  of  the 
runs,  since  earlier  simulations  showed  it  to  be  less  of  a 
factor  than  LOB  system  error.  One  simulation  run  was  made 
in  which  Emitters  A,  B,  and  C  were  affected  by  both 
refraction  and  system  error,  as  noted  in  the  data  (see 
sub-Section  3). 

During  the  simulation,  the  flags  were  turned  on  only 
when  a  comparison  of  system  error  versus  refractivity  error 
was  desired. 


The  algorithm  which  follows  can  be  described  quite 
briefly,  even  though  the  program  itself  is  fairly  large.  An 
aircraft  is  started  on  a  track  with  constant  updating  of  its 
actual  and  reported  positions.  It  selects  an  emitter, 
calculates  the  LOB  to  the  emitter,  and  stores  both  the  LOB 
and  the  position  where  the  LOB  was  received.  In  the  case  of 
Emitter  D,  the  actual  and  reported  aircraft  positions  are 
recorded.  After  each  LOB  is  taken,  position  on  the  track  is 
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checked.  If  the  aircraft  is  at  a  boundary  of  the  corps 
area,  a  subroutine  is  called  in  to  correlate  the  LOBs  into 
fixes,  calculate  an  ’'average'*  fix,  and  plot  the  output  or 
summarize  the  output  statistically. 

In  order  to  validate  the  results,  simulation  runs 
were  made  using  identical  programs,  except  that  variable 
values  were  reported  after  each  step.  A  check  was  then  made 
of  selected  data,  both  qualitative  and  quantitative.  For 
example,  the  LOBs  were  checked  qualitatively  to  ascertain 
whether  their  values  approached  that  of  a  normal  drawn  from 
the  baseline  to  the  emitter,  as  the  aircraft  approached  the 
normal.  They  were  checked  quantitatively  by  using  a  hand¬ 
held  calculator  to  verify  the  angles.  Many  other  checks 
were  made,  including  a  run  with  the  dN/dX  (refractive 
gradient)  value  close  to  zero  in  order  to  determine  ii  any 
fix  error  was  present.  As  a  final  check,  large  variations  in 
gradient  area  and  refractivity  were  introduced  to  insure  they 
had  the  expected  effects  on  the  fix  accuracy. 

Figure  55  illustrates  the  basic  geometry  used  in  the 
simulation.  The  LOB  on  the  left  would  be  placed  into  one 
file  and  that  on  the  right  into  a  second  file.  Later,  after 
the  aircraft  had  reached  the  end  of  the  baseline,  the  two 
files  would  be  checked  and  the  two  LOBs  would  be  correlated 
into  a  fix.  The  correlation  routine  checks  every  possible 
combination  of  LOBs  which  produce  fixes. 


The  program  can  calculate  a  dN/dX  value  from  inputs 
of  temperature,  pressure,  and  hvunidity  for  two  different 
areas  separated  by  a  "gradient"  area,  as  used  in  the 
simulation.  In  addition,  a  dN/dX  value  can  be  inserted  into 
the  program  to  facilitate  sensitivity  analysis.  Different 
types  of  data  can  be  produced  by  use  of  "write"  statements 
located  throughout  the  program.  For  example,  the  program 
will  print  out  all  of  the  LOB  angles  along  with  the  locations 
where  they  were  taken. 

An  outline  of  the  basic  algorithm  follows : 

Initialize  variables. 

Calculate  a  dN/dX  value. 

Use  a  random  number  routine  to  vary  aircraft  start  on 

track. 

Establish  a  loop  to  control  round  trips  of  the  corps  front. 
Update  actual  aircraft  position. 

Check  to  see  if  aircraft  has  completed  left  to  right 
sweep;  otherwise  continue. 

Option  for  random  nav.  error,  update  reported  position. 
Option  for  drift  nav.  error,  update  reported  position. 
Option  for  LOB  error. 

Select  emitter  from  uniform  distribution. 

Check  aircraft  position,  if  past  point  where  a  normal 
drawn  from  the  emitter  to  the  baseline  would  hit  the 
baseline,  set  storage  flag  (Emitter  D  only). 

For  Emitters  A,  B,  and  C  use  a  slightly  different  routine 
to  calculate  the  angles  -  one  which  depends  on  the  aircraft 
location  on  the  baseline. 

Calculate  angle  to  selected  emitter,  in  radians. 


Convert  angle  to  degrees. 

Calculate  location  where  LOB  was  taken  (Emitters  A, 

B ,  and  C ) . 

Check  to  see  if  the  aircraft  position  is  past  the  normal 
points  (A,  B,  and  C). 

Call  a  subroutine  to  store  the  LOB,  actual  and  reported 
aircraft  position  (Emitter  D) . 

Call  a  subroutine  to  store  the  LOB  and  the  position 
where  the  LOB  was  received  (A,  B,  and  C). 

Check  to  see  if  a  right  to  left  sweep  is  complete. 

If  complete,  call  a  correlation  subroutine. 

Otherwise,  continue  flying  and  taking  LOBs. 

Update  aircraft  position. 

Go  to  emitter  selection. 

End  of  loop. 

Subroutine  to  store  LOBs  (4). 

Select  one  of  two  LOB  files  based  on  flag. 

Store  LOB,  actual  position,  and  reported  position  in 
proper  file  (D). 

Store  LOB  and  position  received  in  proper  file  (A,  B,  or  C). 

Return  to  main  program. 

Subroutine  to  correlate  LOBs. 

Establish  2  loops  to  check  every  possible  combination  on 
LOBs  in  the  storage  files  for  each  emitter. 

Establish  a  counter. 

Add  2  LOBs  together  and  make  a  check  for  acceptable 
geometry.  This  is  optional:  for  further  information, 
see  text. 

Determine  vertex  angle  and  coi.vert  to  radians . 
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Calculate  distance  between  both  LOBs  at  baseline. 

Convert  LOB  angles  to  radians. 

Calculate  length  of  one  side  of  two-bearing  cross  from 
vertex  to  baseline. 

Calculate  the  Y  coordinate. 

Add  in  Y  nav.  error  (D). 

Calculate  the  X  coordinate. 

Add  the  X  coordinate  to  reported  aircraft  X  position. 

Keep  a  running  total  of  the  X  and  Y  coordinates  for  all 
fixes . 

End  of  loops. 

Call  plot  routine  for  fixes. 

Calculate  average  X,  Y,  and  average  fix  error. 

Call  plot  routine  for  X,  Y,  and  average  error. 

Return  to  Main  Program. 

3.  Results  and  Conclusions 

The  graphs  and  tables  in  this  section  and  in 
Appendix  C,  represent  the  data  resulting  from  the  multiple 
runs  of  this  program.  Since  one  graph  for  each  emitter 
could  be  produced  on  each  sweep  of  the  corps  area,  it  was 
necessary  to  limit  the  amount  of  output  considerably.  This 
was  done  through  condition  statements  which  allowed  for 
saunpling  of  data  via  collection  from  selected  sweeps.  Due 
to  the  uniform  distributions  used  and  the  geometry,  there 
were  differences  from  sweep  to  sweep,  as  would  be  expected 
in  the  real  world.  Therefore,  data  used  for  comparison 
purposes  were  taken  from  the  same  sweep. 


One  of  the  first  problems  encountered  concerned  "angle 
gating".  This  problem  is  related  to  the  geometry  of  the 
two-bearing  cross  used  to  determine  fixes.  Ideally,  an 
angle  of  ninety  degrees  formed  by  two  LOBs  (known  as  the 
vertex  angle)  would  give  the  best  fix.  However,  because  of 
various  factors,  including  width  of  baseline,  emitter  depth, 
emitter  proximity  to  corps  boundaries ,  and  the  necessity  to 
have  several  fixes  for  better  averaging,  it  was  necessary  to 
specify  an  angle  gate  which  would  provide  acceptable 
geometry  and  yet  allow  for  several  fixes  on  each  sweep  for 
each  emitter.  This  was  done  by  using  a  hand-held  calculator 
to  determine  the  maximum  and  minimum  possible  vertex  angles 
for  each  emitter  at  all  ranges.  Throughout  the  data  runs, 
the  accuracy  was  highly  sensitive  to  the  angle  gates.  This 
results  from  the  fact  that  narrowing  the  gates  eliminates 
some  of  the  outliers  in  the  data. 

The  first  simulation  run  checked  the  program  by  using 
a  dN/dX  value  close  to  zero  to  indicate  negligible  fix  error. 

The  results  showed  that  double  precision  was  required  in  the 
program  (due  to  division  errors  by  numbers  close  to  zero). 

The  results  are  shown  in  data  Table  C-1  in  Appendix  C. 

The  initial  simulation  runs  compared  the  effects  of 
different  amounts  of  refractivity  on  the  fixes.  The  data 
is  presented  in  Graphs  1,  2,  and  3 ,  on  the  following  pages, 
and  summarized  in  Tables  C-2,  C-3,  C-4,  and  C-5,  in  Appendix  C. 
The  size  of  the  gradient  area  was  held  at  a  constant  15  kilometers 
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during  these  runs.  The  output  shows  that  the  refractive 
index  did  have  an  effect  on  fix  accuracy,  but  there  were  no 
cases  where  the  error  was  as  high  as  fixes  produced  with  an 
LOB  system  error  of  2.0  degrees  (Emitter  D).  Even  with 
dN/dX  =  -150.0,  the  errors  for  Emitters  A,  B,  and  C  did  not 
exceed  500  meters.  Part  of  this  is  due  to  a  ’’cancelling" 
phenomenon,  which  will  be  explained  further. 

A  gradient  of  -150  N  units  could  not  occur  over  a  15 
kilometer  area  since  the  range  of  N  values  would  greatly 
exceed  the  norm.  However,  there  is  a  case  where  the  geometry 
could  cause  individual  LOBs  to  travel  for  a  fairly  long 
distance  through  a  strong  gradient.  This  situation  is 
illustrated  in  Figure  56.  The  individual  LOB  departs  the 
emitter  and  travels  in  a  direction  which  roughly  parallels 
the  refractive  gradient.  This  results  in  the  ray  being  bent 
over  a  longer  distance,  thus  introducing  greater  error. 
Although  extensive  changes  to  the  program  would  have  been 
required  to  reproduce  this  particular  geometry  exactly, 
simulation  runs  were  made  with  large  gradients  over  areas  of 
15  to  30  kilometers  in  order  to  investigate  the  effect  of 
strong  gradients  and  longer  distances  on  LOB  accuracy.  The 
results  are  shown  in  Graphs  4,  5,  and  6,  on  the  following 
pages,  and  are  summarized  in  Tables  C-6  and  C-7  (Appendix  C). 
As  shown,  this  situation  produced  the  greatest  error  from 
refractivity ,  particularly  for  the  deepest  Emitter  (C), 
where  the  error  exceeded  700  meters.  The  data  appeared  to 
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Fiigure  56.  Diagram  of  a  Retracted  Line  of  Bearing 
Using  Alternate  Geometry. 
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reflect  some  ’’cancelling”  effect.  This  resulted  when  LOBs 
refracted  in  one  direction  were  correlated  with  LOBs  refracted 
in  the  opposite  direction,  causing  some  lateral  error  to  be 
cancelled  out.  The  error,  therefore,  occurred  mainly  in 
range  (Y-coordinate)  and  is  shown  in  the  fix  data  by  the 
reported  fix  being  deeper  in  range  than  the  actual  location. 

In  order  to  investigate  this  further,  a  run  was  made  where 
refracted  LOBs  from  Emitter  A  were  correlated  with  straight 
LOBs  in  a  situation  like  that  depicted  in  Figure  54  (i.e., 
refraction  occurs  over  a  large  area) .  An  LOB  system  error  of 
1.5  degrees  was  also  introduced.  The  results  are  shown  in 
Graphs  7  and  8  on  the  following  pages,  and  Table  C-8  of 
Appendix  C.  Note  that  Emitter  A’s  fix  error  greatly 
increased  from  previous  runs,  exceeding  Emitter  D’s  error  in 
fifty  percent  of  the  cases. 

The  primary  conclusion  to  be  drawn  from  the  simulation 
is  that  horizontal  refractive  gradients  can  be  a  problem  under 
certain  conditions,  but  that,  generally,  system  errors  will 
induce  greater  LOB  inaccuracy  than  will  refractive  effects. 

The  worst  case  is  the  problem  shown  in  Figure  56,  where  the 
LOB  is  refracted  over  a  long  distance  and  the  refracted  LOBs 
are  correlated  with  straight  LOBs.  However,  if  the 
possibility  of  these  conditions  exists,  measures  can  be 
tcdcen  to  minimize  their  impact.  For  example,  the  correlation 
software  routine  could  ignore  fixes  produced  by  LOBs  which 
came  through  an  area  of  suspected  strong  refractivity ,  or  a 
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correction  factor  could  be  introduced  for  refracted  LOBs. 

Since  good  refractivity  data  would  be  very  difficult  to 
obtain  in  the  emitter  area,  the  software  routine  could  employ 
a  predictor  system  for  LOBs  to  check  for  any  that  consistently 
arrive  at  the  baseline  at  a  different  angle  than  would  be 
expected.  This  would  require  close  operator  interface  and  an 
assumption  that  the  baseline  of  the  aircraft  is  long  enough 
to  allow  sufficient  LOBs  to  be  taken  which  do  not  pass  through 
an  area  of  strong  refraction.  Although  correction  tables  are 
sometimes  used  to  correct  DF  errors  arising  from  known 
refractive  conditions,  they  are  not  a  practical  solution  for 
mobile  DF  systems  CRef.  10:  pp.  4-7]. 

This  simulation  pointed  out  several  areas  which  could 
be  studied  further.  These  include  a  computer  simulation  which 
exactly  reproduces  the  geometry  of  Figure  56  or  one  which 
looks  at  a  time-of -arrival  system  working  in  an  area  of  strong 
refractive  gradients.  The  latter  could  involve  some  extensive 
work,  since  vertical  gradients  could  not  be  considered 
independently.  Both  vertical  and  horizontal  bending  would 
affect  the  time  of  arrival. 

The  magnitude  of  the  effect  of  horizontal  refractive 
gradients  is  dependent  on  many  factors,  including  the 
emitter  depth,  the  strength  of  the  gradient,  and  the  distance 
the  LOB  travels  through  the  gradient.  In  many  cases,  the 
effects  of  inherent  system  errors  will  be  greater  than  any 
refraction  induced  errors.  Since  refractivity  errors  are 


often  random  and  difficult  to  predict,  attempts  to  correct 
or  compensate  for  them  may  be  hard  to  implement.  Although 
improved  software  routines  and  other  techniques  may  be 
effective  in  improving  accuracy,  the  best  procedure  may 
be  to  recognize  possible  areas  of  strong  refractivity  and 
to  use  the  fix  results  accordingly. 


E.  SCENARIO 

The  purpose  of  this  Section  is  to  show,  by  means  of  a 
general  ground  combat  scenario,  the  large  number  and  variety 
of  U.S,  Army  systems  that  could  be  affected  by  atmospheric 
refractivity  and  the  wide  range  of  the  resulting  effects. 

1.  Scenario  Background 

To  provide  an  extreme  but  not  unusual  atmospheric 
refractive  situation,  a  Mid-Eastern  coastal  scenario  will  be 
used.  A  mechanized  division,  with  the  necessary  corps 
supporting  assets,  has  deployed  to  this  region  during  a 
period  of  increased  tension,  but  before  any  actual  combat 
has  occurred. 

The  mechanized  division  has  three  brigades  on  line, 
a  tank  battalion  held  in  reserve,  and  a  strong  covering 
force  forward  of  the  front  line  of  troops  (PLOT).  Since 
this  description  is  not  meant  to  be  a  lesson  in  ground 
combat  operations,  little  further  effort  will  be  spent  on 
describing  the  disposition  of  or  organization  for  combat  of 
the  elements  of  the  division. 
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The  division  has  one  flank  anchored  in  a  coastal 
mountain  range,  adjacent  to  another  U.S.  division,  and  the 
other  flank  is  anchored  in  the  sea.  The  division  is  currently 
located  in  a  blocking  position  across  a  major  high  speed 
avenue  of  approach,  as  shown  in  Figure  57, 

The  division  is  deployed  in  an  arid,  desert  region, 
characterized  by  large  areas  of  sand  and  rock  of  generally 
good  trafficability ,  interspersed  with  impassable  semi- 
forested  mountains,  sand  dunes,  and  ravines.  Very  little 
vegetation  exists  in  the  region  except  on  the  upper  slopes 
of  the  coastal  mountain  range.  The  temperatures  at  this 
time  of  year  range  from  the  70 ’s  at  night  to  110  degrees 
during  the  day.  Days  and  nights  are  often  clear,  but, 
vision  can  quickly  be  obscured  by  violent  sandstorms. 

In  this  scenario,  two  days  and  one  night  will  be 
considered.  On  the  first  day  the  extreme  daytime  tempera¬ 
ture  leads  to  the  creation  of  a  strong  subrefractive 
condition  over  the  coastal  strip  which  is  strengthened  by  an 
inferior  mirage  from  the  coastal  range  to  the  sea.  The  high 
temperatures  and  onshore  breeze  also  lead  to  high 
turbulence,  throughout  the  day,  along  the  coastal  strip. 

That  night,  a  high  pressure  area  begins  moving 
slowly  into  the  region,  enhancing  the  superrefractive 
condition  created  by  the  radiation  of  heat  from  the  desert 
surface.  Some  turbulence  exists  but  not  nearly  so  severe  as 
that  which  existed  during  the  day. 
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During  the  second  day,  the  high  pressure  area  leads 
to  the  creation  of  an  offshore  flow  of  dry,  hot  air  and 
results  in  the  creation  of  a  strong  elevated  layer  and 
surface  duct  that  extends  inland  almost  to  the  coastal  range 
2 .  Atmospheric  Refractive  Effects 
a.  Optical  Effects 

The  presence  of  the  subrefractive  condition, 
inferior  mirage,  and  optical  turbulence  on  the  first  day 
will  significantly  affect  many  of  the  division’s  optical 
surveillance,  fire  control,  and  tracking  systems. 

The  distance  at  which  the  mirage  begins  will 
limit  visual  observation,  by  ground  observers,  of  areas  to 
the  front  of  the  division’s  covering  force.  Mirages  can 
significantly  affect  optical  systems  as  close  as  1000  meters 
In  any  combat,  this  near  boundary  of  the  mirage  can  reduce 
the  detection  and  engagement  ranges  and  accuracy  of 
direct-fire  weapons,  particularly  if  optical  turbulence  is 
also  severe.  In  a  like  manner,  the  presence  of  the  mirage 
condition  and  optical  turbulence  could  also  reduce  the 
effectiveness  and  range  of  laser  rangefinders  and  laser 
designators . 

Some  relief  from  these  adverse  refractive 
effects  can  be  gained  by  placing  observers  as  high  above  the 
surface  of  the  coastal  strip  as  possible.  Those  observers 
located  above  the  top  of  the  mirage,  either  on  the 
surrounding  high  ground  or  in  aircraft,  will  be  unaffected 
by  the  mirage. 


190 


Another  adverse  refractive  effect  would  result 
if  the  division’s  tank  main  guns  are  bores ighted  and  zeroed 
during  the  first  day  and  fired  the  first  night  or  second 
day.  As  explained  in  Section  III.B.,  the  result  of  such  a 
situation  could  be  significant  reduction  in  tank  main  gun 
accuracy  during  night  firing  and  during  the  second  day. 

Some  bending  of  IR  paths  will  occur  at  night  due 
to  the  superrefractive  condition,  but  the  effect  of  this  on 
FLIR  systems  would  probably  be  insignificant. 

The  ducting  on  the  second  day  would  result  in 
a  weakening  or  total  absence  of  the  inferior  mirage  seen  on 
the  first  day.  Thus,  visual  observation  and  laser  designation 
and  rangefinding  systems  would  be  more  effective  on  the 
second  day,  although  optical  turbulence  might  still  be  quite 
strong . 

b.  Effects  on  Communications  and  Electronic  Warfare 

Systems 

Desert  regions,  in  general,  have  an  adverse 
effect  on  RF  communications.  Some  of  the  effects  include 
difficulty  in  adequately  grounding  communications  systems, 
existence  of  a  weak  HF  ground  wave,  and  a  generally  low 
surface  reflectivity,  all  of  which  lead  to  a  significant 
shortening  of  HF,  VHF,  and  UHF  communication  ranges. 

The  presence  of  the  subrefractive  condition  on 
the  first  day,  although  causing  some  bending  of  ground-to- 
ground  communications  signals  and  ECM  jamming  signals,  will 
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not  result  in  any  additional  significant  loss  of  range.  Thi 
is  not  necessarily  the  case  for  air-to-ground  communications 
and  jamming.  The  existence  of  the  subrefract ive  layer  over 
the  surface  of  the  coastal  strip  will  set  a  limit  on  the 
minimum  angle,  from  the  aircraft  to  the  ground  station,  at 
which  communications  will  be  dependable  between  the  two 
stations.  This  latter  effect  will  also  place  some  limit  on 
friendly  airborne  communications  intercept  and  direction 
finding  (DF)  capabilities  for  platforms  located  behind  the 
FLOT  and  attempting  to  acquire  and  DF  hostile  communications 
emitters.  The  height  of  the  aircraft  above  ground,  the 
distance  to  the  ground  communications  system,  the  power 
output  of  the  airborne  system,  and  the  strength,  height,  and 
extent  of  the  subrefractive  condition  are  all  critical  in 
determining  the  overall  effect. 

The  accuracy  of  airborne  and  ground  based  DF 
systems  can  also  suffer  from  a  horizontal  atmospheric 
refractive  gradient  created  by  the  juxtaposition  of  land  and 
sea  air  masses  along  the  coastal  zone.  These  errors  can  be 
quite  significant,  as  shown  in  Section  III.D. 

The  night  time  superrefractive  condition  and  the 
ducting  occurring  during  the  second  day,  both  of  which  are 
surface  based,  will  tend  to  extend  communication  ranges  of 
those  signals  operating  within  these  conditions  and  above  or 
near  the  minimum  trapping  frequency,  f^.  This  will  serve  to 
increase  the  distance  between  ground  based  units,  for  which 
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conununi  cat  ions  will  be  considered  reliable,  and  also  the 
distance  from  jammer-to-target  at  which  jamming  will  be 
effective.  Communication  intercepts  of  distant  hostile 
emitters  will  similiarly  be  enhanced.  In  a  like  manner, 
hostile  Electronic  Warfare  (EW)  units  will  be  able  to  jam, 
intercept,  and  DF  friendly  communications  from  greater 
distances.  Terrain  features  may,  however,  override  the 
effect  of  the  superrefractive  condition,  and,  in  any  event 
the  actual  enhancement  may  be  only  a  few  kilometers  for 
ground  based  systems,  as  shown  in  Section  III.C. 

c.  Effects  on  Noncommunications  Emitters 

The  division’s  noncommunications  emitters, 
primarily  radar,  are  spread  throughout  the  division  and 
serve  a  variety  of  purposes. 

The  accuracy  of  the  division’s  counterbattery/ 
countermortau?  radars  may  be  degraded  by  the  subrefractive , 
superrefractive,  and  ducting  conditions.  The  extent  and 
nature  of  such  atmospheric  refraction  produced  errors  in 
the  Firefinder  radar  is  currently  a  subject  of  interest. 

Ground  surveillance  radar  systems  will  suffer 
very  little  from  atmospheric  refractive  effects,  since  these 
systems  are  operated  over  relatively  short  ranges  and  are 
directed  along  the  surface  of  the  earth. 

Ground  based  air  defense  radars,  however,  can  be 
adversely  affected  by  all  forms  of  atmospheric  refraction. 
The  first  day’s  subrefractive  condition  will  result  in 
errors  in  both  angle  and  range,  but  will  be  unlikely  to 


affect  early  warning  or  tracking  radars  significantly.  The 
first  night's  superrefractive  condition  will  cause  radar 
signals  to  bend  in  the  opposite  direction  from  those 
propagating  during  the  day,  and  will  thus  cause  angle  and 
range  errors  in  air  defense  radars;  but,  again,  because  of 
the  nature  of  the  radar's  purpose,  these  errors  will 
probably  be  insignificant.  The  strong  ducting  occurring  on 
the  second  day  will  result  not  only  in  larger  angle  and 
range  errors  than  those  occurring  at  night,  but  will  also, 
in  all  probability,  result  in  the  development  of  radar 
holes.  These  radar  holes  could  cause  reduced  detection 
ranges  for  early  warning  radars. 

Airborne  radars,  such  as  the  side-looking 
airborne  radar  (SLAR),  will  suffer  a  reduction  in  range  as  a 
result  of  the  subrefractive  conditions  occurring  on  the 
first  day,  if  the  aircraft  is  flying  above  the  condition. 

The  superrefractive  and  ducting  conditions  of  the  first 
night  and  second  day  may  result  in  an  extended  range  for  the 
SLAR  system,  but  this  would  depend  on  the  height  at  which 
the  aircraft  is  flying  and  the  height  of  these  refractive 
conditions. 

3 .  Conclusions 

It  can  be  seen  that  a  multitude  of  combat  systems  as 
well  as  surveillance  and  intelligence  gathering  systems  are 
affected  by  atmospheric  refractive  conditions.  An  ability 
to  predict  the  existence  of  such  conditions  can  go  a  long 


way  towards  explaining  irregular  or  ’’unusual"  effects  on 
army  ground  systems .  The  ability  to  forecast  atmospheric 
refractive  conditions  will  also  provide  an  opportunity  to 
take  advantage  of  those  effects  that  can  enhance  certain 
capabilities  and  to  compensate  for  those  that  will  degrade 
other  capabilities. 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

It  has  been  the  purpose  of  this  thesis  to  highlight  the 
effects  of  atmospheric  refraction  on  ground  force’s  systems 
that  utilize  EM  waves  propagating  through  the  atmosphere. 

To  this  end,  the  various  forms  of  atmospheric  refractivity 
in  the  troposphere  have  been  described,  and  their 
characteristic  effects  on  EM  waves  have  been  examined. 
Additionally,  specific  examples  and  simulations  were  used  to 
expand  and  to  reinforce  some  of  the  points  brought  out 
earlier  in  the  thesis. 

From  the  discussions,  examples,  and  simulations  presented 
herein,  a  number  of  points  should  be  apparent.  First, 
atmospheric  refraction  in  several  related  but  different 
forms  does  affect  the  propagation  of  EM  waves  traveling 
through  the  atmosphere.  Second,  the  specific  effects  of 
atmospheric  refraction  on  the  different  regions  of  the  EM 
spectrum  will  vary,  not  only  with  wavelength,  but  with  the 
form  of  the  refractive  condition.  Third,  to  some  degree  the 
forms  of  atmospheric  refraction  can  be  forecast,  given 
certain  meteorological  parameters.  In  a  like  manner, 
knowledge  of  the  earth  region  of  interest  can  be  used  to 
predict  what  form  or  forms  ^f  atmospheric  refraction  to 
expect.  Fourth,  specific  military  systems  using  the  EM 


spectrum  are  affected  by  refractive  conditions  in  the 
atmosphere,  and  these  effects  result  in  the  associated 
weapons  system’s  performance  also  being  affected.  Fifth, 
the  significance  of  atmospheric  refraction  on  any  given 
system  is  dependent  on  the  form  of  the  refraction,  the 
strength  of  the  effect,  the  form  and  purpose  of  the  affected 
system,  and  the  tactical  scenario.  There  is,  therefore,  no 
easy  answer  to  the  question  of  whether  atmospheric  refraction 
is  significant  to  a  system's  performance  or  whether  it  occurs 
often  enough  to  be  considered  a  significant  effect.  This 
latter  consideration  is  dependent  on  the  system  affected  and 
on  the  region  of  interest.  Sixth,  tactics  can  be  developed 
to  compensate  for  or  to  utilize  certain  atmospheric 
refractive  effects.  Finally,  it  should  be  apparent  that 
atmospheric  refractive  effects  on  ground  warfare  should  be 
given  consideration  by  the  tactical  commander. 

The  ability  to  forecast,  predict,  and  understand  the 
effects  of  atmospheric  refraction  allows  the  commander  to 
avoid  being  surprised  by  many  of  the  "unexpected”  and 
"unexplained"  problems  that  occur  in  systems  using  EM 
radiation.  This  ability  also  provides  the  commander  with 
opportunities  to  use  certain  refractive  effects  to  his  own 
advantage . 

The  tank  gunnery  results  and  the  simulations  in  this 
thesis  indicate  that  the  U.S.  Army's  primary  problems  with 
atmospheric  refraction  will  occur  in  the  electro-otpics 
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spectrum.  There  are  some  cases  in  the  RF  spectrum  where 
refractive  effects  are  fairly  significant,  as  when  lines  of 
bearing  travel  for  long  distances  through  a  refractive 
gradient.  However,  the  effect  on  electro-optics  (and 
possibly  millimeter  waves)  over  relatively  short  distances 
is  more  significant.  Due  to  the  shorter  operational  ranges 
of  most  EO  systems,  as  compared  with  RF  systems,  the 
possibility  of  detecting  and  correcting  refractive  problems 
in  this  region  may  also  be  greater. 

In  short,  it  is  hoped  this  thesis  has  shown  atmospheric 
refraction  to  be  worth  further  consideration  and  study.  In 
the  next  sub-section  a  list  of  recommendatiors  for  the 
direction  of  further  studies,  is  provided. 

B.  RECOMMENDATIONS 

Following  are  recommendations  for  further  study  of  the 
effects  of  atmospheric  refraction  on  ground  warfare,  and 
additional  recommendations  as  to  when  and  where  this 
subject  should  be  considered. 

1.  Further  study  of  the  wavelength  dependence  of 
refractive  effects  is  needed. 

2.  Atmospheric  refraction  should  be  a  consideration  in 
system  design,  particularly  in  system  testing. 

3.  Study  is  needed  concerning  atmospheric  refractive 
effects  on  specific  systems  (e.g.  AH64  FLIR,  Firefinder 
Radar,  JSTARS,  etc.). 


4.  As  a  result  of  studying  the  effects  of  atmospheric 
rerfraction  on  specific  systems,  methods  should  be  developed 
to  compensate  for  and/or  take  advantage  of  these  effects. 

5.  An  effort  should  be  made  to  relate  the  effects  of 
dispersion,  extinction,  reflection  (multipath),  and 
diffraction  into  the  overall  picture  of  atmospheric 
propagation. 

6.  Further  study  is  needed  on  determining  the  causes  of 
atmospheric  refraction  in  its  many  forms ,  and  into  those 
meteorological  parameters  needed  to  forecast  atmospheric 
refractive  conditions  accurately. 

7.  A  greater  understanding  of  the  reflection 
coefficients  for  various  vegetation  and  terrain  types  and 
their  impacts  on  EM  propagation  is  needed  before  a 
prediction  model  can  be  developed. 

8.  Development  of  a  systems  prediction  model  along  the 
lines  of  the  Navy's  IREPS  model  is  a  must. 

9.  Once  a  prediction  model  is  under  development,  a 
study  will  be  needed  to  determine  who,  within  the  Army's 
organization,  will  measure  the  meteorological  parameters 
needed  to  predict  atmospheric  refractive  effects. 
Determinations  must  also  be  made  as  to  whom  the  responsibility 
for  op€*rating  the  prediction  model  will  be  given,  to  whom 

the  data  will  be  disseminated,  and  to  whom  the  data  will  be 


disseminated . 


APPENDIX  A 


IREPS 

The  U.S,  Navy  has  been  aware  of  the  effects  of  the 
atmosphere,  particularly  of  atmospheric  refract ivity,  on  EM 
waves  since  at  least  World  War  II.  However,  until  recently 
the  exact  nature  of  these  atmospheric  refractive  effects, 
and  therefore  an  ability  to  accurately  predict  them,  has 
been  non-existent. 

With  the  introduction  of  an  Integrated  Refractive  Effects 
Prediction  System  (IREPS)  by  the  Naval  Ocean  Systems  Center 
(NOSC)  of  San  Diego,  the  Navy  has  now  developed  the 
capability  to  predict  atmospheric  refractive  conditions  at 
sea.  IREPS  uses  a  desk-top  programmable  computer  (HP  9845) 
with  a  graphic  display  to  provide  evaluations  of  the  effects 
of  refraction  on  surveillance,  communication,  electronic 
warfare,  and  weapon  systems.  This  system  is  currently 
being  installed  on  all  Navy  aircraft  carriers,  to  be  used  as 
an  aid  in  developing  tactics  to  eissess,  compensate  for,  and 
take  advantage  of  refractive  effects  at  sea.  Eventually, 
other  atmospheric  effects  will  also  be  integrated  into  the 
IREPS  program,  along  with  the  PROCAL-PREOS  airborne  FLIR 
performance  model,  in  order  to  develop  a  more  comprehensive 
prediction  system  for  future  use. 


It  IS  important  to  note  that  the  effects  of  concern  in 
IREPS  are  only  significant  at  EM  frequencies  above  100  MHz 
are  considered  only  for  the  lower  atmosphere  (troposphere) 
and  only  over  water. 

1.  Products  of  IREPS  Model 


There  are  eight  basic  products  that  can  be  generated 
by  IREPS,  according  to  the  IREPS  Revision  2.0  User's  Manual 
[Ref.  30:  p.  30].  Each  product  is  dependent  on  a  variety  of 
data  that  must  be  entered  into  the  program.  For  the  first 
product,  only  the  latitude  and  longitude  are  needed,  but  for 
the  remaining  seven  products,  a  variety  of  on  site  environ¬ 
mental  data  are  required.  Each  product  is  displayed  on  the 
minicomputer  screen,  and  can  be  printed  on  an  8  1/2-by-ll-inch 
page.  The  eight  products  are: 

a.  The  Historical  Propagation  Conditions  Summary 
Figure  58  is  an  example  of  the  Historical 

Propagation  Conditions  Summary.  This  product  provides  a 
description  of  the  refractive  average  conditions  for  the 
geographic  location  specified  by  an  input  longitude  and 
latitude.  A  refractive  profile  can  also  be  generated  for 
this  region  by  season  and  by  conditions  of  darkness. 

b.  The  Environmental  Data  List 

Figure  59  is  an  example  of  the  Environmental 
Data  List,  used  primarily  for  checking  values  of  radiosonde 
data  entered  into  the  program,  such  as  pressure,  temperature, 
and  relative  humidity.  Values  that  are  produced  and  printed 
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EVAPORATION  DUCT  HISTOGRAM  IN  PERCENT  OCCURRENCE 
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Figure  58.  Historical  Propagation  Conditions 
Summary  Product. 
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Figure  59.  Environmental  Data  List  Product 


out  in  the  Environmental  Data  List  are  the  dewpoint 
depression,  altitude,  N  unit,  N  unit  gradient,  M  unit  and 
a  description  of  the  refractive  condition. 

c.  The  Propagation  Conditions  Summary 

Figure  60  is  an  example  of  the  Propagation 

Conditions  Summary,  showing  the  existence  of  refractive 
conditions  for  a  specified  region  at  a  given  date  and  time. 
The  summary  provides  a  plot  versus  altitude  diagram  of  N 
units  and  M  units.  The  extent  of  any  duct  present  is  shown 
by  the  use  of  shaded  areas  in  the  block  on  the  right  hand 
side  of  the  summary.  A  brief  statement  on  the  performance 
of  several  types  of  EM  systems  operating  in  the  specified 
region  is  provided  at  the  bottom  of  this  summary. 

d.  Surface-Search  Radar  Range  Tables 

Figure  61  is  the  Surface-Search  Radar  Table 
format.  This  product  provides  detection  range  predictions 
for  naval  surface-search  radars.  An  actual  example  of  the 
output  is  not  provided  since  the  predicted  range  values  are 
classified  CONFIDENTIAL. 

e.  The  Coverage  Display 

Figures  19  and  20  are  examples  of  the  Coverage 
Display  product.  This  product  shows  the  area  of  coverage  of 
a  specific  type  of  radar  or  radio  emitter  on  a  curved  earth 
range-versus-height  plot.  The  shaded  area  is  the  area  of 
detection  or  communication  for  the  system  represented.  The 
impact  of  atmospheric  refraction  on  a  system  can  be  seen 
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Ij.  U  U  h- 


Fi- OFriuftTIOft  CONDITIONS  SiJIlllHP'r 


♦  •  *  « 


location:  31  I>SN  113  3eN 
DATE  time:  ir  JUN  00452 


9K-  \ 


210  240  270  300  330  360  330  430  530  630  730 


REFPflCTIVITY 
N  UNITS 


MODIFIED  REFRACTIVITY 
M  UNITS 


MIND  SPEED*  12,0  KNOTS 


EVAPORATION  DUCT  HEIGHT*  8.5  METRES 

»  28.0  FEET 


:upfai:e-to-surface 

e:;tended  ranges  at  all  frequencies 

SURFACE-TO-AIR 

EXTENDED  RANGES  FOP  ALTITUDES  UP  TO  1,022  FEET 
POSSIBLE  HOLES  FOP  HuTITUDES  ABOVE  I , 022  FEET 

AlP-lO-AIR 

EXTENDED  RANGES  FOR  ALTITUDES  UP  TO  1,072  FEET 
POSSIBLE  HOLES  FOR  ALTITUDES  ABOVE  1,072  FEET 

:uFFAi. t  pefractivity:  341  --set  5PS-48  to  344 


Figure  60.  Propagation  Conditions  Sununary  Product. 
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LOCftTION: 

time: 


*****  SURFfiCE  SEftRCH  RfiDftR  RANGE  TABLE  ***** 


SURFACE  SEARCH  RADAR: 
RADAR  ANTENNA  HEIGHT: 


U  S  SHIP  TYPE/CLASS 

TTTTVH - 

CG/CGN 

DD/DDG 


FEET 

}  RANGE  IN  H 

AVG  1  MAX 


LST 
AE-  AF 

AO 'AOE^'AOR 


SOVIET  SHIP  TYPE/CLASS 

"yiT-'V  CLASS - 

MOSKVA  CLASS 
CLG 

OD  DDG 
FRIGATE 

~'CPffVETTE - 

OSA--STENKA  CLASS 
PRIMORYE  CLASS  AGI 

■  LENTRfl  CLASS  flCl - 

OKEAN  CLASS  AGI 


“UETETT 

MIN 

nJFT  RHTTCI 
AVG 

1 . IK  HI 

MAX 

1 

1 

1 

1 

EVAPORATION  DUCT  HEIGHT* 


METRES. 

FEET 


Figure  61.  Surface  Search  Radar  Range  Table  Format. 


graphically  by  comparing  printouts  of  the  radar  coverage 
with  and  without  refractive  effects  present.  Note  also  that 
the  presence  of  a  radar  hole  can  be  seen  in  this  product  by 
the  shape  of  the  coverage  pattern  above  the  duct. 

f.  The  Loss  Display 

Figure  62  is  an  example  of  a  Loss  Display 
showing  one-way  path  loss  in  dB  versus  range.  The  dashed 
line  on  the  printout  is  the  threshold  for  detection, 
communications,  or  intercept.  The  solid  line  is  the 
system's  signal.  Thus,  a  determination  of  where  the  signal 
would  be  strong  enough  to  be  detected  can  be  estimated  from 
this  product. 

g.  The  Radiosonde  Observation  Analysis 

This  product  is  provided  largely  to  assist 
meteorological  personnel,  and  consists  of  a  radiosonde 
listing,  a  skew  T,  a  log  P  thermodynamic  diagram,  and  an 
encoded  upper  air  observation  message. 

h.  The  ESM  Intercept  Range  Table 

Figure  63  is  the  format  for  the  ESM  Intercept 
Range  Table.  This  table  provides  the  maximum  intercept 
range  for  several  different  Naval  ESM  systems  versus  a 
variety  of  hostile  emitters.  An  actual  example  is  not 
provided  since  much  of  the  information  on  the  hostile 
emitters  and  on  the  predicted  values  is  classified. 
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2.  Limitations  of  the  IREPS  Model 


There  are  a  number  of  limitations  in  the  IREPS 
model,  some  of  which  are  being  corrected.  These  limitations 


include : 


a.  Limited  frequency  range. 

b.  Neglect  of  the  effect  of  land  and  sea  clutter. 

c.  Assumption  of  horizontal  homogeneity  in  the 


refractive  structure . 


d.  Limits  on  the  maximum  and  minimum  ground 


based  antenna  height. 


e.  Exclusion  of  the  effects  of  sea-reflection 


interference  on  airborne  systems. 


polarized . 


f.  Assumption  that  all  EM  systems  are  horizontally 


g.  Absorption  is  not  accounted  for. 

h.  Limits  exists  on  the  consideration  of  low 


level  elevated  ducts. 

i.  The  coverage  diagram  is  not  applicable  to  all 
radar  applications. 

3.  ApDlicabilitv  of  IREPS  to  Over-land  Refractive 


Prediction 

While  IREPS  is  proving  to  be  of  considerable  value 
to  the  Navy  even  with  its  wide  range  of  limitation,  its 
applicability  over  water  does  not  imply  that  it  will  work 
equally  well  over  land.  IREPS  in  its  present  form,  because 
it  was  designed  for  over-water  use,  does  not  account  for 


terrain  features,  vegetative  effects,  soil  type,  or 
variation  in  the  surface  reflective  coefficient,  all  of  which 
have  great  effect  over  land.  In  addition,  the  nature  of 
atmospheric  refractivity  over  land  is  much  more  transient 
and  localized  than  over  water,  and  this  must  also  be  taken 
into  account  in  developing  an  over-land  model.  The  basic 
design  of  IREPS,  however,  may  prove  to  be  viable  in  the 
development  of  an  over-land  model  and,  at  the  very  least,  is 
worthy  of  further  study  in  this  regard. 


APPENDIX  B 


COMPUTER  SIMULATION  1  LISTING  AND  OUTPUT 

This  Appendix  contains  the  program  and  output  for 
Simulation  1,  Each  computer  run  is  summarized  in  four 
pages:  Summary  of  the  Parameters,  a  Transmission  Loss  Table 

A  Loss  Diagrcun,  and  a  Signal  Strength  Table.  The  latter 
contains  the  refractivity  value.  The  type  of  system  (either 
radar  or  communications)  can  be  identified  by  the  frequency 
at  the  top  of  the  first  output  page  for  each  run. 
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WRITE  •«OISTANCE  ANO  SIGNAL  STRENGTHm  HEADING  AND  DATA 
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*••••  FOR  *  RECEIVER  INFUT  LEVEL  OF  -40.33  06N.  DISTANCE  OETNEEN  RNTENNAS  CAN  BE  AT  LEAST  130.00 
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FOR  «  RECEIVER  INPUT  LEVEL  OF  -60.00  UON.  OISTANCE  BETUEcN  ANTENNAS  CAN  BE  AT  LEAST  26.10  KN  666*6 


APPENDIX  C 


COMPUTED  SIMULATION  2  LISTING  AND  OUTPUT 

This  Appendix  contains  the  program  and  output  for 
Simulation  2.  The  output  consists  of  Tables  which  are 
referenced  by  nvimber  in  the  main  text  of  the  thesis,  i 
the  top  of  each  data  Table,  is  an  explanation  of  the 
parcuneters  used  in  the  run.  Each  Table  represents  one 
run  of  the  program. 
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TABLE  C-1 

Simulation  2,  Run  1.  Simulation  to  check  program. 
Gradient  area  equal  to  15  Km.  All  values  are  in 
kilometers.  Emitter  D  has  system  error  of  2.0 
degrees. 


0NDX«  -  1.300 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y-IOO. 001012  AVERAGE  X»  80.000079 
AVERAGE  MISSEO  DISTANCE*  0.001015 


EMITTER  a  AVERAGE  ERROR 

AVERAGE  y*  70.000441  AVERAGE  X*119. 999938 
AVERAGE  MISSEO  DISTANCE*  0  .  000  445 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y*250. 002003  AVERAGE  X*149. 999896 
AVERAGE  MISSEO  DISTANCE*  0.0020  06 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y*  99.999820  AVERAGE  X*  79.999996 
AVERAGE  MISSEO  DISTANCE*  0.000180 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y*100.000d53  AVERAGE  X*  80.00  0006 
AVERAGE  MISSEO  DISTANCE*  0.000653 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y*  70.000413  AVERAGE  X*119. 999964 
AVERAGE  MISSEO  DISTANCE*  0.000415 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  7*2  5  0.0022  6  3  AVERAGE  X*149 .999909 
AVERAGE  MISSED  DISTANCE*  0.  C022  64 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y*  99.999889  AVERAGE  X*  79.999996 
AVERAGE  MISSED  DISTANCE*  0.000111 
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TABLE  C-1 


Simulation  2,  Run  1  (continued). 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y«10a. 000803  AVERAGE  X«  80 

AVERAGE  HISSED  DISTANCE*  0.  COO 8 04 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y«  70.000464  AVERAGE  X»119 

AVERAGE  MISSEO  DISTANCE*  0.  000  4  71 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y*250. 002500  AVERAGE  X*149. 

AVERAGE  MISSEO  DISTANCE*  0.002  5  03 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y*  99.999870  AVERAGE  X*  79, 

AVERAGE  MISSEO  DISTANCE*  0.000130 


EMITTER  A  AVERAGE  ERROR 

YflQQ.OOQBQB.  _AVERAG 

AVERAG 


Mli?i5^8?l?iNCE*^'^0?So0809 


80  < 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y*  70.000431  AVERAGE  X*119. 

AVERAGE  MISSEO  DISTANCE*  0.000436 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y*250.002332  AVERAGE  X*149, 

AVERAGE  MISSEO  DISTANCE*  0.002533 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y*  99.999909  AVERAGE  X*  79, 

AVERAGE  MISSED  DISTANCE*  0.000091 


000032 


999922 


999883 


999997 


000025 


99  9934 


99  9904 


999997 
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TABLE  C-2 

Simulation  2,  Run  2.  Gradient  area  equal  to  15  Kn. 
All  values  are  in  kilometers.  Emitter  D  has  system 
erjx)r  of  2.0  degrees. 


0N0X«  ~  40.000 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y«100. 0433  80  AVERAGE  X*  80.002240 
AVERAGE  MISSED  DISTANCE-  0.043  438 


EMITTER  8  AVERAGE  ERROR 

AVERAGE  V-  70.019390  AVERAGE  X-119 .998662 
AVERAGE  MISSED  DISTANCE-  0.C19436 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y-250. 127474  AVERAGE  X-149. 999134 
AVERAGE  MISSED  DISTANCE-  0.  127477 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y-102. 759927  AVERAGE  X-  81.703879 
AVERAGE  MISSED  DISTANCE-  3.243517 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y-100.040231  AVERAGE  X-  80.001522 
AVERAGE  MISSED  DISTANCE-  0.040260 


EMITTER  8  AVERAGE  ERROR 

AVERAGE  Y-  70.020793  AVERAGE  X-119 .997969 
AVERAGE  MISSED  DISTANCE-  0.020892 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y-250. 12646 8  AVERAGE  X-149. 99  9207 
AVERAGE  MISSED  DISTANCE-  0.126471 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y-  98.64956  6  AVERAGE  X-  79.368710 
AVERAGE  MISSED  DISTANCE-  1.490704 


TABLE  C-2 

Simulation  2,  Run  2  (continued) 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y-103.042220  AVERAGE  X«  83.002135 
AVERAGE  MISSED  DISTANCE*  0.042274 


EMITTER 
AV  ERAGE 
AV  ERAGE 


B  AVERAGE  ERROR 

Y»  70.020391  AVERAGE  X»119. 997753 
MISSED  DISTANCE*  0.020514 


EM  I TTER 
AV  ERAGE 
AVERAGE 


C  AVERAGE  ERROR 

Y*250. 130829  AVERAGE  X*  149. 99  8  735 
HISSED  DISTANCE*  0.130835 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y*102« 725045  AVERAGE  X*  80.692302 
AVERAGE  MISSED  DISTANCE*  2.811610 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y*130.042224  AVERAGE  X*  80.002115 
AVERAGE  MISSED  DISTANCE*  0.042277 


EMITTER 
AV  ERAGE 
AV  ERAGE 


8  AVERAGE  ERROR 

Y*  70.019693  AVERAGE  X*119.998251 
MISSED  DISTANCE*  0.019770 


EMITTER 
AVERAGE 
AV  ERAGE 


C  AVERAGE  ERROR 

Y*250.127309  AVERAGE  X*149. 999286 
MISSED  DISTANCE*  0.127311 


EMITTER  D  AVERAGE  ERROR 

AVERAGE  Y*102.410316  AVERAGE  X*  79.863946 
AVERAGE  MISSED  DISTANCE*  2.414153 


TABLE  C-3 


Simulation  2,  Run  3.  Gradient  area  equal  to  15  Km. 
All  values  are  in  kilometers.  Emitter  D  has  system 
error  of  2.0  degrees. 


0N0X«  -  80.000 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  YalOO.086917  AVERAGE  X»  80.004486 
AVERAGE  NISSEO  DISTANCE*  0.C87032 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y«  70.038854  AVERAGE  X«119  .997338 
AVERAGE  NISSEO  DISTANCE*  0.038945 


EMITTER  C  AVERAGE  ERROR 
AVERAGE  Y*250.2560|3  AVERAGE 
AVERAGE  MISSED  DISTANCE^ 


__  X*  149. 99 8342 
0.256  0  29 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y»102.75992  7  AVERAGE  X*  81.70  3  8  79 
AVERAGE  MISSED  DISTANCE*  3.243517 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y«100. 080602  AVERAGE  X*  80.003  047 
AVERAGE  MISSED  DISTANCE*  0.080659 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y«  70.041659  AVERAGE  X*  119.995953 
AVERAGE  NISSEO  DISTANCE*  0.041855 


emitter  c  average  error 

AVERAGE  Y*250. 253971  AVERAGE  X*149. 998479 
AVERAGE  MISSED  DISTANCE*  0.253  9  75 


EMITTER 
AVERAGE 
AV  ERAGE 


0  AVERAGE  ERROR 

Y*  98.649566  AVERAGE  X*  79.368710 
MISSED  DISTANCE*  1.490704 


TABLE  C-3 

Simulation  2,  Run  3  (continued). 


AVERAGE  Y«lOoJoi4615^°AVERAGE  X*  80.004267 
AVERAGE  HISSED  DISTANCE*  0.084722 


ISiKil  X.U9.995525 

avIrasI  missed  01STANCE>  0.041109 


ISiSill  54^l?^ll4ir Average  x.i49 .997556 
aSerAGE  MitsEO  DISTANCE.  0.262498 


rr 


miSIli  ?4SI??iioir5vERAGE  X.  80.692302 
AVERAGE  MISSED  DISTANCE-  2.811610 


f-Kl^SilsirSvERAGE  X.  80.004232 
AVERAGE  MISSED  DISTANCE-  0.084696 


ISiSIlS  ?.*?l5Sfl4f°5vERACE  X.U9. 996522 
AvIragI  MISSeS  distance-  0.039621 


I! ISIfE  5-f5!J?fi8ir“|vERAGEX.  149 .99  8652 
average  HISSED  DISTANCE*  0.255845 


X*  79.863946 

AVERAGE  HISSED  DISTANCE*  2.414153 
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TABLE  C-4 


Simulation  2,  Run  4.  Gradient  area  equal  to  15  Km. 
All  values  are  in  kilometers.  Emitter  D  has  system 
error  of  2.0  degrees. 


0N0X«  -120.000 


EHITTER  A  AVERAGE  ERROR 

AVERAGE  Y-100.13048  5  AVERAGE  X»  80.006726 
AVERAGE  MISSED  DISTANCE*  0.130658 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y*  70.058320  AVERAGE  X*119. 996020 
AVERAGE  MISSED  DISTANCE-  0.058456 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y-250. 384719  AVERAGE  X-149. 997550 
AVERAGE  MISSEO  DISTANCE-  0.384727 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y-102. 759927  AVERAGE  X-  81.703879 
AVERAGE  MISSEO  DISTANCE-  3.243517 


EM  I TTER 
AVERAGE 
AV  ERAGE 


A  AVERAGE  ERROR 

Y-100. 121006  AVERAGE  X-  80.004569 
MISSEO  DISTANCE-  0.121092 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y-  70.062523  AVERAGE  X-U9. 993945 
AVERAGE  MISSEO  DISTANCE-  0.062  8  15 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y-250. 381620  AVERAGE  X-149. 997752 
AVERAGE  MISSEO  DISTANCE-  0.381627 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y«  98.649566  AVERAGE  X-  79.368710 
AVERAGE  MISSEO  DISTANCE-  1.A90704 
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TABLE  C-4 

Simulation  2,  Run  4  (cc *tinued) 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Vs lOO*  127041  AVERAGE  X*  80.006394 
AVERAGE  MISSED  DISTANCE*  0.127202 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Ys  70.061337  AVERAGE  Xsll9. 993306 
AVERAGE  MISSED  DISTANCE*  0.061702 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  ¥*250.394293  AVERAGE  X«149 .996377 
AVERAGE  MISSED  DISTANCE*  0.394310 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  ¥*102.725045  AVERAGE  X*  80.692302 
AVERAGE  MISSED  DISTANCE*  2.811610 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  ¥*100.126989  AVERAGE  X*  80.006345 
AVERAGE  MISSED  DISTANCE*  0.127147 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y*  70.059245  AVERAGE  X*  119  .994801 
AVERAGE  MISSED  DISTANCE*  0.059473 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  ¥*250.3  8452  0  AVERAGE  X*  149. 99 8019 
AVERAGE  MISSED  DISTANCE*  0.384525 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  ¥*102.410316  AVERAGE  X*  79.863946 
AVERAGE  MISSED  DISTANCE*  2.414153 
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TABLE  C-5 


Simulation  2,  Run  5.  Gradient  area  equal  to  15  Km. 
All  values  are  in  kilometers.  Emitter  D  has  system 
error  of  2.0  degrees. 


0N0X«  -150.000 


EMITTER 
AV  ERA6E 
AVERAGE 


A  AVERAGE  ERROR 

Y«100. 163182  AVERAGE  X«  80. 

MISSED  DISTANCE*  0.163398 


00  8403 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y«  70.072923  AVERAGE  X«119. 995035 
AVERAGE  MISSED  DISTANCE*  0.073092 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y*250. 481339  AVERAGE  X*149. 996957 
AVERAGE  MISSED  DISTANCE*  0.481348 


EMITTER 

AVERAGE 

AVERAGE 


0  AVERAGE  ERROR 
Y*r02.75992  7  AVERAG 
MISSED  DISTANCE*  3. 


E  X*  81 .703879 
243517 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y«100. 151332  AVERAGE  X*  80.005709 
AVERAGE  MISSED  DISTANCE*  0.151440 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y*  70.078171  AVERAGE  X*119 .992447 
AVERAGE  MISSED  DISTANCE*  0.078535 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y*250. 477453  AVERAGE  X*  149  .  99  7206 
AVERAGE  MISSED  DISTANCE*  0.477462 


EMITTER 

AVERAGE 

AVERAGE 


0  AVERAGE  ERROR 

Y*  98.649566  AVERAGE  X*  79.368710 
MISSED  DISTANCE*  1.490704 
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TABLE  C-5 

Simulation  2,  Run  5  (continued) 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  y>130.15d883  AVERAGE  X«  80.007987 
AVERAGE  MISSED  DISTANCE*  0.159084 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y*  70.076691  AVERAGE  X*  119. 991647 
AVERAGE  MISSED  DISTANCE*  0.077144 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  ¥*250.493245  AVERAGE  X*149. 995494 
AVERAGE  MISSED  DISTANCE*  0.493266 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  ¥*102.725045  AVERAGE  X*  80.692302 
AVERAGE  MISSED  DISTANCE*  2.811610 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  ¥*100.158809  AVERAGE  X*  80.00  7927 
AVERAGE  MISSED  DISTANCE*  0.159007 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y*  70.074079  AVERAGE  X*il9. 993514 
AVERAGE  MISSED  DISTANCE*  0.074363 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y*250. 481126  AVERAGE  X*149. 997545 
AVERAGE  MISSED  DISTANCE*  0.481133 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  ¥*102.410316  AVERAGE  X*  79.863946 
AVERAGE  MISSED  DISTANCE*  2.A14153 
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TABLE  C-6 

Simulation  2,  Run  6.  Gradient  area  equal  to  30  Km. 
All  values  are  in  kilometers.  Emitter  D  has  system 
error  of  2.0  degrees. 


0N0X«  -  80.000 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y«100.146546  AVERAGE  X«  80.007539 
AVERAGE  HISSEO  DISTANCE*  0.146740 


EMITTER  6  AVERAGE  ERROR 

AVERAGE  Y«  70.044350  AVERAGE  X*119  .996988 
AVERAGE  MISSED  DISTANCE*  0.044452 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y*250. 494021  AVERAGE  X* 149. 996880 
AVERAGE  MISSED  DISTANCE*  0.494030 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y*102. 759927  AVERAGE  X*  81.703879 
AVERAGE  MISSED  DISTANCE*  3.2435  17 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y*100.135913  AVERAGE  X*  80.005123 
AVERAGE  MISSED  DISTANCE*  0.136009 


EMITTER  6  AVERAGE  ERROR 

AVERAGE  Y«  70.047533  AVERAGE  X*119.99  5420 
AVERAGE  MISSED  DISTANCE*  0.047753 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y*250.490033  AVERAGE  X*149.99  7136 
AVERAGE  MISSED  DISTANCE*  0.490042 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y«  98.649566  AVERAGE  X*  79.368710 
AVERAGE  MISSED  DISTANCE*  1.A90704 
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TABLE  C-6 

Simulation  2,  Run  6  (continued). 


EMITTER 

AVERAGE 

AVERAGE 


A  AVERAGE  ERROR 
Y«I00.142684  AVERAG 
MISSED  DISTANCE*  0. 


E  X*  80.007169 
142864 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y«  70.046624  AVERAGE  X *119. 994932 
AVERAGE  MISSED  DISTANCE*  0.046898 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y*250.506230  AVERAGE  X*I49  .995379 
AVERAGE  MISSED  DISTANCE*  0.506251 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y*10 2. 725045  AVERAGE  X*  80 

AVERAGE  MISSED  DISTANCE*  2.811610 


692302 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y*100. 142623  AVERAGE  X*  80.007113 
AVERAGE  MISSED  DISTANCE*  0.142801 


EMITTER  8  AVERAGE  ERROR 

AVERAGE  Y*  70.045045  AVERAGE  X*  119. 996062 
AVERAGE  MISSED  DISTANCE*  0.045217 


EMITTER 

AVERAGE 

AVERAGE 


C  AVERAGE  ERROR 

Y*250. 493808  AVERAGE  X»149. 997484 
MISSED  DISTANCE*  0.493814 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y*102. 410316  AVERAGE  X*  79.86  3  946 
AVERAGE  MISSED  DISTANCE*  2.414153 
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TABLE  C-7 


Simulation  2,  I^un  7,  Gradient  area  equal  to  3  0  Km. 
All  values  are  in  kilometers.  Emitter  D  has  system 
error  of  2.0  degrees. 


DNOX*  -120.000 


EMITTER 
AV  ERAGE 
AV  ERAGE 


A  ERRriR 

Y>ia0.220026  AVERAGE  X«  80.011285 
MISSED  OlSTANCEa  0.220315 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y>  70.06652  8  AVERAGE  X«119. 995517 
AVERAGE  MISSED  DISTANCE^  0.066679 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y«25a.742346  AVERAGE  X«149. 995360 
AVERAGE  MISSED  DISTANCE*  0.742361 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y»102. 759927  AVERAGE  X*  81.703879 
AVERAGE  MISSED  DISTANCE*  3.243517 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y*100.204084  AVERAGE  X*  80.007670 
AVERAGE  MISSED  DISTANCE*  0.204228 


EM ITTER 
AVERAGE 
AVERAGE 


B  AVERAGE  ERROR 

Y*  70.071281  AVERAGE  X*119  .993181 
MISSED  DISTANCE*  0.  C7I606 


EMITTER 

AVERAGE 

AVERAGE 


C  AVERAGE  ERROR 

Y*250. 736341  AVERAGE  X*149. 995739 
MISSED  DISTANCE*  0.  736353 


EMITTER 

AVERAGE 

AVERAGE 


D  AVERAGE  ERROR 

Y*  98.649566  AVERAGE  X*  79.368710 
MISSED  DISTANCE*  1.490704 


TABLE  C-7 

Simulation  2,  Run  7  (continued) 


ENITTER  A  AVERAGE  ERROR 

AVERAGE  Y«130.214249  AVERAGE  X«  80.010728 
AVERAGE  NISSEO  DISTANCE*  0.214518 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y*  70.069917  AVERAGE  X*  119. 992449 
AVERAGE  MISSED  OIST  .NCE«  0.070324 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y*250« 760544  AVERAGE  X*  149.993117 
AVERAGE  MISSED  DISTANCE*  0.760575 


ENITTER  0  AVERAGE  ERROR 

AVERAGE  Y*102. 725045  AVERAGE  X*  80.69  2302 
AVERAGE  MISSED  DISTANCE*  2.811610 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y*IOO. 214143  AVERAGE  X*  80.010646 
AVERAGE  MISSED  DISTANCE*  0.214408 


ENITTER  B  AVERAGE  ERROR 

AVERAGE  Y«  70.06756  8  AVERAGE  X*  119. 99 4138 
AVERAGE  MISSED  DISTANCE*  0.067822 


EMITTER  C  AVERAGE  ERROR 

AVERAGE  Y*250«742099  AVERAGE  X*149. 996269 
AVERAGE  MISSED  DISTANCE*  0.742109 


EMITTER  0  AVERAGE  ERROR 

AVERAGE  Y*102. 410316  AVERAGE  X*  79.863946 
AVERAGE  MISSED  DISTANCE*  2.A14153 


TABLE  C-8 

Simulation  2,  Run  8.  Gradient  area  of  30  Km.  System 
induced  error  of  1.5  degrees  for  all  LOB’s.  Emitters 
A,  B  and  C  have  refracted  LOB’s  correlated  with 
straight  LOB's.  All  values  are  in  kilometers. 


ONOXs  -100.000 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  Y»  96.136822  AVERAGE  X»  80 .311562 
AVERAGE  MISSED  DISTANCE*  3.875721 


AvIrAGE  Y-*7oJo5l4o*^°AVERAG e  x- 119. 996250 
AVERAGE  MISSEO  DISTANCE-  0.055567 


lUiflil  5.t5!!S!|l!f‘’!ivERAG6  X.14,.996120 
AVERAGE  MISSEO  DISTANCE-  0.618124 


EMITTER  0  AVERAGE  ERROR  ^ 

AVERAGE  Y-101.991913  AVERAGE  X*^81. 272969 
AVERAGE  MISSED  DISTANCE-  2.363930 


X.  80.  1E7291 

AVERAGE  MISSED  DISTANCE-  1.489362 


?.*^i?Sll4SI”Sve»AG6  X.119.994296 
AVERAGE  NISSEO  DISTANCE.  a.C59683 


HJiJUl  f.!!!S?J?ll!2'“!!vERAGEX.l49.996437 
AVERAGE  MISSEO  DISTANCE-  0.613126 


ICIHil  ?.fX!?*7l!lirLERAGE  X.  79.514574 
AVERAGE  MISSEO  DISTANCE-  0.  677645 
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TABLE  C-8 

Simul^ition  2,  Run  8  (continued) 


EMITTER  A  AVERAGE  ERROR 

AVERAGE  98.926232  AVERAGE  X>  80. 

AVERAGE  MISSED  DISTANCE*  1.315779 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  y»  70.058276  AVERAGE^X«119  « 

AVERAGE  MISSED  DISTANCE*  0.058617 


EMITTER  C  AVERAGE  ERROR^,^^^^^  ^ 
AVERAGE  7*250.633317  ^ AVERAGE  X* 149 < 
AVERAGE  MISSED  DISTANCE*  0.633343 


EMITTER  0  AVERAGE  ERROR. ^ 
AVERAGE  Y*10 1.94669 9  AVERAGE  X*^80i 
AVERAGE  MISSED  DISTANCE*  2.013887 


EMITTER  A  AVERAGE  ERROR... 

AVERAGE  Y«101.203482  AVERAG E  X*  p  . 
AVERAGE  MISSED  DISTANCE*  1.247070 


EMITTER  B  AVERAGE  ERROR 

AVERAGE  Y*  70.056309  AVERAGE.^X«119  , 

AVERAGE  MISSED  DISTANCE*  0.056522 


EMITTER  C  AVERAGE  ERROR  ...... 

AVERAGE  7*250.617883  ^AVERAGE  X*  149 
AVERAGE  MISSED  DISTANCE*  0.617891 


EMITTER  0  AVERAGE  ERROR...... ^  _ 

AVERAGE  7*131.649868  AVERAGE  X*  79 
AVERAGE  MISSED  DISTANCE*  1.653102 


760453 


993687 


994248 


515853 


673178 


995098 


996876 


896641 
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